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Abstract  (Short  version) 


As  the  only  low-latitude  connection  between  ocean  basins,  the  Indonesian 
Throughflow  allows  the  direct  transmission  of  heat  and  salinity  between  the  Pacific 
and  Indian  Oceans.  The  Mg/Ca  and  5180  of  calcite  of  Globigerinoides  ruber  ( G .  ruber ) 
were  used  to  estimate  the  sea  surface  temperature  (SST)  and  5180  of  water,  an 
indicator  of  hydrologic  conditions,  over  the  past  20,000  years.  I  also  attempted  to 
estimate  thermocline  structure  using  PuIIeniatina  obliquiloculata,  but  the  Mg/Ca  and 
5180  of  calcite  data  yield  conflicting  interpretations,  indicating  further  work  on  this 
proxy  is  required.  The  G.  ruber  Mg/Ca  results  suggest  that  the  SST  of  the  outflow 
passages  were  influenced  by  high  latitude  Southern  Hemisphere  temperature.  At 
approximately  10,000  years  before  present,  there  was  a  warming  in  the  Makassar 
Strait.  This  local  warming  was  coincident  with  the  flooding  of  the  Sunda  Shelf,  which 
opened  a  connection  between  the  South  China  Sea  and  the  Indonesian  Throughflow. 
Regional  5180  of  seawater  reconstructions  suggest  that  the  mean  position  of  the 
Intertropical  Convergence  Zone  (ITCZ)  was  approximately  the  same  as  modern  at 
the  last  glacial  maximum  and  was  displaced  to  the  south  during  the  Younger  Dryas 
and  Heinrich  Stadial  1,  suggesting  the  ITCZ  responds  to  changes  in  the 
interhemispheric  temperature  gradient. 
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Abstract 

As  the  only  low-latitude  connection  between  ocean  basins,  the  Indonesian 
Throughflow  allows  the  direct  transmission  of  heat  and  salinity  between  the  Pacific 
and  Indian  Oceans.  Despite  its  potential  importance,  the  role  of  the  Indonesian 
Throughflow  in  global  ocean  circulation  and  regional  climate  is  still  not  clear  due  to 
sparse  measurements  and  the  relative  difficulty  of  modeling  the  region.  The  Mg/Ca 
and  5180  of  calcite  of  the  calcitic  planktic  foraminifera  Globigerinoides  ruber  ( G . 
ruber )  were  used  to  estimate  the  sea  surface  temperature  and  5180  of  water,  an 
indicator  of  hydrologic  conditions,  over  the  past  20,000  years.  I  also  attempted  to 
estimate  thermocline  structure  using  the  foraminifera,  Pulleniatma  obliquiloculata, 
but  the  Mg/Ca  and  5180  of  calcite  data  yield  conflicting  interpretations,  indicating 
further  work  on  this  proxy  is  required.  The  G.  ruber  Mg/Ca  results  suggest  that  the 
sea  surface  temperature  of  the  outflow  passages  was  influenced  by  high  latitude 
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Southern  Hemisphere  temperature.  This  connection  is  likely  via  intermediate 
waters  that  upwell  in  the  Banda  Sea.  At  approximately  10,000  years  before  present, 
there  was  a  warming  in  the  Makassar  Strait.  This  local  warming  was  coincident  with 
the  flooding  of  the  Sunda  Shelf,  which  opened  a  connection  between  the  South  China 
Sea  and  the  Indonesian  Throughflow.  Regional  S180  of  seawater  reconstructions 
show  that  during  the  last  glacial  maximum  the  5180  of  seawater  pattern  was  very 
similar  to  modern,  but  there  were  relatively  enriched  values  over  the  equatorial 
IndoPacific  during  high  latitude  Northern  Hemisphere  cold  events  (Heinrich  Stadial 
1  and  the  Younger  Dryas).  From  these  results  we  postulate  that  the  mean  position  of 
the  Intertropical  Convergence  Zone  was  approximately  the  same  as  modern  at  the 
last  glacial  maximum  and  was  likely  displaced  to  the  south  during  the  Younger 
Dryas  and  Heinrich  Stadial  1,  suggesting  the  Intertropical  Convergence  Zone 
primarily  responds  to  changes  in  the  interhemispheric  temperature  gradient.  These 
results  shed  light  on  the  primary  controls  of  the  temperature  and  hydrology  of 
Indonesian  Throughflow  region. 
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Chapter  1 


Introduction 

Approximately  10-15  Sv  (1  Sv  =  a  million  cubic  meters  per  second)  of  water 
makes  its  way  from  the  Pacific  Ocean  into  the  Indian  Ocean,  meandering  through  a 
series  of  straits  known  as  the  Indonesian  Throughflow  [Gordon,  2005;  Sprintall  et 
al,  2009].  The  water’s  path  is  not  straightforward;  there  are  rich  seasonal  and 
interannual  cycles,  topographic  barriers,  and  vigorous  tidal  and  wind  driven  mixing 
[Ffield  and  Gordon,  1996;  Gordon  etal.,  2010;  Koch-Larrouy  etal.,  2008].  There  are 
three  main  entrances  into  the  Indonesian  Throughflow:  the  Makassar  Strait,  the  ]ava 
Sea,  and  the  Eastern  Passages.  Of  these  the  Makassar  Strait  dominates 
volumetrically,  supplying  water  from  the  subtropical  North  Pacific  [Gordon  etal., 
2008].  The  ]ava  Sea  adds  freshwater  from  the  South  China  Sea  [Qu  etal,  2006]  and 
the  Eastern  Passages  contain  the  only  deep  entrance  (the  Lifamatola  Strait,  sill 
depth  2,000m)  [Talley  and  Sprintall,  2005].  There  are  also  three  outflow  passages: 
the  Timor  Strait,  the  Ombai  Strait,  and  the  Lombok  Strait  (Fig  1).  Each  outflow 
passage  has  its  own  characteristics  as  well.  The  smallest,  the  Lombok  Strait,  is  the 
most  influenced  by  the  seasonal  reversal  of  the  winds  (the  monsoons).  During  the 
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boreal  winter  months,  the  Lombok  Strait  occasionally  loses  its  status  as  an  outflow 
passage;  its  flow  is  reversed  by  the  winds,  particularly,  at  the  surface.  The  Ombai 
Strait  also  has  occasional  seasonal  flow  reversals  in  the  surface,  and  the  largest 
passage,  the  Timor  Strait  has  very  little  seasonal  variability  [Sprintall  et  al,  2009]. 

Though  volumetrically  the  most  important,  and  relatively  well  studied,  the 
Makassar  Strait  is  still  poorly  understood.  Initial  measurements  indicated  the  flow 
of  the  surface  layer  reversed  with  the  seasons,  with  a  strong  and  constant 
thermocline  flow.  This  reduced  surface  flow  was  attributed  to  a  plume  of  freshwater 
from  the  Java  Sea  blocking  flow  during  the  boreal  winter  [Gordon  et  al,  2003], 
though  it  could  also  be  explained  by  the  Ekman  transport  caused  by  the  monsoon 
winds  [Sprintall  and  Liu,  2005].  Later  work  showed  that  the  surface  flow  is  inhibited 
throughout  the  year,  and  the  maximum  flow  is  always  in  the  thermocline  [Gordon  et 
al,  2008].  Neither  of  the  proposed  mechanisms  explains  a  year  round  suppression 
of  surface  flow.  The  reason  for  this  year  round  suppression  of  the  surface  flow  is  still 
not  clear.  Here,  we  propose  an  alternative  mechanism:  the  persistently  higher  sea 
level  of  the  southern  Java  Sea  compared  to  the  Flores  Sea  creates  a  pressure 
gradient.  A  surface  northward  geostrophic  flow  is  established,  inhibiting  surface 
flow.  The  thermocline  flow  is  unaffected  since  the  sill  depth  between  the  Java  Sea 
and  the  Indonesian  Throughflow  is  only  40m. 

Oceanography  (Chapter  2): 
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What  can  our  climate  reconstructions  add  to  our  understanding  of  the 
Indonesian  Throughflow?  Quite  a  lot,  it  turns  out.  We  can  test  some  of  the  emerging 
hypotheses  about  the  Indonesian  Throughflow.  We  can  start  with  a  very  basic 
question:  does  the  water  exiting  the  Indonesian  Throughflow  resemble  the  water 
going  in?  Not  at  all.  The  water  that  enters  the  Indonesian  Throughflow  is  warm  and 
salty.  What  comes  out  is  relatively  cool  and  fresh  (Fig  2).  The  fundamental  nature  of 
this  great  transformation  is  the  focus  of  work  in  modern  oceanography  [Fang  et  al., 
2010;  Gordon  etal,  2003;  Koch-Larrouy  etal,  2008],  How  can  the  paleo  record  help 
us  understand  this  issue?  A  current  hypothesis  for  explaining  the  cool,  fresh  nature 
of  the  Indonesian  Throughflow  is  that  water  from  the  Java  Sea  blocks  the  warm 
surface  flow  through  the  Makassar  Strait.  The  flow  of  the  cooler  subsurface  water  is 
unaffected  and  therefore  this  cooler  water  dominates  [Gordon  et  al.,  2003].  In  a 
stroke  of  geological  good  luck,  the  Java  Sea  is  so  shallow  that  20,000  years  ago  when 
the  great  continental  ice  sheets  robbed  the  oceans  of  water  and  sea  level  was  120 
meters  lower,  the  entire  Java  Sea  was  exposed  shelf  area.  Thus  by  reconstructing 
oceanographic  conditions  as  sea  levels  rose  to  modern,  we  can  examine  what,  if  any, 
effect  the  Java  Sea  water  had  on  the  Indonesian  Throughflow. 

Our  long  records  can  also  give  us  insight  into  the  major  controls  on  the 
temperature  of  the  Indonesian  Throughflow.  Here,  we  find  something  unexpected; 
the  opening  of  the  Java  Sea  when  sea  levels  rose  did  not  appear  to  affect  the  sea 
surface  temperatures  of  the  outflow  passages.  Instead,  over  the  past  20,000  years 
the  sea  surface  temperatures  of  the  Indonesian  Throughflow  outflow  paralleled  high 
latitude  Southern  Hemisphere  temperature  variations  recorded  in  ice  cores.  This  is 
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not  a  temperature  pattern  we  see  in  the  inflow  area  to  the  Indonesian  Throughflow. 
How  does  this  temperature  pattern  get  to  the  Indonesian  Throughflow?  We  suggest 
that  there  is  more  deep  water  entering  through  a  poorly  studied  passage,  the 
Lifamatola  Strait,  than  is  currently  estimated.  This  water  is  then  vigorously  mixed  in 
the  Banda  Sea,  which  brings  deep  water  to  the  surface.  It  is  likely  that  this  mixing  is 
the  primary  source  of  the  cool  the  Indonesian  Throughflow  waters  in  the  modern,  as 
our  records  suggest  the  Java  Sea  inflow  does  not  change  the  sea  surface  temperature 
of  the  outflow  passages.  Further  observational  work  in  the  Indonesian  Throughflow 
region  will  allow  us  to  test  this  hypothesis. 

Hydrology  (Chapter  3): 

The  primary  controls  on  the  sea  surface  temperature  of  the  Indonesian 
Throughflow  turned  out  to  be  much  more  complicated  than  expected  and  highly 
dependent  on  the  unique  topography  of  the  region.  As  is  turns  out,  the  major 
changes  in  hydrology  are  much  more  straightforward.  Along  the  equator,  winds 
from  the  north  and  south  meet  over  warm  water.  At  this  point  of  convergence,  air 
rises  and  heavy  precipitation  occurs.  This  area  of  convection  is  the  Intertropical 
Tropical  Convergence  Zone  (ITCZ).  The  ITCZ  migrates  with  the  seasons  and  gets 
tugged  towards  the  warmer  hemisphere.  Due  to  the  current  continental 
configuration,  the  northern  hemisphere  is  slightly  warmer,  and  the  position  of  the 
ITCZ  is  biased  to  the  north  [Koutavas  and  Lynch-Stieglitz,  2005].  What  would  happen 
if  the  interhemispheric  temperature  gradient  were  reduced?  Temperature 
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compilations  show  that  during  times  of  rapid  cooling  in  the  northern  hemisphere 
(for  example  during  Heinrich  Event  1  and  the  Younger  Dryas),  the  southern 
hemisphere  warmed,  decreasing  the  interhemispheric  temperature  gradient 
[Shakun  and  Carlson,  2010].  The  expectation  is  that  the  ITCZ  would  lose  some  of  its 
northward  bias.  Our  hydrologic  proxies  show  exactly  that.  Our  data  suggest  that, 
despite  the  complexity  of  the  Indonesian  Throughflow,  the  ITCZ  responded 
uniformly  to  these  high  latitude  northern  hemisphere  cooling  events  and  its  mean 
position  shifted  southward. 

Proxies  (Chapter  4): 

In  order  to  reconstruct  temperature  and  hydrology  on  timescales  longer  than 
the  instrumental  record,  we  have  turned  to  proxies.  Proxies  can  give  us  estimates  of 
past  oceanographic  conditions.  The  data  in  this  thesis  is  based  on  the  geochemistry 
of  planktic  foraminifera.  Foraminifera  are  single  celled  organisms  that  build  shells, 
specifically  of  calcium  carbonate  (CaCOs),  and  we  utilize  the  chemical  variation  in 
these  shells.  The  amount  of  180  compared  to  160  (or  5180)  that  is  incorporated  into  a 
foraminifera’s  shell  depends  on  both  the  temperature  and  the  5180  of  the  water  in 
which  the  shell  grows  [Emiliani,  1970].  The  5180  of  the  water  can  tell  us  about 
hydrologic  processes,  but  we  need  an  independent  temperature  estimate  to  isolate 
the  5180  of  seawater  component.  Here  we  take  advantage  of  the  observation  that 
foraminifera  growing  at  higher  temperatures  incorporate  more  Mg  into  their  shells. 
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Thus,  the  Mg/Ca  ratio  in  a  shell  can  be  used  estimate  past  temperatures  [Rosenthal 
etal,  1997]. 

We  investigate  two  types  of  foraminifera.  Globigerinoides  ruber  ( G .  ruber)  has 
photosynthetic  symbionts  and  is  primarily  found  in  the  mixed  layer.  PuIIeniatina 
obliquiloculata  (P.  obliquiloculata)  lives  in  the  upper  thermocline  at  ~80-100m.  By 
measuring  both  species,  the  hope  was  to  reconstruct  thermocline  structure.  This 
was  a  very  exciting  prospect.  Firstly,  early  observational  results  from  the  Makassar 
Strait  showed  that  thermocline  temperature,  but  not  surface  temperature  was 
correlated  to  Makassar  Strait  transport  [Ffield  etal.,  2000],  Secondly,  temperature 
anomalies  associated  with  the  Indian  Ocean  Dipole  are  often  larger  in  the 
subsurface  than  in  the  surface  [Shinoda  et  al,  2004].  Finally,  modeling  simulations  of 
the  Last  Glacial  Maximum  suggest  that  thermocline  temperatures  in  the  Western 
Equatorial  Pacific  may  be  the  best  diagnostic  for  past  changes  in  the  Walker 
Circulation  [DiNezio  et  al.,  2011].  Thus,  the  depth  habitat  of  P.  obliquiloculata  makes 
this  foraminifera  well  suited  to  address  several  important  climatic  questions. 

While  our  5180Caidte  and  Mg/Ca-based  temperature  results  from  G.  ruber  are 
self-consistent,  and  can  be  reconciled  with  modern  conditions  and  other  climate 
reconstructions  (and  are  the  basis  for  the  results  described  in  Chapters  2  and  3),  our 
P.  obliquiloculata  results  are  difficult  to  interpret.  We  use  the  difference  between  the 
Mg/Ca-based  temperatures  of  G.  ruber  -  P.  obliquiloculata  to  infer  the  thermocline 
structure.  Using  the  5180Caidte  of  G.  ruber  and  P.  obliquiloculata  we  make  a  similar 
estimate.  In  the  coretop  data,  these  two  methods  yield  qualitatively  similar  results. 

In  the  downcore  records,  they  predict  a  different  sign  of  change  in  the  thermocline 
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structure,  making  an  oceanographic  interpretation  difficult.  While  previous  authors 
have  suggested  that  the  unique  oceanographic  conditions  at  their  study  sites  cause 
this  inconsistency  [Steinke  et  al.,  2010;  Xu  et  ah,  2008],  we  suggest  that  the  difficulty 
of  using  the  combination  of  G.  ruber  and  P.  obliquiloculata  to  reconstruct 
thermocline  conditions  is  pervasive  throughout  the  IndoPacific  Warm  Pool.  We  urge 
caution  in  using  the  geochemistry  of  P.  obliquiloculata  to  reconstruct  water  column 
structure  and  suggest  more  work  is  required  before  this  proxy  can  be  reliably  used. 

Summary 

Not  only  is  the  Indonesian  Throughflow  at  an  oceanic  crossroads,  this  thesis 
demonstrates  it  is  also  a  climatic  balancing  point  between  the  Northern  and 
Southern  Hemisphere.  High  latitude  Southern  Hemisphere  temperature  variation 
appears  to  influence  the  sea  surface  temperature  of  the  outflow  passages.  On  the 
other  hand,  the  hydrology  of  the  entire  region  responds  in  concert  with  changes  in 
the  interhemispheric  temperature  gradient.  This  suggests  that  the  geometry  of  the 
Indonesian  Seas  has  an  influence  on  local  sea  surface  temperatures,  but  the 
hydrologic  changes  are  primarily  influenced  by  extratropical  forcings. 
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Figure  1:  Map  of  Indonesian  Throughflow.  Water  flows  from  the  Pacific  into  the 
Indian  Ocean  through  the  Indonesian  Throughflow.  Water  from  the  three  major 
inflow  passages  (the  Makassar  Strait,  the  Java  Sea,  and  the  Eastern  Passages)  is 
mixed  within  the  Banda  Sea  and  flows  out  via  the  Timor  Sea,  the  Ombai  Strait,  and 
the  Lombok  Strait.  The  circulation  schematic  is  based  on  [Fang  et  al.,  2010;  Gordon 
etal,  2010;  Tomczak  and  Godfrey,  1994]. 
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Figure  2:  Hydrographic  data  from  within  the  Indonesian  Seas.  Data  are  from 
CTD  casts  from  two  cruises,  BJ8-03  and  S0184,  both  of  which  took  place  during  the 
boreal  summer.  The  inflow  from  both  the  North  Pacific  (red)  and  the  South  Pacific 
(teal)  have  subsurface  salinity  maxima.  Vigorous  mixing  and  freshwater  influx  (both 
from  the  Java  Sea  and  local  precipitation)  erode  these  salinity  maxima  and  water 
exits  the  Indonesian  Throughflow  fresher  than  it  enters. 
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Chapter  2 


Southern  Hemisphere  Influence  on  the  Indonesian  Throughflow 


Abstract 

In  the  modern  ocean,  vigorous  vertical  mixing  in  the  Banda  Sea  combines 
North  Pacific  and  South  China  Sea  surface  and  thermocline  waters  with 
intermediate  water  from  the  southern  hemisphere  to  create  a  distinctive  Indonesian 
Throughflow  water  mass  that  flows  from  the  tropical  Pacific  into  the  Indian  Ocean 
[Gordon  et  ah,  2003;  Koch-Larrouy  et  al,  2008].  The  Indonesian  Throughflow 
contributes  to  the  warm  water  return  path  of  the  meridional  overturning 
circulation,  but  the  controls  on  and  the  influence  of  the  Indonesian  Throughflow  on 
deglacial  and  Holocene  sea  surface  temperature  evolution  are  unknown.  Here,  we 
present  Mg/Ca-based  sea  surface  temperature  reconstructions  from  the  Makassar 
Strait,  the  main  path  for  North  Pacific  water  to  enter  the  Indonesian  Throughflow, 
and  from  the  Savu  Sea,  where  Indonesian  Throughflow  exits  into  the  eastern  Indian 
Ocean.  We  show  that  during  the  last  deglaciation  surface  temperature  trends  in  both 
the  Makassar  Strait  and  the  Indonesian  Throughflow  outflow  resemble  high-latitude 
southern  hemisphere  temperature  trends.  In  contrast  to  the  high-latitude  southern 
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hemisphere,  the  Makassar  Strait  warmed  ~10,000  years  ago,  when  sea  level  rise 
permitted  the  inflow  of  South  China  Sea  water  across  the  previously  exposed  Sunda 
Shelf  [Linsley  et  al.,  2010],  which  apparently  isolated  the  Makassar  Strait  from 
southern  hemisphere  influence.  We  argue  that  sea  surface  temperatures  of  the 
outflow  passages  of  the  Indonesian  Throughflow  were  unaffected.  We  suggest  that  a 
sustained  oceanic  connection  between  the  southern  hemisphere  and  the  Indonesian 
Throughflow  through  the  mixing  in  the  Banda  Sea  provides  a  mechanism  for 
Southern  Source  Water  to  influence  tropical  Indian  Ocean  sea  surface  temperatures. 

Introduction: 

The  Indonesian  Throughflow  (ITF)  transfers  approximately  10-15  Sv  (1  Sv  = 
106  m3/s)  of  water  from  the  Pacific  to  the  Indian  Ocean  through  the  Indonesian  Seas, 
and  plays  an  important  role  in  the  global  oceanic  circulation.  Water  enters  the  ITF 
via  three  routes:  the  Sulawesi  Sea,  the  Java  Sea,  and  the  Eastern  Passages  (Figure  1) 
[Gordon  etal,  2003;  Koch-Larrouy  etal,  2008],  North  Pacific  surface  and 
thermocline  water  flows  through  the  Sulawesi  Sea  and  into  the  Makassar  Strait. 
From  the  Makassar  Strait,  some  flows  into  the  Indian  Ocean  via  the  Lombok  Strait, 
with  the  rest  flowing  to  the  Flores  Sea,  and  then  eastward  into  the  Banda  Sea 
[Gordon  et  al,  2003;  Koch-Larrouy  etal.,  2008].  Relatively  fresh  surface  water  from 
the  South  China  Sea  is  advected  across  the  shallow  Sunda  Shelf  and  through  the  Java 
Sea,  and  into  the  southern  Makassar  Strait  and  Flores  Sea  [Fang  etal.,  2010]. 

Approximately  2.5  Sv  of  Southern  Source  Waters  enter  the  ITF  through  the 
Eastern  Passages  between  Sulawesi  and  Irian  Jaya  [Van  Aken  etal.,  2009].  The 
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Southern  Source  Water  is  composed  of  both  of  Antarctic  (0  =  -1.3 1°C,  salinity  = 
34.62)  and  Subantarctic  water  (0  =  5.09°C,  salinity  =  34.25)  (potential  temperature 
and  salinity  estimates  are  from  Gebbie  and  Huybers  [2011]).  Water  from  the 
Eastern  Passages,  Makassar  Strait  and  Java  Sea  enter  the  Banda  Sea,  where  vigorous 
vertical  mixing,  both  tidal  and  wind  driven,  combines  the  water  masses,  mixing 
down  relatively  fresh  surface  water  throughout  the  column  [Koch-Larrouy  et  al, 
2008],  and  creating  characteristic  cool,  low-salinity  ITF  water  [Fang  etal,  2010; 
Gordon  etal.,  2003].  This  ITF  water  then  flows  into  the  Indian  Ocean  via  the  Ombai 
Strait  and  the  Timor  Sea. 

Though  a  relatively  small  net  volumetric  contribution  [Fang  et  al.,  2010],  the 
flow  from  the  South  China  Sea  plays  a  key  role  in  setting  the  character  of  ITF  water 
and  transport  [Gordon  etal.,  2003;  Koch-Larrouy  et  al.,  2008].  It  has  been  proposed 
that  the  advection  of  freshwater  from  the  South  China  Sea  (Figure  SI)  creates  a 
northward  pressure  gradient  within  the  Makassar  Strait,  increasing  the  residence 
time  of  surface  water  in  the  Strait,  resulting  in  thermocline-intensified  flow  [Gordon 
etal.,  2003],  We  suggest  an  alternate  explanation  is  required  since  the  reduced 
surface  flow  relative  to  the  thermocline  occurs  throughout  the  year  [Gordon  et  al., 
2008],  even  in  boreal  summer  when  the  surface  flow  is  towards  the  South  China  Sea. 
In  all  seasons,  the  sea  level  of  the  southern  Java  Sea  is  higher  than  the  Banda  Sea 
[Carton  and  Giese,  2008]  (Figure  Sl-3).  The  resulting  pressure  gradient  may  be 
associated  with  a  northward  geostrophic  flow  that  opposes  southward  transport 
through  the  Makassar  Strait.  Since  the  sill  between  the  Java  Sea  and  the  ITF  is  only 
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40m,  this  pressure  gradient  is  not  present  at  thermocline  depths,  explaining  the 
relatively  enhanced  subsurface  flow. 

Previous  studies  have  shed  light  on  the  role  of  the  South  China  Sea 
Throughflow  by  focusing  on  the  effect  of  the  flooding  of  the  shallow  Sunda  Shelf 
~10,000  years  before  present  (yr  BP).  This  flooding  opened  a  connection  for  South 
China  Sea  Throughflow  to  enter  the  Makassar  Strait,  freshening  the  ITF  relative  to 
the  Western  Pacific  Warm  Pool  [Linsley  etal.,  2010].  Additionally,  cooling  in  the 
thermocline  of  the  ITF  outflow  was  attributed  to  thermocline  intensified  flow  that 
resulted  from  this  fresh  South  China  Sea  inflow  [Xu  etal.,  2008].  Previous  work  has 
also  shown  that  the  deglacial  sea  surface  temperature  (SST)  changes  in  the 
Makassar  Strait  and  Antarctica  were  synchronous  [Visser  et  al,  2003].  The  link  was 
attributed  to  a  switch  between  an  El  Nino-like  climate  during  glacial  times  to  a  La 
Nina-like  climate  during  deglacial  times. 

Here,  we  evaluate  the  link  between  high  latitude  southern  hemisphere 
temperatures  and  the  ITF  in  greater  detail,  and  re-evaluate  the  influence  of  the 
South  China  Sea  on  SST.  We  present  core-top  and  downcore  Mg/Ca-based  SST 
estimates,  using  the  surface  mixed-layer  planktic  foraminifera  Globigerinoides  ruber 
[G.  ruber],  to  study  past  changes  in  the  SST  of  the  ITF  region. 

Materials  and  Methods 

Approximately,  40  G.  ruber  (white,  morphotype  sensu  stricto ;  212-300pm) 
were  used  for  Mg/Ca  analysis.  Samples  were  gently  crushed;  cleaned  using  a 
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method  that  included  reductive  followed  by  oxidative  steps;  and  run  on  an  Element 
XR ICP-MS  at  Rutgers  University  or  an  Element  2  ICP-MS  at  Woods  Hole 
Oceanographic  Institution. 

The  age  models  for  all  sediment  records  presented  here  are  based  on  the 
radiocarbon  ages  of  planktic  foraminifera  (Appendix  A).  All  radiocarbon  ages  were 
converted  to  calendar  ages  using  the  Marine09  calibration  [Reimer  et  al,  2009]  and 
no  local  reservoir  age  adjustment  (see  also  supplementary  discussion). 

Results: 

In  agreement  with  a  previous  study  [Arbuszewski  et  al,  2010],  and  in  contrast 
with  another  [Mathien-Blard  and  Bassinot,  2009],  our  new  core  top  data,  in 
conjunction  with  data  from  other  regions  [Arbuszewski  etal,  2010;  Dahl  and  Oppo, 
2006;  Ferguson  etal.,  2008;  McConnell  and  Thunell,  2005;  Mohtadi  etal.,  2011a; 
Regenberg  etal,  2006],  demonstrate  that  there  is  no  discernable  salinity  influence 
on  the  Mg/Ca-based  temperature  estimates  of  G.  ruber  at  low  salinities.  Arubszewski 
etal.  [2010]  suggest  this  threshold  is  35.  With  more  data  from  fresh  regions,  we 
suggest  the  threshold  is  36  (Figure  2  and  supplementary  discussion);  thus,  we  can 
interpret  the  Mg/Ca  data  strictly  in  terms  of  temperature  with  no  correction  for 
salinity  effects. 

Our  new  high-resolution  SST  record  from  the  Savu  Sea  (GeoB10069-3,  here 
69-3)  is  bathed  in  water  supplied  by  the  ITF  entering  the  Indian  Ocean  (Figure  1)  via 
the  Ombai  Strait.  In  addition  we  generated  a  high-resolution  SST  record  from  the 
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Flores  Sea  (BJ8-03  136GGC,  here  136GGC)  and  a  low-resolution  SST  record  from  the 
Makassar  Strait  (BJ8-03-23GGC,  here  23GGC). 

We  group  temperature  reconstructions  of  our  new  records  and  published 
records  into  three  regions:  the  Western  Pacific  Warm  Pool  [Rosenthal  et  al,  2003; 
Stott  et  al,  2007],  the  Southern  Makassar  Strait  [Linsley  et  al,  2010;  Visser  et  al, 
2003],  and  the  outflow  of  the  ITF  [Levi  et  al,  2007;  Xu  et  al,  2008]  (Figure  1,  Table 
SI),  and  reconstruct  regional  averages  (Figure  3  and  4).  We  include  136GGC  in  our 
Makassar  Strait  group  since  it  is  geographically  close,  though  it  is  in  the  nearby 
Flores  Sea. 

Western  Pacific  Warm  Pool  SSTs  warmed  continuously  throughout  the 
deglacial  as  noted  previously  [Kiefer  and  Kienast,  2005].  This  continuous  warming 
lacked  the  characteristic  millennial  scale  variability  of  Antarctic  records.  Western 
Pacific  Warm  Pool  SSTs  peaked  at  ~10,000  yr  BP,  and  cooled  slightly  during  the 
remainder  of  the  Holocene  (Figures  3,4). 

Deglacial  SST  evolution  in  the  Southern  Makassar  Strait  differed  from  the 
Western  Pacific  Warm  Pool.  The  three  longer  records  suggest  a  plateau  in  warming 
that  began  at  ~14,000  yr  BP  and  lasted  ~2,000  yr.  These  sites  then  warmed  until 
~8,000  yr  BP  (Figure  3,4).  SSTs  in  the  Makassar  Strait  cool  only  slightly  after 
peaking  in  the  mid-Holocene. 

In  the  ITF  outflow  sites,  SSTs  rose  from  ~20,000  yr  BP  (Figures  3,4)  to 
~15,000  yr  BP,  remained  constant  between  ~15,000  yr  BP  and  ~12,000  yr  BP, 
peaked  at  ~  10,000  yr  BP  and  then  cooled  by  about  1°C  into  the  late  Holocene.  The 
temperature  pattern  of  the  outflow  records  differs  from  the  pattern  in  other  areas  of 
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the  eastern  Indian  Ocean.  For  example,  SST  reconstructions  from  sites  to  the  north 
of  the  outflow,  off  the  coasts  of  western  Java  and  central  Sumatra,  differ  in  that  they 
show  continuous  warming  from  20,000  yr  BP  until  ~5,000  yr  BP  and  lack  a 
Holocene  cooling  [Mohtadi  etal,  2010]. 

Discussion: 

A  strong  Antarctic-type  signal  [e.g.  Kawamura  et  ah,  2007]  is  seen  in  the 
outflow  sites  (Figure  4).  Not  only  did  warming  at  each  of  the  sites  cease  during  the 
Antarctic  Cold  Reversal  (Figure  3),  but  the  sites  reached  peak  temperature  at 
approximately  the  same  time  as  Antarctica,  and  then  cooled  in  concert  with  the 
Indian  Ocean  sector  of  Antarctica.  In  fact,  the  outflow  passage  SST  stack  correlates 
more  strongly  with  the  Dome  Fuji  record  than  it  does  with  either  the  WPWP  or 
Makassar  Strait  stack  (Table  S3).  This  deglacial  and  Holocene  temperature  pattern  is 
also  suggested  in  a  SST  reconstruction  from  the  Chatham  Rise,  east  of  New  Zealand 
[Pahnke  etal,  2003]  (Figure  4).  The  close  correspondence  between  all  these  records 
suggests  a  strong  high-latitude  southern  hemisphere  imprint  on  SSTs  of  the  outflow 
region. 

We  propose  that  the  similar  deglacial  SST  trends  in  the  ITF  outflow  and  high- 
latitude  southern  hemisphere  are  due  to  the  direct  oceanic  connection  between 
Southern  Source  Water  and  the  ITF  provided  by  the  Eastern  Passages  [Van  Aken  et 
al.,  2009].  Additionally,  while  small,  the  deep  flow  in  the  Timor  Sea  is  from  the 
Indian  Ocean  [Sprintall  etal.,  2009],  which  may  contribute  Southern  Source  Water 
to  the  ITF.  The  Southern  Source  Water  is  mixed  throughout  the  water  column  in  the 
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Banda  Sea,  and  thus  affects  ITF  outflow  SSTs.  Although  the  inflow  of  Southern 
Source  Water  through  the  deep  Eastern  Passages  only  makes  up  20%  of  the  ITF 
[Van  Aken  etal,  2009],  the  apparently  much  larger  reconstructed  temperature 
variations  at  the  higher  latitudes  (as  seen  at  the  Chatham  Rise  [Pahnke  et  al,  2003]) 
would  cause  water  sourced  from  these  high  southern  latitude  regions  to  have  a 
disproportionate  influence  in  the  southern  Indonesian  region. 

We  test  whether  our  hypothesis  of  a  Southern  Source  Water  influence  at  the 
outflow  passages  is  consistent  with  modern  observations  of  the  ITF.  We  estimate 
the  portion  of  Southern  Source  Water  (SSW)  required  to  cause  the  Holocene  SST 
decrease  in  the  outflow  passages  using  the  following  equations: 

Equation  1:  TiTFoutflow(ti)  —  ccsswTssw(ti)  +  2ctiTi(ti) 

Equation  2:  TiTFoutflow(t2)  —  cisswTssw(t2)  +  2ctiTi(t2) 

Where  TiTFoutflow(ti)  is  the  outflow  temperature  at  time  q.  This  temperature  is 
determined  by  the  water  temperature  at  all  source  regions  and  the  portion  of  water 
(a)  that  comes  from  each  source  region.  Using  the  Chatham  Rise  record  to  estimate 
the  temperature  of  the  Southern  Source  Water  (Tssw),  we  estimate  the  portion  of 
SSW  (oissw)  required  to  cause  the  Holocene  cooling  in  the  outflow  passages.  The  SSTs 
at  all  other  source  regions  are  assumed  unchanged.  Equation  2  is  subtracted  from 
Equation  1  to  obtain: 

ATitf  outflow  —  ctsswAT  SSW 

The  temperature  change  from  the  Early  to  Late  Holocene  is  ~1.2°C  at  the 
outflow  passages  (ATitf  outflow),  and  ~3°C  at  the  Chatam  rise  (ATSSW)  [Pahnke  etal, 
2003].  This  would  require  the  ITF  to  be  composed  of  40%  Southern  Source  Water, 
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which  is  much  higher  than  the  20%  (2.5  Sv  of  a  total  13  Sv  transport)  suggested  by 
the  INSTANT  monitoring  program  [Gordon  etal,  2010;  VanAken  etal,  2009].  We 
note,  however,  that  a  Holocene  TEXs6  record  off  the  coast  of  West  Antarctica  shows 
a  6°C  Holocene  cooling  [Shevenell  etal.,  2011].  Weighting  the  Chatham  Rise  record 
and  the  West  Antarctica  record  equally  suggests  a  4.5°C  cooling  in  Southern  Source 
Water.  This  reduces  the  portion  of  Southern  Source  Water  required  to  ~25%.  If  the 
other  source  regions  for  the  ITF  in  fact  cooled  slightly  (the  WPWP  average  cools  by 
0.6°C,  Fig  4),  then  even  less  Southern  Source  Water  would  be  required. 

In  summary,  if  the  Chatham  Rise  is  representative  of  Southern  Source  Water 
temperature  entering  the  ITF,  and  the  temperature  at  the  rest  of  the  source  regions 
to  the  ITF  did  not  change  over  the  course  of  the  Holocene,  our  proposed  mechanism 
is  not  consistent  with  the  most  recent  modern  observations  of  the  ITF.  However,  if 
we  consider  additional  temperature  records  from  Southern  Source  regions  and  the 
apparent  cooling  of  the  WPWP,  then  our  proposed  mechanism  is  consistent  with 
modern  estimates  of  ITF  transport. 

Though  Makassar  Strait  SSTs  appear  to  follow  an  Antarctic  pattern  through 
the  deglacial,  during  the  Holocene,  Makassar  Strait  SST  warms  compared  to  both  the 
outflow  passages  and  Antarctic  trends.  We  hypothesize  that  the  warming  of  the 
Makassar  Strait  is  due  to  a  change  in  the  influence  of  Banda  Sea  water.  In  the 
modern,  some  Banda  Sea  water  enters  the  Makassar  Strait  during  the  southeast 
monsoons  [Fallon  and  Guilderson,  2008;  Gordon  etal.,  2003],  We  propose  that  prior 
to  the  flooding  of  the  Sunda  Shelf  there  was  a  greater  influx  of  Banda  Sea  water. 
After  the  flooding  of  the  Sunda  Shelf,  the  influence  of  the  South  China  Sea  increased 
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the  residence  time  of  the  surface  waters  within  the  Makassar  Strait,  which  warmed 


the  Makassar  Strait  relative  to  the  ITF  outflow  passages.  We  suggest  that  the 
reduction  of  surface  transport  through  the  Makassar  Strait  trapped  a  portion  of  the 
heat  that  was  previously  exported,  leading  to  a  surface  warming  in  the  Makassar 
Strait. 

We  draw  on  our  new  data,  modern  observations,  and  a  modeling  study  to 
suggest  that  there  was  no  discernable  effect  of  the  flooding  of  the  Sunda  Shelf  on  the 
SST  of  the  ITF  outflow.  We  note  that  both  before  and  after  the  Sunda  Shelf  flooding, 
outflow  temperatures  closely  follow  Antarctic  temperature  trends. 

As  described  above,  today,  the  South  China  Sea  Throughflow  into  the 
southern  Makassar  Strait  reduces  surface  transport  through  the  Makassar  Strait. 
Results  from  an  Ocean  General  Circulation  Model  suggest  that  this  reduction  in 
surface  transport  leads  to  a  0.18  PW  (1  PW  =  1015  Watts)  decrease  in  heat  transport 
through  the  Makassar  Strait,  lowering  the  Makassar  Strait  transport-weighted 
temperature  by  2°C  [Tozuka  et  al.,  2007].  Here  we  hypothesize  that  this  decrease  in 
heat  transport  does  not  result  in  cooler  outflow  temperatures  as  previously 
suggested  [Gordon  etal,  2003].  Modern  observations  show  that  the  South  China  Sea 
Throughflow  is  a  source  of  0.2  PW  of  heat  to  the  ITF  [Qu  etal,  2006],  enough  to 
compensate  for  the  decreased  heat  transport  through  the  Makassar  Strait,  implying 
a  zero  net  balance  for  the  heat  budget  of  the  ITF  outflow.  Therefore,  we  propose  that 
the  opening  of  the  Java  Sea  that  occurred  with  the  flooding  of  the  Sunda  Shelf  at 
~9,500  yr  BP  did  not  influence  ITF  outflow  surface  temperatures.  Rather,  we 
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postulate  that  the  deglacial  and  Holocene  trends  in  the  ITF  outflow  SSTs  are  due  to 
the  influence  of  high  latitude  southern  hemisphere  temperature. 

We  explore  alternate  hypotheses  to  explain  our  temperature  reconstructions, 
some  of  which  we  are  able  to  rule  out,  while  others  require  further  work.  These 
hypotheses  include  an  ENSO-like  teleconnection,  changes  in  local  insolation, 
changes  in  wind  driven  mixing,  and  changes  in  tidal  mixing. 

It  is  possible  that  a  direct  atmospheric  connection  with  Antarctica  [Visser  et 
al.,  2003]  is  the  source  of  the  deglacial  SST  variability  in  the  ITF  outflow  records. 
However,  the  absence  of  a  significant  correlation  between  ENSO  and  SSTs  in  the 
outflow  region  during  the  instrumental  era  (Figure  S4)  strongly  suggests  that  an 
ENSO-Antarctica  teleconnection  is  not  the  source  of  the  connection  between  the  ITF 
and  Antarctica. 

Changes  in  radiative  forcing  via  either  mean  annual  insolation  or  CO2  are 
unlikely  to  be  the  source  of  the  SST  variability  in  the  outflow  region.  Firstly,  we  note 
that  all  three  regions  (WPWP,  the  Makassar  Strait,  and  the  Outflow  Passages)  would 
experience  similar  radiative  forcing  as  they  are  geographically  close.  Additionally, 
peak  SSTs  in  our  outflow  region  occur  at  10,000  BP,  coincident  with  a  minimum  in 
mean  annual  insolation  (Fig  S5)  [Laskar  et  al.,  2004],  Similarly,  CO2  values  rise 
during  most  of  the  Holocene  [. Monnin  etal,  2004],  so  CO2  greenhouse  warming 
would  not  explain  the  cooling  in  our  outflow  passages. 

A  change  in  wind-driven  upwelling  is  another  potential  source  of  variability 
in  the  outflow  passages.  The  percentage  of  Globigerina  bulloides  ( G .  bulloides) 
relative  to  the  total  number  of  planktic  foraminifera  was  used  as  a  proxy  for 
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upwelling,  with  a  higher  percent  of  G.  bulloides  corresponding  to  more  upwelling 
[Prell,  1984].  The  percent  of  G.  bulloides  in  both  69-3  and  a  nearby  site  off  the 
southern  coast  of  Java  (10053-7)  [Mohtadi  etal,  2011b]  suggest  there  is  a  peak  in 
upwelling  at  ~  10,000  BP,  with  a  decrease  in  upwelling  throughout  the  Holocene 
(Fig  S5).  If  upwelling  were  the  primary  control  on  SSTs  in  the  outflow  passages,  we 
would  expect  coolest  SSTs  at  10,000  BP  and  increasing  SSTs  throughout  the 
remainder  of  the  Holocene.  However,  our  SST  reconstructions  show  the  opposite 
trend. 

The  transformation  of  water  properties  within  the  Indonesian  Sea  has  been 
attributed  to  tidal  mixing  [Ffield  and  Gordon,  1996],  Modeling  results  suggest  that 
tides  cause  a  majority  of  the  mixing  that  occurs  within  the  ITF  [Koch-Larrouy  et  al, 
2008],  Changes  in  sea  level  [e.g.  Clark  and  Mix,  2002]  likely  altered  the  amount  of 
tidal  mixing,  and  in  particular  how  tidal  dissipation  was  partitioned  between  the 
continental  shelves  and  the  deep  ocean  [Egbert  et  al.,  2004;  Wunsch,  2005].  The 
degree  to  which  changes  in  sea  level  may  have  changed  the  mixing  within  the  ITF  is 
unknown.  Even  today,  estimates  of  tidal  dissipation  within  the  ITF  region  suffer 
from  large  uncertainties  [Ray  et  al.,  2005].  In  the  modern,  much  of  the  tidal  mixing 
occurs  when  the  main  flowpath  of  the  ITF  is  over  shallow  sills  or  narrow  passages 
[Hautala  etal,  1996;  Ray  etal,  2005].  Since  the  vertical  mixing  at  these  choke 
points  is  influenced  by  the  details  of  sill  topography  [Hatayama,  2004]  past  sea  level 
differences  may  have  changed  the  location  and/or  amount  of  tidal  mixing.  The 
degree  to  which  mixing  within  the  Indonesian  Seas  changed  in  the  past  could  be 
empirically  evaluated  with  benthic  foraminifera  temperature  reconstructions  along 
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depth  transects.  Within  the  modern  Banda  Sea,  the  entire  water  column  is  mixed, 
but  the  intensity  of  local  insolation  restratifies  the  water  column  with  respect  to 
temperature.  In  the  modern,  temperatures  are  uniform  only  below  1,000m 
[Tomczak  and  Godfrey,  1994].  A  change  in  depth  of  this  temperature  gradient  would 
suggest  a  change  in  mixing  intensity.  This,  along  with  modeling  the  influence  of 
reduced  sea  level  on  the  mixing  that  occurs  at  the  ITF's  choke  points  would  give 
insight  into  the  potential  role  of  changes  in  tidal  mixing  on  ITF  water  properties. 

Conclusions: 

Our  study  demonstrates  a  persistent  link  between  SSTs  in  the  eastern  ITF 
outflow  region  and  high  latitude  southern  hemisphere  temperature  during  the  last 
20,000  years,  which  we  suggest  is  through  an  oceanic  pathway.  Southern  Source 
Water  enters  the  Indonesian  Seas  primarily  through  the  Eastern  Passages  and  to  a 
lesser  extent  through  the  Timor  Stait.  Subsequently,  it  is  mixed  upward  in  the  Banda 
Sea,  where  it  enters  the  main  ITF.  This  direct  oceanic  connection  between  Antarctic 
temperature  and  tropical  SSTs  would  allow  the  transmission  of  southern  high 
latitude  climate  signals  to  the  ITF  outflow  and  into  the  Indian  Ocean.  Our  proposed 
mechanism  of  a  direct  oceanic  connection  appears  as  a  testable  hypothesis  of  the 
modern  controls  on  the  ITF  outflow. 

Our  results  also  highlight  the  importance  of  the  South  China  Sea  Throughflow 
for  sea  surface  temperatures  in  the  ITF  region.  Prior  to  the  establishment  of  the 
South  China  Sea  Throughflow,  southern  hemisphere  climate  had  a  greater  impact  on 
Makassar  Strait  SST.  The  South  China  Sea  Throughflow,  however,  changed  the 
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surface  circulation  so  that  the  Makassar  Strait  was  no  longer  strongly  influenced  by 
the  southern  hemisphere.  Makassar  Strait  SSTs  warmed,  and  remained  elevated 
throughout  the  Holocene. 
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Figure  1: 

Map  of  the  Indonesian  Throughflow  area.  The  primary  flow  path  of  surface  water 
is  shown  in  gray  arrows.  The  freshwater  flow  from  the  Java  Sea  is  indicated  by  a 
dashed  arrow.  The  flow  of  deep  water  through  the  Eastern  Passages  is  shown  by 
blue  arrows.  Coretops  are  indicated  by  filled  circles,  new  downcore  records  are 
indicated  with  diamonds,  and  triangles  indicate  downcore  data  from  previously 
published  records.  Black  triangles  correspond  to  Western  Pacific  Warm  Pool  sites, 
red  triangles  and  diamonds  with  Maskassar  Strait  sites,  and  blue  triangles  and 
dimonds  with  outflow  passage  sites.  This  color  coding  corresponds  with  the  records 
shown  in  figures  3-4.  Background  color  shows  mean  annual  SST  from  the  World 
Ocean  Atlas  05  dataset  [Locarnini  etal.,  2006]. 
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Mg/Ca  excess  and  salinity 
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Figure  2: 

Coretop  excess  Mg/Ca  versus  salinity.  Excess  Mg/Ca  [Arbuszewski  etal.,  2010]  is 
calculated  by  subtracting  the  expected  Mg/Ca  based  on  observed  sea  surface 
temperatures  from  measured  Mg/Ca.  These  data  suggest  that  the  influence  of 
salinity  on  Mg/Ca  is  small  below  36. 
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Figure  3: 

Mg/Ca-based  sea  surface  temperature  estimates  from  this  paper,  along  with 
previously  published  data.  All  records  are  converted  to  temperature  using  a 
multispecies  equation  [Anand  et  al.,  2003],  except  for  MD41,  for  which  we  used  the 
published  sea  surface  temperatures  [Rosenthal  et  al.,  2003]  (see  also  the 
supplementary  discussion).  Thin  lines  connect  raw  data,  bold  lines  are  500  year 
averages  from  the  Western  Pacific  Warm  Pool  (black),  the  Makassar  Strait  (red),  and 
the  outflow  passages  (blue).  Mean  annual  sea  surface  temperatures  measured  by 
satellite  are  indicated  for  each  core  site  (+)  and  modern  sea  surface  temperature 
estimates  from  multi-core  samples  are  indicated  by  diamonds. 
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Figure  4: 

Comparison  of  sea  surface  temperature  records  to  Southern  Hemisphere 
records.  Mg/Ca-based  sea  surface  temperature  records,  placed  into  500-year,  non¬ 
overlapping  bins  and  averaged,  with  the  shading  around  the  bins  representing  the 
propagated  error  that  includes  both  the  analytical  error  and  the  error  in  the  Mg/Ca - 
temperature  calibration  (see  supplementary  discussion  for  a  detailed  discussion).  A 
Mg/Ca-based  sea  surface  temperature  record  from  the  Chatham  Rise  [Pahnke  etal, 
2003]  is  shown  in  purple  (vertical  bar  is  uncertainty,  see  supplementary  discussion 
for  details)  and  5180  from  the  Indian  Ocean  sector  of  Antarctica  (Dome  Fuji)  is 
shown  in  black  [Kawamura  etal.,  2007].  Black  rectangles  denote  the  Antarctic  Cold 
Reversal  (timing  inferred  from  Dome  Fuji  data)  and  the  timing  of  the  Sunda  Shelf 
flooding  [Sathiamurthy  and  Voris,  2006]. 
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Supplementary  material: 

I.  Coretop  study  -  The  influence  of  salinity  on  the  Mg/Ca-temperature 
relationship 

All  our  samples  for  the  coretop  study  were  taken  with  a  multi-coring  device 
on  three  cruises  in  the  Indonesian  Throughflow  region  and  the  Eastern  Indian 
Ocean  (Figure  1,  dots;  also  see  Appendix  B). 

The  influence  of  salinity  on  Mg/Ca  was  explored  by  combining  our  new 
coretop  data  with  previously  published  G.  ruber  data.  All  our  samples  were  prepared 
using  a  full  trace  metal  cleaning,  including  reductive  and  oxidative  steps  [Boyle  and 
Keigwin,  1985/6;  Rosenthal  et  ah,  1997].  However,  some  of  the  previously  published 
data  were  from  samples  that  were  not  cleaned  with  a  reductive  step.  To  account  for 
this,  we  added  10%  to  the  Mg/Ca  of  all  samples  cleaned  using  both  oxidation  and 
reduction  following  Abruszewski  etal.  [2010]. 

We  evaluate  how  well  the  Mg/Ca  measurements  predict  temperature  by 
calculating  the  "excess  Mg/Ca "  [Arbuszewski  etal.,  2010]. 

Excess  Mg/Ca  =  measured  Mg/Ca  -  predicted  Mg/Ca 

In  order  to  determine  predicted  Mg/Ca  we  used  the  Anand  etal.  [2003] 
multispecies  calibration  in  conjunction  with  modern  data  from  either  World  Ocean 
Atlas  2005  (WOA05)  temperature  data  [Locarnini  etal,  2006],  or  the  temperature 
data  presented  in  the  original  publication. 
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Anand  multispecies  equation:  =  0.38<?009/ 

Ca 

When  excess  Mg/Ca  is  compared  with  salinity  (again  either  WOA05  data  [Antonov  et 
al,  2006]  or  the  salinity  from  the  original  publication  was  used),  there  is  no 
correlation  below  36  (Figure  2,  r2  =  0.001).  This  finding  supports  a  recent  paper 
which  found  a  salinity  threshold  of  35  [Arbuszewski  et  al,  2010],  but  conflicts  with 
another,  which  suggests  that  the  influence  of  salinity  on  Mg/Ca  is  a  concern  at  all 
salinities  [Mathien-Blard  and  Bassinot,  2009].  This  discrepancy  is  likely  due  to  two 
reasons.  First  there  were  very  few  Mg/Ca  data  from  relatively  fresh  sites  in  previous 
publications.  The  new  data  presented  here  address  this  issue.  Secondly,  one  must 
use  caution  if  using  518Ocaicite  and  518Oseawater  measurements  to  estimate  calcification 
temperature. 

We  did  not  consider  calculating  the  calcification  temperature  based  on 
5180Caicite  and  5180Seawater  because  of  the  very  limited  518Oseawater  data  in  our  study 
area.  We  suggest  the  gridded  database  of  518Oseawater  [LeGrande  and  Schmidt,  2006] 
is  not  appropriate  to  use  in  this  case  because  it  uses  both  P04*  and  salinity  to 
estimate  518Oseawater  in  areas  with  few  518Oseawater  measurements.  Therefore,  the 
most  robust  manner  to  evaluate  salinity’s  effect  on  Mg/Ca  is  to  calculate  excess 
Mg/Ca  with  modern  temperature  data  and  not  estimated  isotopic  temperature.  This 
avoids  the  circularity  of  including  a  salinity-based  5180  seawater  estimate  in  a 
comparison  with  salinity. 

We  also  explored  the  influence  of  core  depth  on  our  samples  by  examining 
the  relationship  between  Mg/Ca  excess  and  core  depth  (Figure  S6,  R2  =  0.005).  The 
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lack  of  correlation  between  depth  and  Mg/Ca  suggests  that  preferential  dissolution 
is  not  the  source  of  the  Mg/ Ca  excess. 

Modern  SST 

Modern  SST  data  were  downloaded  from  NASA’s  Ocean  Color  Radiometry 
Online  Visualization  and  Analysis  Global  Monthly  Products 
(http://reason.gsfc.nasa.gov/Giovanni/). 

Our  Mg/Ca-based  modern  temperature  estimates  for  the  Makassar  Strait  and 
ITF  outflow  regions  are  consistently  cooler  than  satellite-derived  mean  annual 
temperature,  though  they  are  all  within  the  seasonal  range.  This  could  be  due  to  a 
growth  preference  during  the  upwelling  season  and/or  because  G.  ruber  may  calcify 
throughout  the  upper  30  m  [Wang,  2000],  even  in  regions  with  a  shallow  mixed 
layer. 

II.  Downcore  records 

Temperature  calibration  (cleaning  method  adjustments) 

Previous  work  [Arbuszewski  et  al,  2010;  Barker  et  al,  2003;  Rosenthal  et  al, 
2004]  has  shown  that  the  reductive  step  in  Mg/Ca  cleaning  reduces  the  Mg/Ca  of  a 
sample  by  approximately  10%.  The  Mg/Ca-temperature  calibration  [Anand  etal, 
2003]  we  used  was  based  on  samples  that  were  cleaned  without  a  reductive  step.  To 
account  for  this  we  add  10%  to  the  Mg/Ca  of  samples  that  were  cleaned  with  the  full 
(with  oxidative  and  reductive  steps)  cleaning  method  [Boyle  and  Keigwin,  1985/6; 
Rosenthal  etal,  1997]. 
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Original  Anand  etal.  (2003)  multispecies  equation  [Anand  etal.,  2003]: 


Mg 

Ca 


0.38e°°9r 


Revised  equation  for  the  full  cleaning  method: 


Mg 

Ca 


+  0.1 


Mg 

Ca 


03SeOWT 


or 


Mg 

Ca 


0.345e009r 


Errors  on  our  temperature  estimates 

There  are  two  primary  sources  of  error  in  converting  Mg/Ca  to  temperature. 
The  first  is  within  the  Mg/Ca-temperature  calibration  equation,  and  the  second  is 
the  reproducibility  of  Mg/Ca  measurement. 

We  used  the  values  and  associated  errors  on  the  Anand  multispecies  equation 
[Anand  et  al,  2003]: 

— ^  =  beaT 
Ca 

b  =  0.38  (+/-  0.02);  a  =  0.090  (+/-  0.003)  °C  1 
or  for  the  full  cleaning  equation: 
b  =  0.345(+/-  0.02);  a  =  0.090  (+/-  0.003)  T'1 

Our  estimate  of  the  reproducibility  of  our  Mg/Ca  measurements  is  based  on 
comparing  samples  from  the  same  depth.  The  foraminifera  from  these  duplicates 
were  not  homogenized  prior  to  running,  so  this  should  be  a  conservative  estimate  of 
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reproducibility  based  on  both  the  natural  variability  that  a  sample  may  contain  and 
instrumental  error.  The  standard  deviation  of  the  Mg/Ca  difference  between  17 
such  duplicates  was  0.21. 

These  two  sources  of  error  were  propagated  [Bevington  and  Robinson,  1992]  as 
follows: 

2  dSST  ,  dSST  ,  ,  dSST  , 

<rSST -(—aj  ♦<— 


where 


dSST 

da 

dSST 

db 


1  M§/ 

— ln(  Ca)  , 
a  b 

J_ 

ab 


and 


dSST  1  1 

,Mg/  ~  a  Mg/ 
/Ca  /Ca 


For  each  core,  if  there  was  a  single  data  point  in  the  500  year  bin,  the  error 
was  simply  osst  ,  as  determined  above.  If  there  was  more  than  one  temperature 
value  within  the  500  year  bin,  the  errors  were  calculated  in  the  following  manner 
[Bevington  and  Robinson,  1992],  where  n  is  the  number  of  data  points  in  each  bin, 


1  " 

CtbinSST  =  2  STi 


The  errors  on  each  500  year  bin  in  the  regional  stack  were  determined  by 
summing  the  square  of  the  errors  (o2binssi)  from  each  individual  core. 
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^ averageSST 


1  n 


binSSTi 


i= 1 


The  error  on  the  Chatham  Rise  core  (MD97-2120)  was  determined  by 
propagating  the  reported  analytical  relative  uncertainty  (5.8%)  [Pahnke  etal, 

2003],  and  the  reported  calibration  error  (±0.8°C)  [Mashiotta  etal,  1999]  which 
gives  a  total  error  of  ±1°C. 

Our  approach  to  error  propagation  is  conservative  in  that  a  portion  of  the 
calibration  error  is  due  to  the  analytical  uncertainty  in  the  Mg/Ca  measurements.  By 
also  including  our  analytical  uncertainty  in  our  error  propagation  we  are  likely 
overestimating  the  errors. 

Radiocarbon  calibration  and  age  model  construction 

All  radiocarbon  dates  were  converted  into  calendar  ages  using  the  Calib  6.0 
program  and  the  Marine  09  calibration  [Reimer  etal.,  2009],  No  local  reservoir 
correction  was  applied  as  modern  measurements  suggest  the  regional  AR  is  close  to 
zero  [Bowman,  1985;  Southon  etal.,  2002].  The  median  probability  was  used  for 
constructing  age  models.  All  dates  used,  and  their  associated  errors,  can  be  found  in 
Appendix  A. 

In  cases  where  the  first  standard  deviation  of  calendar  ages  had  two  modes, 
we  reported  a  range  of  calendar  ages  that  encompasses  both  peaks,  as  a 
conservative  estimate.  In  cases  where  we  used  a  polynomial  or  linear  fit,  we 
estimated  error  at  the  age  control  points  by  summing  the  squares  of  both  the 


51 


radiocarbon  calibration  errors  (first  standard  deviation)  and  an  estimate  of  the 
polynomial’s  error  at  these  points  (also  first  standard  deviation).  These  values  are 
included  in  Appendix  A. 

For  cores  23GGC,  70GGC,  136GGC,  69-3,  MD62,  and  MD65,  age  models  were 
constructed  by  linearly  interpolation  between  dates.  The  age  models  for  MD41, 
MD78,  and  MD81  were  constructed  by  utilizing  the  equations  listed  below. 

Polynomial  for  MD97-2141 

y  =  -0.00000000 124x5  +  0.000001404x4  -  0.0007356x3  +  0.1782x2  +  31.90x+4314 

R2  =  0.979 

Linear  fit  for  MD01-2378  (the  intercept  was  set  at  1050  years,  to  produce  a  realistic 
coretop  age) 
y  =  48.884x+  1050 
R2  =  0.991 

Polynomial  for  MD98-2181 

y  =  -0.0000000022  62x4+0.000004298x3+0.006586x2+4.412x+10  7 
R2  =  0.997 

There  is  also  uncertainty  in  the  age  model  of  the  Dome  Fuji  record  of 
approximately  +/-  1,000  years  [Kawamura  etal,  2007].  For  a  detailed  explanation 
of  the  errors  associated  in  that  age  model,  refer  to  that  paper’s  supplementary 
material. 
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Figure  SI:  Sea  surface  salinity  and  height  of  the  ITF.  January  through  March  and 
July  through  September  climatologies  were  creating  using  the  SODA  reanalysis 
product  [Carton  and  Giese,  2008].  We  note  caution  is  required  when  using  a 
reanalysis  product  and/or  climatologies  as  these  products  create  the  appearance  of 
complete  data  coverage,  which  is  not  the  case,  particularly  in  an  area  such  as  the 
Indonesian  Throughflow,  which  has  relatively  sparse  measurements. 
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Figure  S2:  Boxes  used  to  calculate  sea  surface  height  difference  for  Fig  S3. 
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Ocean  Data  View 


Time 

Figure  S3:  Difference  in  sea  level  height.  The  difference  in  sea  level  height  [Carton 
and  Giese,  2008]  between  the  southern  Java  Sea  and  the  Banda  Sea  (red  line)  and  the 
southern  Java  Sea  and  the  Flores  Sea  (blue  line).  The  sea  level  of  the  southern  Java 
Sea  is  nearly  always  higher  than  the  Flores  Sea.  Boxes  used  to  average  the  sea  level 
height  are  shown  in  Fig  S2. 
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Figure  S4:  SST  ENSO  correlation.  The  linear  correlation  coefficient  between  the 
NIN03  index  [Kaplan  et  al.,  1998;  Reynolds  et  al,  2002]  and  global  reconstructed 
monthly  SST  anomalies  [Smith  et  al,  2008;  Xue  et  al.,  2003]  from  January  1900  until 
December  2010.  Tan  coloring  indicates  no  significant  correlation  at  the  95%  level. 
The  lack  of  correlation  between  the  NIN03  index  and  SST  anomalies  in  the  ITF 
outflow  region  suggests  that  an  ENSO  teleconnection  cannot  be  the  source  of  SST 
variability  in  the  region. 
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Figure  S5:  Alternate  hypotheses  to  explain  the  cooling  in  the  outflow  passages. 

An  upwelling  indicator  (%  G.  bulloides )  from  a  core  off  the  southern  coast  of  Java 
[Mohtadi  etal,  2011]  and  69-3  suggest  either  a  no  change  in  upwelling  or  a 
maximum  upwelling  when  SSTs  were  warmest  at  the  outflow  passages.  Mean 
annual  insolation  at  the  equator  has  minimum  during  the  early  Holocene,  which 
cannot  explain  the  peak  in  temperatures  at  the  outflow  passages,  and  CO2  [Monnin 
etal.,  2004]  rises  throughout  the  Holocene,  while  SSTs  at  the  outflow  passages 
decrease.  The  change  in  sea  level  [Clark  and  Mix,  2002]  may  influence  tidal  mixing 
throughout  the  deglaciation. 
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Excess  Mg/Ca  versus  depth 
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Figure  S6:  Excess  Mg/Ca  versus  depth.  Excess  Mg/Ca  is  calculated  by  subtracting 
climatological  [Anand  etal.,  2003]  SSTs  from  Mg/Ca  based  SST  estimates.  The  lack  of 
correlation  between  core  depth  and  excess  Mg/Ca  suggests  that  there  is  no 
preferential  dissolution  in  our  coretop  data. 
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Core  ID,  location 

Core  ID  in 

text 

Lat./Long. 

Water 

depth 

(m) 

Sample 

interval 

(yrs) 

Reference 

BJ8-03  23GGC,  S 

Makassar  Strait 

23GGC 

5°  12'S,  117“ 
29'E 

1189 

180 

this  study 

BJ8-03  70GGC,  SW 
Sulawesi,  Makassar  Strait 

70GGC 

3“  34'S,  119“ 
23 'E 

482 

125 

Linsley  et  al.,  2010 

BJ8-03  136GGC,  Flores 

Sea 

136GGC 

5“  56'S,  120“ 
14'E 

526 

200 

this  study 

GeoB  10069-3  Savu  Sea 

69-3 

9“  36'S,  120“ 
55'E 

1,250 

190 

this  study 

MD97-2141,  Sulu  Sea 

MD41 

8“  47'N, 

121“  17'E 

3,633 

90 

Rosenthal  et  al.,  2003 

MD98-2162,  S.  Makassar 

St. 

MD62 

4“  4 l'S,  117“ 
54'E 

1,855 

415 

Visser  et  al.,  2003 

MD98-2165,  Sumba  St. 

MD65 

9“  39'S,  118“ 
20'E 

2,100 

200 

Levi  et  al.,  2007 

MD01-2378,  Timor  Sea 

MD78 

13“  05'S, 

121°  47'E 

1,783 

125 

Xu  et  al.,  2008; 
Sarnthein  et  al.,  2011 

MD98-2181,  W.  Pacific, 
Mindanao 

MD81 

6“  27'N, 

125“  50'E 

2,114 

50 

Stott  et  al.,  2007 

Table  SI:  Corelocations. 
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Core 

Reference 

Cleaning  Method 

Equation  used 

BJ8-03  23GGC 

this  study 

oxidative  and 
reductive 

modified  Anand 

BJ8-03  70GGC 

Linsley  etal,  2010 

oxidative  and 
reductive 

modified  Anand 

BJ8-03 

136GGC 

this  study 

oxidative  and 
reductive 

modified  Anand 

GeoB10069-3 

this  study 

oxidative  and 
reductive 

modified  Anand 

MD97-2141 

Rosenthal  etal,  2003 

oxidative 

SST  from  original 
publication 

MD98-2162 

Visser  etal,  2003 

oxidative 

Anand 

MD98-2165 

Levi  et  al.,  2007 

oxidative 

Anand 

MD01-2378 

Xu  etal.,  2008 

oxidative  and 
reductive 

modified  Anand 

MD98-2181 

Stott  etal,  2007 

oxidative 

Anand 

Table  S2:  Summary  of  cleaning  method  and  Mg/Ca  -  temperature  equation 
used  on  each  core. 
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WPWP 

Mak 

Outflow 

Dome  Fuji 

WPWP 

1 

0.9626 

0.9582 

0.9717 

Mak 

0.9626 

1 

0.9414 

0.9335 

Outflow 

0.9582 

0.9414 

1 

0.9825 

Dome  Fuji 

0.9717 

0.9335 

0.9825 

1 

Table  S3:  Correlation  matrix  of  SST  stacks  and  Dome  Fuji  6180.  All  correlations 
are  significant  (n  =  45;  p  «  0.01). 
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Chapter  3 


The  deglacial  hydrologic  variability  of  the  Tropical  Pacific  and  Indian  Oceans 
inferred  from  8180  of  seawater  reconstructions 


Abstract 

The  hydrologic  cycle  of  the  tropics  can  change  rapidly  and  dramatically, 
affecting  the  amount  of  water  vapor,  and  potentially  temperatures,  globally. 
Previous  studies  have  demonstrated  variability  in  the  tropical  hydrologic  cycle 
during  the  last  deglaciation  and  both  high  and  low  latitude  origin  of  the  variability 
has  been  proposed.  Here  we  use  the  Mg/Ca-based  temperature  and  the  5180  of 
calcite  of  planktic  foraminifera  to  solve  for  the  5180  of  seawater.  We  combine  new 
data  with  previously  published  5180  of  seawater  reconstructions  from  throughout 
the  tropical  Eastern  and  Western  Pacific  and  Eastern  Indian  Oceans  in  order  to 
evaluate  the  possible  mechanisms  controlling  hydrologic  changes.  Our  synthesis  of 
the  5180  of  the  published  and  new  seawater  reconstructions  suggests  that 
hydrologic  conditions  at  the  last  glacial  maximum  and  late  Holocene  were  similar. 
However,  5180  of  seawater  enrichments  occurred  at  most  of  the  sediment  sites 
during  the  Younger  Dryas  and  Heinrich  Stadial  1.  Similar  enrichments  in  the  East 
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and  West  tropical  Pacific  seem  to  rule  out  forcing  by  an  El  Nino-Southern 
Oscillation-like  mechanism.  The  spatial  pattern  of  change  is  instead  consistent  with 
a  migration  of  the  Intertropical  Convergence  Zone. 

1.  Introduction 

The  temperature  of  the  tropics  appears  to  have  been  relatively  stable  over 
the  last  20  kyr,  compared  to  the  high  latitudes,  which  were  characterized  by  large 
temperature  fluctuations  [e.g.  Alley  and  Clark,  1999;  Dansgaard  etal,  1989; 
Ruddiman  and  Mcintyre,  1981].  The  large,  rapid  temperature  fluctuations 
documented  in  the  high  latitudes  either  do  not  have  temperature  counterparts  in 
the  tropics  [Kiefer  and  Kienast,  2005;  Stott  et  al,  2007]  or,  if  they  do,  their  amplitude 
is  diminished  [Kiefer  and  Kienast,  2005;  Lea  etal.,  2003;  Shakun  and  Carlson,  2010]. 
On  the  other  hand,  recent  work  has  highlighted  dramatic  hydrologic  responses  that 
occurred  in  concert  with  high  latitude  climate  change  (e.g.  [Lynch-Stieglitz  et  al, 
2011;  Partin  etal.,  2007;  Stager  et  al.,  2011;  Tierney  et  al.,  2008;  Yuan  etal.,  2004]). 

During  the  last  deglacial,  two  prominent  Northern  hemisphere  cooling  events 
occurred  —  Heinrich  Event  1  (HI;  ~17.5  to  16  kyr  BP)  and  the  Younger  Dryas 
(~12.9  to  11.7  kyr  BP)  [Alley  and  Clark,  1999].  Heinrich  Events,  first  identified  by 
large  amounts  of  ice  rafted  debris  in  sediment  cores  [Heinrich,  1988;  Hemming, 
2004],  appear  to  coincide  with  1)  the  end  of  cooling  stadials  as  recorded  by  the  6180 
of  Greenland  ice  and  2)  very  cold  Atlantic  sea  surface  temperatures  (SSTs)  [Bond  et 
al,  1993].  Similarly,  the  Younger  Dryas  coincides  with  ice  rafted  debris  in  North 
Atlantic  sediments  and  cooler  temperatures  [Bond  et  al,  1993].  The  cause  of  these 
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abrupt  cooling  events  in  the  North  Atlantic  is  still  under  debate,  but  the  leading 
hypothesis  is  that  they  resulted  from  freshwater  input  at  the  sites  of  deepwater 
formation,  which  altered  the  Meridional  Overturning  Circulation  [Broecker,  1991]. 
The  source  of  freshwater  was  likely  iceberg  discharge  for  Heinrich  Events  [Bond  et 
ah,  1992].  Evidence  of  the  fresh  water  source  is  more  elusive  for  the  Younger  Dryas 
[Broecker,  2006],  though  recent  work  suggests  freshwater  may  have  been  routed 
through  the  Arctic  [e.g.  Murton  etal.,  2010]. 

Terrestrial  records  indicate  that  tropical  hydrological  changes  associated 
with  both  HI  and  the  Younger  Dryas  were  large  and  pervasive  [e.g.  Stager  et  ah, 
2011;  Tierney  et  al,  2008;  Wang  etal.,  2007;  Yuan  etal,  2004].  Whereas  hydrologic 
changes  over  land  have  been  attributed  in  large  part,  but  not  entirely,  to  southward 
shifts  of  the  Intertropical  Convergence  Zone  (ITCZ)  during  North  Atlantic  cold 
events  (See  for  example,  Stager  et  al.  [2011]  and  Oppo  and  Curry  [2012]  for  recent 
compilations  of  HI  terrestrial  responses),  there  is  not  a  consensus  in  studies  of 
marine  sediments.  Some  have  inferred  that  changes  in  the  mean  Pacific  state  (e.g. 
east-west  SST  gradient  changes)  were  a  dominant  control  on  tropical  Pacific 
hydrology  during  these  events  [e.g.  Levi  etal,  2007;  Rosenthal  etal,  2003;  Stott  et 
al,  2002],  while  others  have  suggested  that  the  data  are  consistent  with  a  migration 
of  the  ITCZ. 

The  5180  of  seawater  (5180  seawater  )  can  be  used  to  inform  us  as  to  past 
hydrologic  changes.  The  5180  seawater  and  salinity  often  have  a  positive  linear 
relationship  as  both  quantities  are  similarly  influenced  by  evaporation  and 
precipitation,  with  more  enriched  5180  seawater  values  corresponding  to  saltier 
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conditions.  This  linear  relationship  varies  by  region  [LeGrande  and  Schmidt,  2006] 
as  sea  ice  formation,  atmospheric  distillation,  and  river  runoff  affect  salinity 
differently  than  they  do  the  5180  of  water  [Craig  and  Gordon,  1965].  Moreover,  the 
exact  relationship  between  5180  seawater  and  salinity  has  likely  varied  in  the  past 
[LeGrande  and  Schmidt,  2011;  Oppo  etal,  2007].  Thus,  the  5180  seawater  can  best 
inform  us  about  large-scale  changes  in  the  hydrologic  cycle,  and  should  not  be 
viewed  as  a  simple  paleosalinity  proxy  [Craig  and  Gordon,  1965;  LeGrande  and 
Schmidt,  2011]. 

Here,  we  re-examine  the  influence  of  abrupt  North  Atlantic  events  on  tropical 
hydrology  from  a  marine  perspective.  We  present  a  new  high-resolution  record  of 
the  5180  of  calcite  (5180Caidte)  of  the  mixed  layer  planktic  foraminifera 
Globigerinoides  ruber  ( G .  ruber )  from  the  Savu  Sea,  which  is  supplied  by  Indonesian 
Throughflow  (ITF)  water  entering  the  Indian  Ocean  through  the  Ombai  Strait 
(Figure  1).  Using  submitted  (20ky-present)  [Gibbons  etal.,  in  revision]  and  new 
(35kyr  -  20kyr)  Mg/Ca-based  temperature  reconstructions,  and  following  previous 
work  [e.g.  Lea  et  al,  2000]  we  calculated  the  518Oseawater  at  this  site  over  the  past  35 
kyr.  Recent  studies  have  suggested  that  the  opening  of  the  passages  between 
marginal  basins  and  seas  as  sea  level  rose  during  the  last  deglaciation,  influenced 
5180seawater  and  freshwater  gradients  in  the  region  [Ding  et  al,  2002;  Linsley  et  al, 
2010;  Rosenthal  etal,  2003;  Xu  etal,  2010].  However  in  this  study,  we  focus  on 
larger  scale,  basin-wide  variations  and  the  information  that  518Oseawater  can  provide 
about  past  changes  in  the  tropical  hydrologic  cycle.  We  combine  our  new  record 
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with  published  518Oseawater  reconstructions  from  throughout  the  IndoPacific  Warm 
Pool  and  the  Eastern  Equatorial  Pacific  Ocean  (Table  1)  to  determine  the  most 
important  mechanism(s)  of  deglacial  tropical  hydrologic  change. 

2.  Materials  and  Methods 

The  new  data  we  present  was  generated  on  sediment  core  GeoB10069-3  (9° 
36’S,  120°  55’E,  1,250m;  referred  to  as  69-3),  retrieved  from  the  Savu  Basin  in  the 
Eastern  Indian  Ocean.  The  Savu  Basin  is  in  the  path  of  one  of  the  outflow  passages  of 
the  ITF.  Today,  minimum  temperatures  at  69-3  are  ~26°C  and  occur  in  the  boreal 
summer  (all  seasons  are  with  respect  to  the  Northern  Hemisphere)  and  maximum 
temperatures  are  ~31°C. 

Stable  isotope  measurements  from  69-3  were  made  at  Bremen  University  on 
a  Finnigan  MAT251.  Long-term  precision  of  5180Caidte  measurements  is  0.07%o, 
based  on  repeated  measurements  of  a  suite  of  standards.  Approximately  5-20 
individual  tests  were  used  for  each  measurement. 

Approximately  40  G.  ruber  tests  (212-300pm)  were  used  for  each  Mg/Ca 
measurement.  Samples  were  cleaned  with  a  full  trace  metal  cleaning  procedure, 
including  oxidative  and  reductive  steps  [Boyle  and  Keigwin,  1985/6;  Rosenthal  et  al, 
1997],  Samples  were  run  on  either  an  Element  ICP-MS  at  Wood  Hole  Oceanographic 
Institution  or  an  Element  XR  ICP-MS  at  Rutgers  University.  Internal  reproducibility 
based  on  repeated  measurements  of  consistency  standards  is  ~1%. 

Mg/Ca  values  were  converted  to  temperature  using  the  Anand  multispecies 
equation  [Anand  et  al.,  2003].  The  equation  was  adjusted  [Gibbons  et  al,  in  revision] 
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to  account  for  the  fact  that  the  Anand  calibration  was  based  on  foraminifera  that  had 


been  cleaned  only  with  an  oxidative  step: 

Original  Anand  Multi-species  equation:  Mg/Ca  =  0.38exp(0.090T) 

Modified  Anand  equation:  Mg/Ca  =  0.345exp(0.090T) 

Mg/Ca-based  temperature  estimates  and  5180Caidte  data  from  the  same  depth 
horizons  were  used  to  calculate  518Oseawater.  Following  Lea  etal.  [2000],  the  low- 
light  0.  universa  calibration  [Bemis  et  al.,  1998]  was  used  to  obtain  5180  seawater! 

Bemis  low-light  0.  universa  equation:  T(°C)  =  16.5  -  4.80(518Ocaidte  -  518Oseawater)  -  0.27%o 
For  20-0  kyr  before  present  (BP),  the  change  in  S180  seawater  due  to  continental 
ice  decay  was  removed  using  the  sea  level  curve  of  Clark  and  Mix  [2002]  and  the 
corresponding  change  in  the  global  5180  seawater  of  1.0%o  [Schrag  etal.,  2002] 
allowing  us  to  estimate  the  5180  of  seawater  corrected  for  ice  volume  (618Oseawater- 
icevoiume).  For  35-20kyr  BP  we  used  the  sea  leve  curve  of  Waelbroeck  et  al.  [2002]. 

We  also  compiled  previously  published  Mg/Ca  and  5180Caidte  data  from 
planktic  foraminifera  from  cores  throughout  the  Indian  and  Pacific  tropics  (Fig  1, 
Table  1).  All  records  are  based  on  data  from  G.  ruber,  with  the  exception  of  V2 1-30, 
which  utilized  Globigerinoides  sacculifer  ( G .  sacculifer).  In  our  compilation  we 
examine  the  past  20  kyr.  We  have  attempted  to  treat  the  data  from  the  cores  in  the 
same  manner  in  order  to  better  allow  us  to  compare  datasets. 

Where  possible,  Mg/Ca  was  converted  to  temperature  using  either  the 
original  Anand  multispecies  or  the  modified  version  based  on  cleaning  method  (see 
Table  1).  For  cores  TR163-22  [Lea  et  al.,  2006]  and  MD41  [Rosenthal  et  al.,  2004], 
we  used  the  SSTs  from  the  original  publication  in  order  to  account  for  the  potential 
effects  of  dissolution.  All  518Oseawater  values  were  calculated  as  described  above. 
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We  evaluate  the  518Oseawater  reconstructions  by  comparing  them  to  modern 
water  measurements.  In  each  core,  we  average  518Oseawater  samples  from  2-0  kyr  BP. 
We  plot  these  "modern"  averages  against  the  mean  annual  World  Ocean  Atlas  05 
salinity  [Antonov  etal.,  2006]  from  each  core  location.  We  compare  the  resulting 
5180seawater/salinity  relationship  with  regional  518Oseawater  and  salinity  data  compiled 
by  LeGrande  and  Schmidt  [2006]  (Fig  2).  For  the  IndoPacific  coretops,  the  slopes  of 
the  reconstructed  5180Seawater/salinity  (m  =  0.48)  and  the  modern  water 
5180seawater/salinity  (m  =  0.44)  are  within  error.  For  the  Eastern  Equatorial  Pacific, 
the  slopes  are  also  very  similar  (0.22  for  the  water  samples  and  0.31  for  the 
foraminifera  samples). 

All  age  models  are  based  on  previously  published  planktic  foraminifera 
radiocarbon  dates  with  the  exception  of  the  4  oldest  dates  from  69-3  which  are 
presented  here  (Appendix  A).  All  radiocarbon  dates  were  converted  to  calendar 
ages  using  the  Calib  6.0  program,  the  Marine  09  calibration,  and  the  global  reservoir 
age  [Reimer  et  ah,  2009].  The  calendar  age  with  the  median  probability  was  selected. 
Several  of  the  age  models  were  published  previously:  69-3,  MD62,  MD81,  MD41, 
MD65,  and  MD78  [Gibbons  et  al,  in  revision].  Age  models  for  V21-30,  43]C,  TR163- 
22,  MD65,  29-4,  38-4,  MD76,  MD70,  MD90,  and  MD38  were  created  by  linearly 
interpolating  between  derived  calendar  ages. 

We  computed  the  arithmetic  average  of  our  data  using  500-year  non¬ 
overlapping  bins.  For  each  time  interval  we  only  averaged  the  cores  with  data;  cores 
without  data  were  simply  excluded  from  the  time  bin.  Principal  component  analysis 
[ JoIIiffe ,  2002]  was  also  used  to  explore  the  common  features  of  the  S180  seawater- 
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icevolume-  All  data  were  placed  into  500-year  non-overlapping  bins.  Missing  values 
were  replaced  with  the  average  value  of  the  record  [following  Marchal  et  al.,  2002]. 
We  computed  the  principal  components  for  the  time  interval  from  1.25  kyr  BP  to 
19.75  kyr  BP  (Fig  5). 

Maps  of  5180Seawater-icevoiume,  averaging  data  over  1,000  years,  and  centered  at 
6.5  kyr  BP,  12.5  kyr  BP,  14.5  kyr  BP,  15.5  kyr  BP,  17.5  kyr  BP,  andl9.5  kyr  BP, 
provide  insight  into  the  spatial  pattern  of  changes  through  time  (Figs  6-9),  discussed 
below. 

3.  Results 

GeoB10069-3 

The  5180Caicite  from  69-3  increased  by  ~0.5%o  from  35  kyr  BP  to  20  kyr  BP, 
consistent  with  the  change  in  global  ice  volume  [Waelbroeck  et  al,  2002].  5180Caidte 
began  to  decrease  at  18  kyr  BP,  dropping  by  ~1.75%o  until  stabilizing  at  ~8kyr  (Fig 
3).  Approximately,  l%o  of  this  decrease  can  be  attributed  to  the  global  change  in  ice 
volume  [Clark  and  Mix,  2002;  Schrag  etal,  2002]. 

The  Mg/Ca  varies  between  3.3  and  3.6  mmol/mol  (or  ~25°C  and  26°C)  from 
35-20kyr  BP.  Mg/Ca-based  temperatures  from  20kyr  BP  until  present  were 
previously  published  [Gibbons  et  al.,  in  revision],  but  are  described  again  briefly 
here.  Mg/Ca  values  increase  to  3.8  mmol/mol  (~26.5°C)  from  20kyr  until  15kyr  BP. 
Mg/Ca  values  then  remain  constant  until  ~12kyr  BP.  Mg/Ca  values  peak  at  10  kyr 
BP  (4.5  mmol/mol,  28.5°C),  and  then  decrease  by  ~0.5  mmol/mol  (1.25°C)  by  the 
late  Holocene. 
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The  5180Seawater-icevoiume  estimates  suggest  several  notable  features.  The  glacial 
(35kyr-20kyr  BP)  portion  averages  0.1%o.  There  are  several  enrichment  events 
during  this  portion  of  the  record.  The  largest  is  an  ~0.2%o  enrichment  from  32.5- 
31kyr  BP,  approximately  coincident  with  some  estimates  of  the  age  of  Heinrich 
Event  3  [Hemming,  2004],  Another  small  enrichment  event  occurred  at  ~  28kyr  BP. 
The  most  pronounced  events,  however,  occurred  during  the  deglacial,  from  19.5kyr 
BP  until  15kyr  BP,  and  from  ~12.5kyr  until  11.2kyr  BP.  During  the  Holocene, 
5180Seawater-icevoiume  values  decreased  from  0%o  in  the  early  Holocene  to 
approximately  -0.2%o  by  5kyr  BP.  5180Seawater-icevoiume  values  then  remained 
approximately  constant  for  the  remainder  of  the  Holocene. 

Regional  data 

The  records  that  we  have  compiled  from  the  IPWP  and  Eastern  Equatorial 
Pacific  Ocean  show  a  wide  range  of  deglacial  temperature  histories  (Fig  4b). 
Although  5 18 Oseawater-icevoiume  trends  also  vary,  most  records  show  generally  enriched 
518Oseawater-icevoiume  values  during  the  deglacial  (Fig  4c).  The  relationship  between  the 
518 Oseawater-icevoiume  records  was  explored  using  principal  component  analysis.  The 
first  principal  component  (PCI;  Fig  5)  explains  53%  of  the  variance  in  the  records. 
There  is  a  longitudinal,  though  not  latitudinal,  coherence  to  how  strongly  the 
records  contribute  to  the  first  principal  component  (eigenvector  elements  or 
loadings  are  shown  in  Fig  6).  Cores  in  the  IndoPacific  Warm  Pool  and  Eastern 
Equatorial  Pacific  have  high  loadings  and  contribute  approximately  equally.  The 
cores  with  the  smallest  loadings  for  the  first  principal  component  are  the  northern 
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most  cores  (MD41  and  MD90)  and  MD78,  MD70,  and  MD76.  We  will  refer  to  the 
latter  group  (MD78,  MD70,  and  MD67)  as  the  "southern  cores.”  Strictly  speaking 
they  are  not  the  southernmost  cores  (see  Table  1),  but  they  appear  to  group 
together  hydrologically.  The  second  and  third  principal  components  explain  9%  and 
8%,  respectively,  and  there  is  no  regional  coherence  to  the  loadings,  nor  does  their 
temporal  pattern  appear  meaningful  (not  shown).  We  thus  infer  that  these  principal 
components  are  dominated  by  noise  rather  than  real  oceanographic  variability. 

The  arithmetic  average  of  5180Seawater-icevoiume  from  all  records  results  in  a 
stacked  record  that  contains  all  the  same  major  features  as  the  first  principal 
component,  providing  confidence  that  both  represent  important  basin-wide 
changes.  Average  518Oseawater-icewater  values  were  similar  in  the  glacial  and  late 
Holocene,  with  values  of  approximately  0%o.  As  with  PCI,  there  were  two  large 
enrichment  events,  of  over  0.3%o  from  17. 75-15. 75yr  BP  and  at  12.25  kyr  BP  (Fig 
4). 

4.  Discussion 

Deglacial  Changes  in  hydrology 

Throughout  most  of  the  IndoPacific  Warm  Pool  and  Eastern  Equatorial 
Pacific,  the  5180Seawater-icevoiume  was  similar  at  the  LGM  and  Holocene  with  an 
approximately  0.3%o  enrichment  from  17.75-12kyr  BP.  This  near  basin-scale 
coherency  is  surprising  in  light  of  the  large  range  in  deglacial  temperature  patterns 
(Fig  4b).  Both  PCI  and  the  average  of  the  5180Seawater-icevoiume  records  reveal  that  this 
enriched  period  can  be  divided  into  two  intervals,  which  correspond  to  HI  and  the 
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Younger  Dryas.  We  note  the  enriched  period  corresponding  to  HI  is  much  longer 
than  the  traditional  HI  time  interval  (~17.5-16  kyr  BP)  [Alley  and  Clark,  1999],  and 
spans  the  full  HI  stadial  or  extended  HI  interval  from  ~19  to  14.6  kyr  BP  [Hall  etal, 
2006;  McManus  etal,  2004;  Stanford  etal.,  2006;  Stanford  etal.,  2011]. 

The  first  principal  component  further  suggests  that  this  variability  is  a 
shared  feature  among  most  of  the  cores.  The  loadings  from  PCI  help  us  explore  the 
regional  coherence  in  the  5180Seawater-icevoiume  pattern  (Fig  6).  Cores  from  the  Eastern 
Equatorial  Pacific  and  IndoPacific  contribute  similarly  to  PCI.  The  southernmost 
cores  (MD70,  MD76,  and  MD78)  in  our  compilation  contribute  the  least.  In  fact,  not 
only  are  these  cores  minor  contributors  to  the  first  PC,  their  5180Seawater-icevoiume 
records  show  relatively  little  variation  throughout  the  last  20  kyr.  The 
northernmost  cores  (MD41  and  MD90)  are  also  minor  contributors  to  PCI,  though 
in  both  cores,  there  are  slightly  enriched  values  corresponding  to  HI  and  the 
Younger  Dryas,  respectively.  The  reason  that  MD41  does  not  contribute  to  the 
loading  of  PCI  is  likely  due  to  local  factors  at  MD41  (discussed  below)  as  well  as  a 
short  HI  enrichment.  At  MD90,  the  5180Seawater-icevoiume  increases  throughout  the 
Holocene,  which  has  been  attributed  to  increased  freshwater  (decreased  salt) 
export  from  the  South  China  Sea  due  to  the  flooding  of  the  Sunda  Shelf  [Linsley  et  al, 
2010].  Additionally,  MD90  has  a  relatively  muted  expression  of  the  Younger  Dryas 
and  HI.  The  presence  of  a  latitudinal  gradient  and  lack  of  a  longitudinal  gradient 
provides  additional  insights  on  possible  mechanisms  for  the  dominant  hydrologic 
changes  in  the  tropical  Indo-Pacific  during  the  deglacial,  captured  by  518Oseawater- 

icevolume- 
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In  the  modern,  climate  variability  described  by  El  Nino-Southern  Oscillation 
(ENSO)  results  in  large  precipitation  anomalies  [e.g.  Wang  and  Schimel,  2003]  and 
several  authors  have  attributed  past  changes  in  518Oseawater  to  ENSO  or  changes  in 
Pacific  mean  state  variability  [Lea  et  al,  2000;  Levi  et  al,  2007;  Rosenthal  et  al, 

2003;  Stott  et  al,  2002].  For  example,  the  enriched  5180  seawater  values  in  MD81 
during  high  latitude  cool  events  were  attributed  to  "El  Nino-like"  conditions  [Stott  et 
al.,  2002].  Though  earlier  work  had  suggested  a  reduced  zonal  salinity  gradient  in 
the  Pacific  during  cold  intervals  [Lea  etal,  2000;  Rosenthal  et  al.,  2003],  more  recent 
work  using  an  Eastern  Equatorial  Pacific  518Oseawater  record  derived  from  alkenone- 
based  temperature  and  S18Ocaicite,  suggests  synchronous  deglacial  5180  seawater 
changes  in  the  eastern  and  western  Pacific,  arguing  against  ENSO-like  changes 
[Leduc  etal.,  2009].  Our  results  expand  on  this  initial  observation  using  3  cores  from 
the  Eastern  Equatorial  Pacific  and  11  cores  from  the  IndoPacific  Warm  Pool. 
Moreover,  we  focus  only  on  5180  seawater  reconstructions  based  entirely  on 
foraminifera,  which  helps  to  reduce  age  uncertainty. 

Indeed,  the  data  we  have  compiled  suggest  that  the  East-West  5180  seawater 
gradient  remained  approximately  constant  through  time  (Figs  7-9).  As  ENSO 
variability  is  characterized  by  opposite  precipitation  anomalies  in  the  eastern  and 
western  tropical  Pacific,  the  data  are  not  consistent  with  ENSO-type  changes. 
Moreover,  speleothem  records  from  two  locations  in  Brazil,  which  today  experience 
opposite  rainfall  anomalies  associated  with  ENSO,  show  synchronous  millennial- 
scale  changes  in  precipitation  during  the  deglacial,  also  suggesting  no  mean  ENSO 
change  from  the  LGM  to  present  [Wang  etal,  2007].  It  has  been  hypothesized  that  if 
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the  tropical  Pacific  becomes  "locked”  in  either  an  El  Nino  or  La  Nina  state,  the 
resulting  teleconnections  could  control  global  climate  [Cane,  1998;  Clement  et  ah, 
1999],  Our  compilation,  however,  does  not  provide  any  support  for  the  hypothesis 
that  the  tropical  Pacific  is  the  source  of  high  latitude  millennial-scale  climate  change, 
though  we  cannot  rule  out  that  there  are  some  changes  in  the  Walker  Cell  that  our 
data  do  not  resolve.  Currents  can  advect  and  mute  518Oseawater  signals,  and  vapor 
transport  over  the  Isthmus  of  Panama  can  change  with  the  migration  of  the  ITCZ, 
potentially  masking  a  change  in  the  Walker  Cell  [LeGrande  and  Schmidt,  2011;  Oppo 
etal,  2007]. 

The  variations  in  this  dataset  are  most  easily  explained  by  changes  in  the 
position  of  the  Intertropical  Convergence  Zone  (ITCZ).  We  first  examine  the  LGM  to 
Holocene  difference.  Most  records  suggest  there  was  relatively  little  518Oseawater- 
icvoiume  change  between  the  LGM  and  the  Holocene.  A  notable  difference  between  the 
LGM  and  Holocene  in  one  record  can  be  attributed  to  local  factors;  the  relatively 
depleted  5180Seawater-icevoiume  conditions  in  the  Sulu  Sea  (MD41)  may  be  due  the  lower 
sea  level  at  the  LGM  which  allowed  relatively  more  fresh  South  China  Sea  water  than 
equatorial  Pacific  water  to  enter  the  Sulu  Sea,  decreasing  its  5180  seawater  values 
[Rosenthal et al,  2003].  Speleothem  records  from  China  [Yuan  etal.,  2004]  and 
Borneo  [Partin  etal.,  2007]  may  suggest  a  decrease  in  precipitation  during  the  LGM. 
However,  as  noted  by  Partin  [2007]  once  global  518Oseawater  changes  are  removed, 
the  speleothem  5180  changes  are  smaller,  suggesting  a  modest  decrease  in  the 
monsoons  during  the  LGM,  which  can  be  at  least  in  part  attributed  to  cooler  SST. 
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Thus  the  enriched  5180  values  in  Borneo  are  not  inconsistent  with  5180Seawater-icevoiume 


evidence  suggesting  the  ITCZ  position  during  the  LGM  was  similar  to  the  Holocene. 

Though  there  is  relatively  little  5180Seawater-icevoiume  change  from  the  LGM  to 
modern,  the  data  suggest  dramatic  hydrologic  changes  in  parallel  with  Northern 
Hemisphere  high  latitude  rapid  temperature  change.  Others  have  suggested  a 
southward  migration  of  the  ITCZ  during  the  Younger  Dryas  and  HI  [Benway  et  al., 
2006;  Leduc  et  al.,  2009],  and  the  greater  spatial  coverage  of  the  cores  we  have 
compiled  allows  us  to  evaluate  this  hypothesis  in  greater  detail.  The  two  northern 
most  cores  (MD41  and  MD90)  show  some  enrichment  in  the  5180Seawater-icevoiume 
values,  the  equatorial  cores  (MD62,  MD81,  MD41,  43JC,  V21090,  TR163-22,  69-3, 
MD65,  29-4,  and  38-4)  show  larger  5180Seawater-icevoiume  enrichments,  and  the 
southern  cores  (MD78,  MD76,  and  MD70)  show  no  change.  To  the  south  of  our  study 
area,  wetter  conditions  are  documented  in  a  peat  bog  during  HI  (though  not  during 
the  Younger  Dryas)  [Muller  et  al.,  2008].  As  discussed  in  more  detail  below,  this 
spatial  pattern  is  consistent  with  a  southward  displacement  of  the  ITCZ. 

Enriched  5180  values  in  speleothems  from  Borneo  during  HI  have  been 
interpreted  as  a  decrease  in  local  rainfall  caused  by  a  southward  displacement  of  the 
ITCZ  [Partin  etal,  2007].  Though  Partin  etal.  [2007]  claim  5180  values  are  not 
enriched  during  the  Younger  Dryas,  we  suggest  a  Younger  Dryas  enrichment  is 
masked  by  the  global  decrease  in  5180  seawater.  Removing  the  global  ice  volume- 
related  5180  change  highlights  not  only  the  relatively  small  changes  from  the  LGM  to 
late  Holocene  and  the  relatively  enriched  period  during  HI  (a  1.0%o  enrichment), 
but  also  shows  a  0.5%o  enrichment  during  the  Younger  Dryas.  We  also  note  that  two 


78 


of  the  three  speleothems  that  make  up  the  spliced  Borneo  record  experience  a 
growth  hiatus  coincident  with  the  Younger  Dryas  (13.0  kyBP  -  11.7  kyr  BP  for 
SSCOl  and  12.7  kyr  BP  to  11.9  kyr  BP  for  BA04)  which  is  consistent  with  drier 
conditions.  The  Borneo  record  displays  a  remarkable  correspondence  with  our  PCI 
(Fig  5).  The  difference  between  the  LGM  and  HI  in  the  speleothem  record  at  Borneo 
is  l%o.  This  is  over  twice  as  large  as  is  observed  in  nearby  5180Seawater-icevoiume 
reconstructions  (Fig  4c,  Fig  5)  [Rosenthal  et  al,  2003;  Steinke  etal,  2008;  Stott  etal, 
2004;  Visseretal.,  2003],  confirming  that  the  speleothem  record  reflects  some 
degree  of  decreased  precipitation,  and  not  just  a  change  in  source  water.  Support  for 
a  5180  response  to  reduced  precipitation  in  Borneo  during  North  Atlantic  cold 
events  is  also  provided  by  a  water-isotope-enabled  fully  coupled  model  simulations 
[Lewis  etal,  2010]. 

The  5180Seawater-icevoiume  data  uniformly  suggest  a  southward  migration  of  the 
ITCZ  during  HI  and  the  Younger  Dryas,  but  not  during  the  LGM.  In  the  modern 
ocean,  the  ITCZ  is  biased  towards  the  northern  hemisphere.  This  is  generally 
attributed  to  the  global  continental  configuration,  and  the  resulting  warm  bias  of  the 
northern  hemisphere,  but  the  geometry  of  the  coastlines  and  oceanic/atmospheric 
feedbacks  may  also  play  a  role  [Koutavas  and  Lynch-Stieglitz,  2005].  The  warm  bias 
of  the  Northern  Hemisphere  is  accentuated  by  the  Atlantic  Meridional  Overturning 
Circulation,  which  exports  heat  northward  near  the  surface  at  every  latitude  in  the 
Atlantic  Ocean  [Houghton  etal,  1996].  A  reduction  in  the  Atlantic  Meridional 
Overturning  Circulation  could  decrease  this  heat  transfer,  and  therefore  reduce  the 
interhemispheric  temperature  gradient  [Broecker,  1998], 
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Changes  in  the  interhemispheric  temperature  gradient  associated  with 
millennial  scale  climate  events  consistent  with  the  predictions  of  Broecker  [1998] 
have  been  documented.  Using  a  global  compilation  of  temperature  records  Shakun 
and  Carlson  [2010],  show  that  both  HI  and  the  Younger  Dryas  were  characterized 
by  cooling  in  the  Northern  Hemisphere  and  warming  in  the  Southern  Hemisphere, 
in  contrast  to  the  LGM,  when  cooling  was  global. 

Whether  the  Atlantic  Meridional  Overturning  Circulation  is  actually  the 
source  of  the  millennial  scale  changes  in  the  interhemispheric  temperature  gradient 
is  still  under  debate.  While  many  proxies  suggest  a  change  in  water  mass 
distribution  during  the  LGM  [e.g.  Boyle  and  Keigwin,  1987;  Curry  and  Oppo,  2005; 
Robinson  et  ah,  2005],  it  is  not  obvious  that  this  change  was  associated  with  a 
dramatic  change  in  the  strength  of  the  overturning  circulation  [ Curry  and  Oppo, 
2005;  Hall  etal.,  2006;  McManus  etal,  2004;  Praetorius  etal,  2008;  Stanford  etal., 
2006;  Yu  et  al.,  1996],  and  is  consistent  with  the  lack  of  change  in  interhemispheric 
temperature  gradient.  Thus,  one  possible  explanation  for  the  lack  of  a  LGM  minus 
late  Holocene  5180Seawater-icevoiume  difference  in  our  tropical  study  area  is  that  the  ITCZ 
was  close  to  its  modern  position.  In  addition,  the  5180Seawater-icevoiume  should  not  be 
expected  to  preserve  modest  differences  in  the  5180  of  precipitation,  because 
changes  in  the  ocean  are  muted  due  to  ocean  mixing  and  advection  by  currents  [e.g. 
Oppo  etal,  2007]. 

On  the  other  hand,  several  lines  of  evidence  suggest  a  significant  reduction  in 
the  Atlantic  Meridional  Overturning  Circulation  during  HI  and  the  Younger  Dryas 
[e.g.  Hall  etal,  2006;  McManus  etal,  2004;  Oppo  and  Curry,  2011;  Praetorius  etal, 
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2008;  Stanford  et  al,  2006].  Thus,  the  decrease  in  the  northward  heat  transport 
during  HI  and  the  Younger  Dryas  associated  with  a  reduction  in  the  Atlantic 
Meridional  Overturning  Circulation  may  have  been  more  significant  than  during  the 
LGM. 

Modeling  sensitivity  studies  that  add  freshwater  to  the  North  Atlantic  (e.g. 
"hosing"  experiments),  show  a  decrease  in  the  Atlantic  Meridional  Overturning 
Circulation  and  a  southward  shift  in  the  ITCZ  [Dahl  etal,  2005;  Zhang  and  Delworth, 
2005],  and  specifically  a  drying  in  the  Indo  Pacific  Warm  Pool  region  [Zhang  and 
Delworth,  2005],  A  freshwater  hosing  experiment  using  a  water  isotope-enabled 
fully  coupled  model  allows  us  to  directly  compare  the  618Oseawater  data  with  model 
output  [Lewis  et  al,  2010].  This  model  experiment  simulates  a  southward  shift  in 
the  ITCZ.  The  5180  of  precipitation  is  enriched  in  the  northern  part  of  the  IndoPacific 
and  depleted  in  the  south.  Dividing  the  two  regions  is  a  band  that  experiences  no 
change  in  the  5180  of  precipitation.  According  to  the  model  simulations,  this  region 
of  no  change  experiences  an  increase  in  summer  precipitation  (corresponding  to 
depleted  5180)  and  a  decrease  in  winter  precipitation  (corresponding  to  enriched 
5180).  While  this  band  is  slightly  to  the  north  of  the  cores  that  show  no  change  in  the 
518Oseawater-icevoiume,  we  suggest  that  this  change  in  seasonality  is  the  most  likely 
explanation  for  the  lack  of  5180Seawater-icevoiume  change  in  core  MD70,  MD76,  and 
MD78  during  HI  and  the  Younger  Dryas.  Unfortunately,  5180Seawater-icevoiume  changes 
in  the  entire  tropical  Pacific  and  Indian  Oceans  during  the  simulated  hosing 
experiment  are  not  significant  at  the  95%  level,  so  we  cannot  directly  compare  our 
seawater  reconstruction  to  the  modeled  518Oseawater  response.  However,  the  similar 
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spatial  pattern  to  the  5180  seawater  of  precipitation  lends  support  to  our 
interpretations. 

Conclusion 

Eastern  Equatorial  Pacific  and  IndoPacific  Warm  Pool  changes  in  6180  seawater- 
icevolume  suggest  the  migration  of  the  ITCZ  controls  major  hydrologic  changes  from 
the  LGM  to  present.  The  618Oseawater-icevoiume  values  at  the  sites  we  study  during  the 
LGM  were  similar  to  today’s,  suggesting  relatively  small  changes  in  the  IndoPacific 
tropical  hydrologic  cycle  compared  to  HI  and  the  Younger  Dryas,  when  the 
5180Seawater-icevoiume  pattern  suggests  a  southward  migration  of  the  ITCZ.  The  spatial 
patterns  of  618Oseawater-icevoiume  values  during  HI  and  the  YD  are  consistent  with 
simulated  changes  in  the  5180  of  precipitation  in  a  hosing  experiment.  Utilizing  our 
greater  spatial  resolution,  we  reject  an  ENSO-like  change  from  the  LGM  to  present, 
and  instead  concur  with  the  findings  of  those  previous  studies  which  suggest  that 
there  was  a  southward  ITCZ  migration  during  HI  and  the  Younger  Dryas. 
Additionally,  by  focusing  on  records  that  derive  5180  seawater  entirely  from 
foraminifera,  rather  than  from  two  sediment  constituents  (e.g.  foraminifera  and 
alkenones),  we  increase  our  confidence  in  the  timing  of  these  millennial-scale 
events.  These  results,  combined  with  our  findings  of  an  ITCZ  position  at  the  LGM 
that  is  similar  to  modern,  suggest  that  the  migration  of  the  ITCZ  was  due  to  changes 
in  the  interhemispheric  temperature  gradient. 
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Core  ID,  location 

Core 

ID  in 

text 

Lat./Long. 

Water 

depth 

(m) 

Average 

Sample 

interval 

(yr) 

Mg/Ca  equation 

Reference 

GeoB  10029-4, 

Eastern  Indian  Ocean 

29-4 

r  30's, 

100“  8'E 

962 

600 

Anand  et  al., 

2003 

[Mohtadi  et  al.,  2010] 

GeoB10038-4, 

Eastern  Indian  Ocean 

38-4 

5“  56'S, 

103“  15'E 

1,819 

560 

Anand  et  al., 

2003 

[Mohtadi  et  al.,  2010] 

GeoB  10069-3  Savu 

Sea 

69-3 

9“  36'S, 

120“  55'E 

1,250 

190 

modified  Anand 

this  study 

MD97-2141,  Sulu 

Sea 

MD41 

8“  47'N, 

121“  17'E 

3,633 

90 

SSTs  from 
publication 

[. Rosenthal  et  al.,  2003] 

MD98-2162,  S. 
Makassar  St. 

MD62 

4“  41'S, 

117“  54'E 

1,855 

415 

Anand  et  al., 

2003 

[  Visser  et  al.,  2003] 

MD98-2165,  Sumba 

St. 

MD65 

9“  39'S, 

118“  20'E 

2,100 

200 

Anand  et  al., 

2003 

[Levi  et  al.,  2007] 

MD98-2170,  Timor 

Sea 

MD70 

10“  36'S, 
125“23'E 

832 

300 

Anand  et  al., 

2003 

[Stott  et  al.,  2007] 

MD98-2176,  Seram 

Sea 

MD76 

5“  00'S, 

133“  27'E 

2,382 

75 

Anand  et  al., 

2003 

[Stott  et  al.,  2007] 

MD98-2181,  W. 
Pacific,  Mindanao 

MD81 

6“  27'N, 

125“  50'E 

2,114 

50 

Anand  et  al., 

2003 

[Stott  et  al.,  2007] 

MD01-2378,  Timor 

Sea 

MD78 

13“  05'S, 
121°  47'E 

1,783 

125 

modified  Anand 

[Sarnthein  et  al.,  201 1; 

Xu  et  al.,  2008] 

MD01-2390,  South 
China  Sea 

MD90 

6“  38'N, 

113°  25'E 

1,545 

200 

Anand  et  al., 

2003 

[Steinke  et  al.,  2008] 

ME0005A-43JC 
Eastern  Eq.  Pacific 

43JC 

7“  51'N, 

83“  37'W 

1,368 

240 

SSTs  from 
publication 

[Benway  et  al.,  2006] 

TR1 63-22  Eastern 

Eq.  Pacific 

TR163- 

22 

0“  31'N, 

92“  24'W 

2,830 

270 

SSTs  from 
publication 

[Lea  et  al. ,  2006] 

V21-30  Eastern  Eq. 
Pacific 

V21-30 

1“  13'S, 

89“  41'E 

617 

430 

Anand  et  al., 

2003 

[Koutavas  et  al.,  2002] 

Table  1:  Cores  used  in  this  study.  All  cores  use  data  from  G.  ruber,  except  V21-30 
which  uses  G.  sacculifer. 
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Fig  1: 

Core  locations  and  mean  annual  sea  surface  salinity  from  the  World  Ocean  Atlas  05 
[Antonov  etal.,  2006]. 
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Figure  2:  Foraminifera-based  reconstructed  "modern"  618Oseawater  estimates  plotted 
against  modern  WOA05  salinity  data  [Antonov  etal.,  2006],  Modern  water 
5180Seawater  and  salinity  data  [LeGrande  and  Schmidt,  2006]  are  also  shown.  The 
IndoPacific  plot  includes  all  water  data  from  90°E-140°E  and  10°N-10°S.  The 
Eastern  Equatorial  Pacific  plot  includes  data  from  90°W  to  100°W  and  11°N  to  10°S. 
All  water  samples  are  from  100m  or  above  and  most  are  from  the  surface.  For  the 
IndoPacific  samples,  the  slope  for  the  foraminifera-based  line  is  0.48  (R2  =  0.76)  and 
for  the  water  samples  it  is  0.44  (R2  =  0.79).  For  the  Eastern  Equatorial  Pacific  the 
slope  for  the  modern  is  0.22  (R2  =  0.70)  and  the  reconstructed  518Oseawater  values 
tend  to  fall  within  the  scatter  of  the  data. 


36 


85 


-3.5  — 


0  10  20  30 

Thousands  of  years  BP 


Figure  3:  Data  from  69-3.  Top  panel  is  5180  of  calcite.  Next  panel  is  the  Mg/Ca-based 
temperature.  These  records  are  combined  below  [Bemis  etal.,  1998]  to  derive  the 
5180Seawater.  Finally,  we  show  the  5180  seawater  with  the  global  ice  volume  effect 
removed.  Northern  hemisphere  cold  events  (the  Younger  Dryas  and  Heinrich  Stadial 
1),  which  correspond  to  periods  of  enriched  618Oseawater  are  indicated  by  blue  boxes. 
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Figure  4a:  518Ocaicite  data  from  all  cores  used  in  this  study.  The  cores  are  grouped  by 
region. 
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Figure  4b:  Mg/Ca-based  temperature  data  from  all  cores  used  in  this  study.  Note  the 
wide  range  of  deglacial  temperature  histories. 
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Figure  4c:  5180Seawater-icevoiume  reconstructions  from  all  cores  used  in  this  study.  Note 


the  generally  similar  LGM  and  late  Holocene  values  for  the  618Oseawater-icevoiume.  With 
the  exception  of  the  southern  cores  (bottom  panel),  many  of  the  records  show 
enrichment  events  during  HI  and  the  Younger  Dryas. 
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Figure  5:  A  comparison  of  climate  records.  Top  panel  is  a  5180  records  from  Dome 
Fuji  [Kawamura  etal.,  2007]in  the  Indian  Ocean  sector  of  Antarctica.  Second  panel  is 
GISP  5180  from  Greenland  [Grootes  and  Stuiver,  1997].  The  third  panel  is  the  first 
principal  component  (PCI)  of  all  IndoPacific  and  Eastern  Equatorial  Pacific 
518Oseawater-icevoiume  records  included  in  this  study  (see  table  1  for  included  cores). 
Dots  (blue)  are  an  average  of  all  the  5180Seawater-icevoiume  records.  Note  the  inverted 
axes  for  both  PCI  and  the  average  5180Seawater-icevoiume  record.  The  final  panel  is  a 
spliced  record  of  5180  from  three  Borneo  speleothems  [Partin  et  ah,  2007],  with  the 
global  change  in  S180  due  to  ice  volume  removed  [Clark  and  Mix,  2002;  Schrag  et  al., 
2002], 
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Figure  6:  Top  panel  is  eigenvector  elements  (loadings)  for  PCI.  Following  panel  is 
S180Seawater-icevoiume  anomalies  relative  to  present.  For  MD41,  the  modern  518Oseawater 
estimate  from  the  gridded  data  of  Le  Grade  and  Schmidt  [2006]  was  used  to 
calculate  the  anomaly.  For  all  other  cores  anomalies  were  calculated  with  respect  to 
the  last  2,000  years.  Bottom  panel  is  an  average  of  1,000  years  and  is  centered  on 
6.5  kyr  BP 
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Figure  7:  5180Seawater-icevoiume  anomalies  as  described  for  Figure  6.  Top  panel  is 
centered  on  12.5kyr  BP  and  bottom  panel  is  centered  on  14.5  kyr  BP.  There  are 
generally  enriched  conditions  during  the  Younger  Dryas  (12.5  kyr  BP)  and  a  partial 
recovery  to  late  Holocene  conditions  during  the  Bplling-Allerpd  (14.5  kyr  BP).  Note 
the  lack  of  change  in  the  East-West  gradient. 
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Figure  8:  5180Seawater-icevoiume  anomalies  as  described  for  Figure  6.  Top  panel  is 
centered  on  15.5  kyr  BP  and  bottom  panel  is  centered  on  17.5  kyr  BP.  Enriched 
conditions  occur  throughout  Heinrich  Stadial  1,  and  are  not  only  confined  to 
Heinrich  Event  1. 
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Figure  9:  518Oseawater-icevoiume  anomalies  as  described  for  Figure  6.  Panel  is  centered 
on  19.5kyr  BP  (the  LGM).  Relatively  small  anomalies  suggest  the  position  of  the  ITCZ 
at  the  LGM  was  very  similar  to  the  Late  Holocene. 
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Chapter  4 


The  thermocline  structure  of  the  IndoPacific  Warm  Pool  from  the  Last  Glacial 

Abstract 

The  thermocline  structure  of  the  IndoPacific  Warm  Pool  is  influenced  by  both 
global  and  local  processes.  Several  studies  suggest  that  it  responds  to  changes  in  the 
Walker  Circulation,  and  others  suggest  that  it  is  sensitive  to  the  inflow  of  fresh 
water  from  the  South  China  Sea.  However,  observations  and  paleoceanographic 
reconstructions  are  sparse  throughout  the  IndoPacific  Warm  Pool,  and  there  is  still 
considerable  uncertainty  related  to  the  influence  of  the  atmospheric  circulation 
variability  and  the  South  China  Sea  inflow  on  the  thermocline  structure,  especially  in 
the  Indonesian  Throughflow  region.  Here  we  present  thermocline  reconstructions 
using  the  5180  of  calcite,  Mg/Ca-based  temperature  estimates,  and  5180  of  seawater 
estimates  using  the  mixed  layer  foraminifera  Globigerinoides  ruber  ( G .  ruber )  and 
the  thermocline-dwelling  foraminifera  PuIIeniatina  obliquiloculata  (P. 
obliquiloculata ).  We  present  data  from  both  coretops  and  a  suite  of  downcore 
records  from  throughout  the  IndoPacific  Warm  Pool.  Our  results  suggest  that  in  the 
modern  ocean,  the  Mg/Ca-based  temperature  of  P.  obliquiloculata  is  a  valid  proxy 
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for  upper  thermocline  temperature,  whereas  some  caution  is  required  in  its  use 
downcore.  If  the  Mg/ Ca-based  thermocline  temperatures  for  the  Last  Glacial 
Maximum  (LGM)  are  reliable,  our  data  suggest  the  temperature  difference  between 
surface  and  subsurface  waters  in  IndoPacific  Warm  Pool  was  uniformly  smaller  at 
the  LGM,  perhaps  consistent  with  a  more  vigorous  Walker  Circulation  [DiNezio  et  ah, 
2011].  Further  our  suite  of  data  argue  against  the  hypothesis  that  the  flooding  of  the 
Sunda  Shelf,  which  allowed  South  China  Sea  water  to  enter  the  Indonesian  Seas, 
caused  the  cool  thermocline  temperatures  observed  in  the  outflow  passages  of  the 
Indonesian  Throughflow  after  ~10kyr  before  present  (BP).  We  note,  however,  there 
is  difficulty  in  reconciling  the  5180Caidte  and  Mg/Ca-based  temperature  data  from  P. 
obliquiloculata,  and  suggest  further  calibration  work  is  required. 

1.  Introduction 

The  Indonesian  Throughflow  (ITF)  connects  the  Pacific  to  the  Indian  Ocean 
through  a  series  of  passages  between  Indonesian  islands.  It  forms  the  only  low- 
latitude,  inter-ocean  connection,  allowing  direct  heat  and  salinity  exchanges 
between  the  Pacific  and  Indian  Oceans  [Gordon,  2005].  ITF  water  has  been 
identified  in  Aghulas  eddies  within  the  South  Atlantic  Ocean,  and  it  is  believed  to 
contribute  to  the  warm  water  return  limb  of  the  global  overturning  circulation 
[Gordon,  1986],  A  relatively  higher  sea  level  in  the  western  equatorial  Pacific  Ocean 
drives  the  ITF  into  the  Indian  Ocean  [Wyrtki,  1987].  Approximately  85%  (8  Sv)  of 
the  water  transported  through  the  ITF  enters  from  the  Sulawesi  Sea  and  passes 
through  the  Makassar  Strait  [Gordon  et  al,  2003].  From  there,  most  flows  into  the 
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Flores  Sea,  and  then  to  the  Banda  Sea  (Figure  1).  Within  the  Banda  Sea,  vigorous 
mixing  with  local  freshwater  creates  characteristic  low-salinity  ITF  water  [Fieux  et 
al,  1994;  Gordon,  2005].  There  are  three  main  exits  for  the  ITF:  the  Lombok  Strait, 
the  Ombai  Strait,  and  the  Timor  Sea,  and  they  respectively  account  for 
approximately  20%,  30%,  and  50%  of  the  total  ITF  outflow  [Sprintall  etal.,  2009]. 

Estimates  of  the  total  ITF  transport  vary  widely  due  to  sparse  direct 
observations,  large  seasonal  and  interannual  variability,  and  the  difficulty  of 
modeling  the  region  [reviewed  in  Godfrey,  1996].  The  most  recent  measurements 
suggest  that  the  transport  is  between  9  and  15  Sv  [Gordon,  2005;  Sprintall  et  al, 
2009]  with  the  monsoon  [Gordon  etal,  2003;  Sprintall  et  al,  2009],  El  Nino- 
Southern  Oscillation  (ENSO)  [Ffield  etal,  2000;  Gordon  etal,  2008],  and  the  Indian 
Ocean  Dipole  (IOD)  [Sprintall  etal,  2009]  all  contributing  to  the  observed 
variability. 

During  the  boreal  winter  monsoon  (all  seasons  are  discussed  with  respect  to 
the  Northern  Hemisphere),  fresh  surface  water  from  the  South  China  Sea  is  advected 
through  the  ]ava  Sea  into  the  southern  Makassar  Strait  and  Flores  Sea  (dashed 
arrow  in  Figure  1)  [Fang  et  al,  2010].  This  buoyant  fresh  water  reduces  the  sea 
surface  salinity  of  the  Makassar  Strait  in  relation  to  the  Western  Pacific  Warm  Pool. 
Previous  work  has  suggested  that  this  "freshwater  plug"  reduces  the  density 
gradient  along  the  strait,  inhibits  the  flow  of  surface  water  through  the  Makassar 
Strait,  and  intensifies  the  thermocline  flow  [Gordon  etal,  2008].  While  this 
mechanism  may  contribute  to  the  thermocline  dominated  flow  observed  in  the 
Makassar  Strait,  it  cannot  be  the  only  explanation.  The  "freshwater  plug"  is  only 
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present  during  the  winter  monsoon,  but  the  flow  is  greater  in  the  thermocline  than 
in  the  surface  in  the  Makassar  Strait  throughout  the  year  [Gordon  et  al,  2008],  A 
complementary,  or  perhaps  alternative,  mechanism  involves  the  persistently  higher 
sea  level  in  the  southern  Java  Sea  compared  to  the  Flores  Sea.  The  pressure  gradient 
that  results  from  this  should  be  associated  with  a  northward  geostrophic  flow  in  the 
Makassar  Strait  surface  layer,  but  not  in  the  thermocline,  which  is  below  the  sill 
depth  of  Java  Sea  [Gibbons  et  al.,  in  revision] .  Regardless  of  mechanism,  the  South 
China  Sea  inflow  has  a  profound  affect  on  the  net  volume  and  heat  transport  of  the 
Makassar  Strait,  and  contributes  to  the  fresh  nature  of  the  ITF  [Gordon,  2005]. 

The  extent  to  which  this  thermocline  dominated  flow  in  the  Makassar  Strait 
contributes  to  characteristics  of  the  ITF  is  a  topic  of  investigation.  It  has  been 
proposed  that  the  thermocline  dominated  flow  of  the  Makassar  Strait  reduces  the 
transport-averaged  temperature  of  the  Makassar  Strait  by  2°C  [Tozuka  et  al.,  2007] 
and  has  the  net  effect  of  cooling  the  Indian  Ocean  [Gordon,  2005].  Evidence  that  the 
ITF  cools  comes  from  the  subsurface  temperature  of  the  Indian  Ocean.  The  ITF  can 
be  identified  as  a  cool  tongue  at  ~200m  depth  that  stretches  across  the  equatorial 
Indian  Ocean  [Gordon,  2005].  However,  it  has  been  shown  recently  that  surface 
water  from  the  Java  Sea  is  a  source  of  heat  to  the  ITF  [Fang  etal.,  2010],  which  likely 
compensates  for  the  decreased  surface  heat  transport  through  the  Makassar  Strait. 
The  apparently  cool  nature  of  the  ITF  is  more  likely  due  to  the  intense  upwelling  in 
the  Banda  Sea,  which  brings  cool  deep  water  to  the  surface  [Gibbons  et  al,  in 
revision;  Koch-Larrouy  et  al,  2008]. 
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Much  of  the  shallow  Java  Sea  was  exposed  during  the  Last  Glacial  Maximum 
(LGM),  which  prevented  South  China  Sea  waters  from  entering  the  ITF.  The  early 
Holocene,  therefore,  represents  a  unique  opportunity  to  evaluate  whether  and  how 
the  inflow  from  the  South  China  Sea  influences  the  ITF,  and  in  turn,  the  influence  of 
the  ITF  on  the  region’s  salinity  and  temperature  distribution.  A  recent  study 
documented  a  relative  freshening  (lighter  5180  of  seawater)  in  the  surface  of  the 
southern  Makassar  Strait  relative  to  the  Western  Pacific  Warm  Pool  when  the 
connection  between  the  South  China  Sea  and  the  Makassar  Strait  was  established  at 
~9-10  kyrBP  [Linsley  etal,  2010].  Initial  work  suggested  that  the  flooding  of  the 
Sunda  Shelf  did  not  affect  regional  sea  surface  temperature  (SST)  gradients  [Xu  et 
ah,  2008;  Xu  et  al,  2010].  More  recent  work  has  shown  a  distinct  warming  in  the 
Makassar  Strait  at  the  time  of  the  flooding  of  the  Sunda  Shelf,  but  not  in  the  outflow 
passages,  suggesting  that  the  opening  of  the  Java  Sea  between  the  South  China  Sea 
and  the  ITF  inhibited  surface  flow  through  the  Makassar  Strait,  trapping  heat  that 
was  previously  exported  [Gibbons  et  al,  in  revision].  Furthermore,  it  has  been 
suggested  that  the  flooding  of  the  Sunda  Shelf  cooled  thermocline  temperatures  and 
intensified  thermocline  flow  [Xu  etal,  2008],  a  hypothesis  we  examine  in  this  study. 

Here,  we  present  5180Caidte,  Mg/Ca-based  temperature,  and  5180  seawater  ofP. 
obliquiloculata  from  coretops  and  downcore  records  throughout  the  IndoPacific 
Warm  Pool.  These  new  records,  in  conjunction  with  published  data  from  the  mixed 
layer  dweller  G.  ruber  from  the  same  cores  [Gibbons  et  al,  in  revision;  Linsley  et  al, 
2010;  Rosenthal  et  al,  2003],  and  published  surface  and  thermocline  records  from  a 
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single  Timor  Sea  core  [Xu  etal,  2008]  (Table  1)  allow  us  to  examine  regional 
changes  in  thermocline  structure  that  occurred  throughout  the  last  deglaciation. 

2.  Materials  and  Methods 

All  samples  for  the  coretop  study  were  taken  with  a  multi-coring  device,  on 
three  cruises  in  the  Indonesian  Throughflow  region  and  the  Eastern  Indian  Ocean 
(Figure  1,  dots).  We  develop  downcore  P.  obliquiloculata  records  from  the  Makassar 
Strait  (BJ8-03  70GGC  and  23GGC),  the  Ombai  Strait  (GeoB10069-3),  the  Flores  Sea 
(BJ8-03  136GGC),  and  the  Sulu  Sea  (MD97-2141).  These  core  locations  are 
indicated  by  diamonds  in  Figure  1. 

The  core  sites  are  influenced  to  varying  degree  by  monsoon  winds.  SSTs  in 
the  Makassar  Strait,  Flores  Sea,  and  Savu  Basin  are  highest  in  winter  and  coldest  in 
summer,  reflecting  the  influence  of  the  monsoon  winds  on  upwelling.  Maximum 
SSTs  are  ~31°C  at  70GGC  on  the  Sulawesi  margin,  ~31°C  at  69-3  in  the  Savu  Basin, 
and  ~1°C  cooler  at  23GGC  in  the  south  central  Makassar  Strait  and  136GGC  in  the 
Flores  Sea.  Minimum  SSTs  are  ~28°C  at  70GGC,  27°C  at  23GGC  and  136GGC,  and 
26°C  at  69-3.  The  mixed  layer  is  ~50m  deep  in  the  Makassar  Strait  (23GGC  and 
70GGC)  and  Flores  Sea  (136GGC)  [Wyrtki,  1961].  The  Savu  Sea  (69-3)  has  a  very 
steep  thermocline  and  a  mixed  layer  that  is  no  more  than  20m  deep  [Molcard  et  ah, 
2001;  Wyrtki,  1961].  MD41,  in  the  Sulu  Sea,  is  the  only  sediment  core  in  this  study 
from  the  northern  hemisphere.  SSTs  at  MD41  peak  during  summer  (~30°C)  and 
reach  a  minimum  of  ~27°C  during  the  winter.  The  mixed  layer  is  ~30  m  deep 
[Wyrtki,  1961]. 
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Approximately,  20  P.  obliquiloculata  (355-425pm)  were  used  for  Mg/Ca 
analysis.  Samples  were  gently  crushed  prior  to  cleaning  in  order  to  open  the 
chambers.  Approximately  40-70  pg  was  set  aside  for  stable  isotope  analysis 
(described  below).  Samples  were  cleaned  using  a  full  trace  metal  cleaning  method, 
including  oxidative  and  reductive  cleaning  steps  [modified  from  Boyle  and  Keigwin, 
1985/6;  Rosenthal  etal,  1999]  and  run  on  an  Element  XR ICP-MS  at  Rutgers 
University  or  an  Element  2  ICP-MS  at  the  Woods  Hole  Oceanographic  Institution 
(WHOI).  To  construct  downcore  records,  Mg/Ca  data  were  converted  to 
temperatures  using  the  Anand  [2003]  multispecies  equation,  modified  to  account 
for  cleaning  technique  [Gibbons  et  ai,  in  revision]: 

Anand  multispecies  equation:  Mg/Ca  =  0.38exp(0.090T) 

Modified  Anand  equation:  Mg/Ca  =  0.345exp(0.090T) 

We  note  that  there  is  a  wide  discrepancy  in  temperature  reconstructions 

depending  on  the  temperature  calibration  used.  Others  have  suggested  that  the 
equation  Mg/Ca  =0.78xp(0.052T)  [Cleroux  et  al,  2008;  Mohtadi  etal.,  2011a]  is 
more  appropriate  for  P.  obliquiloculata.  This  issue  is  discussed  below. 

Stable  isotope  measurements  were  made  at  WHOI  on  a  Finnigan  MAT253. 
Long-term  precision  of  5180Caidte  measurements  is  0.07%o,  based  on  repeated 
measurements  of  a  suite  of  standards. 

All  age  models  are  based  on  radiocarbon  dates  from  planktic  foraminifera. 
Radiocarbon  years  were  covered  to  calendar  ages  using  Calib  6.0  and  Marine09 
[Reimer  et  al,  2009].  Age  models  are  based  on  linearly  interpreting  between 
radiocarbon  dates  (23GGC,  70GGC,  136GGC,  and  69-3),  a  linear  fit  through  the 
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radiocarbon  dates  (MD78),  or  a  polynomial  fit  (MD41).  Details  of  the  age  model 
construction  and  original  dates  for  23GGC,  136GGC,  70GGC  and  69-3  can  be  found  in 
Gibbons  etal.  [in  revision].  See  Sarnthein  etal.  [2011]  (MD78)  and  de  Garidel- 
Thoron  et  al.  [2001]  (MD41)  for  the  rest  of  the  radiocarbon  dates. 

We  explored  the  potential  for  an  age  bias  between  radiocarbon  dates  from 
surface  dwelling  foraminifera  (a  mix  of  G.  ruber  and  G.  sacculifer )  and  P. 
obliquiloculata,  as  there  has  been  a  documented  age  offset  between  dissolution 
prone  [G.  sacculifer )  and  dissolution-resistant  (such  as  P.  obliquiloculata  or  N. 
dutertrei )  species  [Barker  et  al.,  2007;  Broecker  and  Clark,  2011],  We  selected 
several  depth  intervals  from  70GGC  and  69-3  to  compare  the  radiocarbon  ages  of 
these  species.  In  all  cases  we  examined,  the  offset  between  species  was  small  (Table 
2).  These  results  support  the  conclusion  of  Broecker  and  Clark  [2011]  that  there  is 
minimal  age  offset  in  high  sedimentation  rate  cores.  We  therefore  applied  no  offset 
for  our  age  models  and  used  only  the  radiocarbon  values  that  were  measured  on  G. 
ruber  and  G.  sacculifer,  and  not  on  P.  obliquiloculata. 

3.  Results 

3.1  Coretop  results 

We  reconstructed  the  Mg/Ca-based  temperature  [Anand  etal,  2003;  Gibbons 
et  al.,  in  revision]  of  P.  obliquiloculata  from  22  coretops  from  throughout  the  ITF 
region  and  the  Eastern  Indian  Ocean.  We  estimated  depth  habitats  by  comparing 
our  Mg/Ca  based  temperature  estimates  with  the  CTD  data  collected  on  the  coring 
cruises.  These  cruises  were  all  during  the  upwelling  season,  and  the  CTD  data  only 
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provide  us  with  a  snapshot  of  subsurface  temperature  conditions.  However,  the 
paucity  of  subsurface  data  in  the  World  Ocean  Atlas  in  the  ITF  and  the  Eastern 
Indian  Ocean  could  make  CTD  data  the  most  reliable  estimate  of  subsurface 
temperature.  The  depth  habitat  estimates  suggest  that  P.  obliquiloculata  calcifies 
within  the  upper  thermocline,  at  approximately  80-100  m  depth  (Fig  2).  Core-top  P. 
obliquiloculata  temperature  estimates  at  23GGC,  70GGC,  and  69-3  are  also 
consistent  with  an  80-100  m  depth  habitat. 

3.2  Downcore  Results 

The  downcore  P.  obliquiloculata  5180Caidte  records  are  shown  in  Figure  3-7.  At 
69-3  (Fig  3)  the  glacial  5180Caidte  values  were  approximately  0%o.  The  deglacial 
decrease  in  518Ocaidte  began  at  ~19  kyr  BP,  with  a  plateau  from  approximately  14  kyr 
BP  to  11  kyr  BP.  The  518Ocaidte  continued  to  decrease  until  ~8  kyr  BP,  reaching  - 
1.6%o,  and  then  remained  approximately  constant  until  the  late  Holocene.  The 
70GGC  core  (Fig  4)  did  not  reach  the  LGM  and  so  the  record  begins  during  the 
deglacial.  There  is  no  change  in  the  518Ocaidte  from  the  start  of  the  record  at  14  kyr 
BP  until  ~11.5  kyr  BP;  during  this  period  5180Caidte  was  — 0.5%o.  After  11.5  kyr  BP, 
the  5180Caidte  continued  to  decrease  until  ~7  kyr  BP,  reaching  a  value  of  -2%o,  where 
it  remained  for  the  rest  of  the  Holocene.  The  23GGC  record  (Fig  5)  spans  10  kyr  BP 
to  the  late  Holocene  and  618Ocaidte  values  were  very  similar  to  70GGC.  The  518Ocaidte 
record  of  MD41  (Fig  6)  shows  considerable  variability,  likely  due  to  the  low 
sampling  resolution,  from  22  kyr  BP  until  12  kyr  BP,  with  an  average  value  of  - 
0.5%o.  5180Caidte  values  reached  a  minimum  of  approximately  -2%o  at  7  kyr  BP  and 
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then  increased  by  over  l%o  by  4  kyr  BP  (the  top  of  the  core).  At  136GGC  (Fig  7) 
§180Caicite  values  were  approximately  0%o  from  22  kyr  BP  until  16  kyrBP.  6180Caidte 
then  decreased  to  -2%o  by  10  kyr  BP,  and  remained  constant  until  modern. 

Our  downcore  records  of  Mg/Ca-based  temperature  estimates  are  shown  in 
Figures  3-7.  Generally,  thermocline  temperatures  were  nearly  as  warm  during  the 
LGM  as  in  the  late  Holocene.  Thermocline  temperatures  peaked  at  ~10.5  ky  BP  in 
the  Makassar  Strait  (23GGC  and  70GGC)  (Fig  4  and  5)  and  Flores  Sea  (136GGC)  (Fig 
7),  remained  warm  until  5.5  ky  BP,  and  then  gradually  cooled  by  ~1°C.  In  the  Savu 
Basin  (69-3)  (Fig  3),  where  glacial  and  interglacial  thermocline  temperatures  were 
similar,  temperatures  peaked  at  10.5  ky  BP,  but  the  cooling  trend  began  earlier  and 
was  approximately  1°C  larger  than  in  the  Makassar  Strait.  The  temperature  trends 
from  the  Sulu  Sea  (MD41)  (Fig  6)  are  similar  to  those  observed  in  the  Savu  Sea. 
Thermocline  temperatures  peak  at  10.5  kyr  BP  and  cool  by  ~3°  C  until  4.3  kyr  BP, 
when  the  record  ends. 

4.  Discussion 

4.1  Geochemistry  ofP.  obliquiloculata  as  a  thermocline  indicator 

Numerous  lines  of  evidence,  including  data  from  plankton  tows  [Fairbanks  et 
al,  1982],  sediment  traps  [Anand  etal.,  2003;  Mohtadi  etal,  2009],  and  coretops 
[Anand  et  al.,  2003;  Sadekov  etal.,  2009;  Xu  etal,  2006;  Xu  etal.,  2008;  Zuraida  etal, 
2009],  suggest  that  within  the  modern  tropics,  P.  obliquiloculata  has  an  upper 
thermocline  depth  habitat.  Our  data  are  consistent  with  these  studies. 
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An  issue  that  requires  exploration  is  the  use  of  the  difference  between  the 
Mg/Ca-based  temperature  of  G.  ruber  and  P.  obliquiloculata  in  order  to  derive  a 
thermal  gradient.  Previous  studies  have  interpreted  this  difference  in  opposite 
ways.  Xu  et  al.  [2008]  suggest  that  a  small  thermal  gradient  is  indicative  of  a  deep 
thermocline.  Using  data  from  the  same  core,  Sarnthein  et  al.  [2011]  claim  that  a 
small  thermal  gradient  indicates  a  shallow  thermocline  with  more  intense 
upwelling.  Steinke  et  al.  [2010]  interpret  a  large  difference  to  indicate  a  smaller 
mixed  layer  depth,  which  is  qualitatively  equivalent  to  the  Xu  et  al.  [2008] 
interpretation.  While  it  is  difficult  to  determine  thermocline  structure  with  only  two 
points,  we  suggest  that  the  coretop  data  can  guide  us  to  the  appropriate 
interpretation.  In  upwelling  regions,  such  as  along  the  southern  coast  of  Java,  our 
coretop  data  show  a  larger  difference  between  the  Mg/Ca  derived  temperature  of  G. 
ruber  and  P.  obliquiloculata  (Fig  2).  To  highlight  this,  we  combine  our  data  with  the 
extensive  coretop  work  of  Mohtadi  etal.  [2011a]  and  bin  all  data  into  regional 
averages  (Fig  8).  Thus,  coretop  evidence  strongly  suggests  that  the  most  valid 
interpretation  of  the  reconstructed  thermal  gradient  is  consistent  with  that  of  Xu  et 
al.  [2008]  and  Steinke  etal.  [2010],  and  is  not  consistent  with  Sarnthein  etal. 

[2011].  A  larger  difference  between  the  temperature  of  G.  ruber  and  P. 
obliquiloculata  is  found  in  areas  of  upwelling,  while  a  smaller  thermal  gradient  is 
found  in  areas  with  a  more  stratified  water  column. 

Previous  work  has  examined  the  relationship  of  the  difference  in  518Ocaicite 
between  G.  ruber  and  P.  obliquiloculata  to  infer  upper  water  column  structure. 
Consistent  with  the  Mg/Ca-based  temperature  difference  described  above,  coretop 
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results  from  the  South  China  Sea  suggest  that  a  larger  difference  between  G.  ruber 
and  P.  obliquiloculata  occurs  in  areas  with  more  upwelling  [Tian  etal,  2005].  This 
result  has  also  been  generally  found  when  comparing  the  5180Caidte  of  other  surface 
and  thermocline  dwelling  foraminifera  [Mohtadi  etal.,  2011a;  Ravelo  and  Fairbanks, 
1992;  Ravelo  and  Andreasen,  1999]. 

Although  the  coretop  Mg/Ca  and  5180Caidte  differences  between  G.  ruber  and 
P.  obliquiloculata  are  consistent  in  the  water  column  structure  they  suggest  (Fig  8), 
if  we  attribute  the  518Ocaidte  changes  entirely  to  temperature,  the  results  are 
quantitatively  different.  We  use  pooled  Eastern  Indian  Ocean  data  from  this  paper, 
Mohtadi  et  al.  [2011a]  and  Gibbons  et  al.  [in  revision],  to  examine  the  predicted 
water  column  structure  difference  between  an  area  of  intense  upwelling  (the  ]ava 
Basin  or  ]B  in  Fig  8)  and  an  area  with  a  stratified  water  column  (the  Nias  Basin  or 
NB).  The  G.  ruber  -  P.  obliquiloculata  Mg/Ca-based  temperature  in  the  Java  Basin  is 
8.4°C  and  at  the  Nias  basin  it  is  6.2°C.  The  G.  ruber  -  P.  obliquiloculata  difference  in 
5180Caiidte  is  -1.4%o  at  the  ]ava  Basin  and  -0.7%o  at  the  Nias  basin.  If  we  convert  the 
5180Caidte  difference  to  temperature  (using  the  low-light  0.  universa  equation  from 
Bemis  et  al.  [1998])  this  is  equivalent  to  a  6.7°C  temperature  difference  in  the  Java 
Basin  and  a  3.4°C  temperature  difference  in  the  Nias  Basin.  In  order  to  reconcile 
these  differences,  the  subsurface  water  must  have  more  depleted  518Oseawater  values 
than  the  surface  waters.  Modern  observations  indicate  that  the  subsurface  water  in 
the  Eastern  Indian  Ocean  is  saltier  than  the  surface  water  (eg.  WOA05  [Antonov  et 
al.,  2006]).  Thus  if  subsurface  waters  actually  have  more  depleted  518Oseawater,  lower 
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5180Seawater  values  are  associated  with  saltier  conditions.  Generally,  however,  salinity 
and  5180Seawater  are  positively  correlated  [LeGrande  and  Schmidt,  2006]. 

Due  to  the  scarcity  of  5180  seawater  measurements  in  the  region,  particularly  at 
the  subsurface,  we  cannot  directly  examine  whether  5180  seawater  values  are  more 
depleted  with  depth  in  the  region.  We  examine  this  problem  indirectly  by  utilizing 
the  conservative  nature  of  518Oseawater  and  taking  advantage  of  the  relative 
abundance  of  surface  518Oseawater  measurements  compared  to  subsurface  518Oseawater 
measurements.  At  several  core  sites  in  the  IndoPacific  Warm  Pool,  we  use  the  Total 
Matrix  Intercomparison  method  developed  by  Gebbie  and  Huybers  [2010].  Total 
Matrix  Intercomparison  uses  six  water  properties  to  trace  back  the  surface  source 
region(s)  of  a  water  parcel.  This  method  allows  for  multiple  source  regions  to 
contribute  to  each  water  parcel.  We  identify  the  major  source  regions  at  each 
location  and  use  an  average  5180  seawater  surface  value  from  the  gridded  data  of 
LeGrande  and  Schmidt  [2006]  to  predicted  the  518Oseawater  at  varying  depths  in  the 
IndoPacific  Warm  Pool.  This  approach  suggests  that  518Oseawater  does  become  more 
depleted  with  depth  in  this  region,  with  an  approximate  change  of  -0.06%o  per 
100m.  While  this  may  contribute  to  the  difference  between  the  5180Caidte  and  Mg/Ca- 
based  temperature  gradients,  it  is  not  large  enough  to  resolve  the  discrepancy 
between  the  thermal  gradient  based  on  Mg/Ca  versus  5180Caidte. 

Previous  authors  [Steinke  etal,  2010]  have  noted  the  apparent  difficulty  in 
reconciling  the  518Ocaidte  difference  of  G.  ruber  -  P.  obliquiloculata  and  the  Mg/Ca- 
based  temperature  difference.  While  Steinke  etal.  [2010]  attributed  this  to  the  local 
hydrology  of  the  South  China  Sea,  our  data  suggest  this  discrepancy  may  be 
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widespread.  Others  [Cleroux  etal,  2008;  Mohtadi  etal,  2011a]  have  used  the 
predicted  temperature  based  on  modern  oceanographic  data  and  an  assumed 
5180seawater/salinity  relationship  to  recalculate  the  Mg/Ca-temperature  relationship. 
This  approach  forces  a  reconciliation  of  the  Mg/Ca-based  temperature  and  5180Caidte 
and  suggests  that  the  sensitivity  of  the  Mg/Ca  of  P.  obliquiloculata  to  temperature 
change  is  approximately  half  that  found  by  other  studies  in  areas  where  the 
oceanographic  conditions,  including  the  5180Seawater/salinity  relationship  were  better 
known  [Anand  etal.,  2003].  Further  work,  specifically  on  the  modern  518Oseawater 
relationship  to  water  mass  properties,  the  depth  habitat  of  P.  obliquiloculata  via 
plankton  tows,  and  the  sensitivity  of  P.  obliquiloculata  Mg/Ca  and  5180Caidte  to 
temperature  changes  via  culturing  work  are  required  to  clarify  the  situation. 

4.2  Changes  in  the  upper  water  column  -  The  LGM  and  deglacial 

We  interpret  our  downcore  data  assuming  that  our  temperature  estimates 
using  Mg/Ca  and  the  modified  Anand  etal.  [2003]  multispecies  calibration  are 
correct.  This  allows  us  to  compare  our  results  with  previous  work  that  used  the 
Mg/Ca  of  G.  ruber  and  P.  obliquiloculata  to  infer  changes  in  thermocline  structure. 

We  compare  these  interpretations  to  the  thermocline  structure  derived  from  the 
5180Caidte-  Finally,  we  examine  possible  reasons  for  the  discrepancy  between  these 
two  methods. 

4.2.1  The  LGM  and  deglacial  -  climatic  implications  based  on  Mg/Ca 
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70GGC  and  23GGC  do  not  contain  LGM  sediment,  but  in  all  our  cores  that 


have  data  from  the  LGM,  reconstructed  thermocline  temperatures  are  no  more  than 
1°C  cooler  than  modern.  This  is  a  striking  regional  feature  that  is  observed  in  the 
Sulu  Sea  (MD41,  Fig  6),  the  Flores  Sea  (136GG,  Fig  7),  and  the  Savu  Basin  (69-3,  Fig 
3).  It  has  also  been  documented  in  the  Timor  Sea  [Xu  et  al.,  2008],  the  Timor  Strait 
[Holbourn  etal,  2011],  and  the  South  China  Sea  [Steinke  etal,  2010;  Yang  and 
ZhiMin,  2009].  The  only  Mg/Ca-based  P.  obliquiloculata  record  that  shows  a  larger 
cooling  (nearly  2°C)  is  MD67  from  the  Mindanao  dome  [Bolliet  et  al,  2011]. 

The  warm  LGM  thermocline  temperatures  create  a  smaller  vertical  thermal 
gradient  between  the  surface  and  the  thermocline  (Fig  9),  which  our  coretop  data 
suggest  indicate  more  stratified  oceanic  conditions  (Fig  8).  Greater  LGM 
stratification  is  consistent  with  previous  work  from  the  Western  Pacific  Warm  Pool. 
Planktic  foraminifera  abundances  suggest  that  the  thermocline  was  deeper 
throughout  the  Western  Pacific  Warm  Pool  at  the  LGM  [Andreasen  and  Ravelo, 

1997].  Using  a  combination  of  Mg/Ca-based  SST  reconstruction  (using  G.  ruber], 
alkenone  (Uk’37),  and  foraminferal  transfer  functions,  de  Garidel-Thoron  etal. 

[2007]  suggest  LGM  SSTs  were  ~2.5°C  cooler  and  subsurface  temperatures  were 
1.8°C  cooler,  again  implying  more  stratified  conditions. 

Andreasen  and  Ravelo  [1997]  and  de  Garidel-Thoron  [2007]  suggest  their 
data  are  consistent  with  an  enhanced  Walker  circulation  at  the  LGM.  Support  for  this 
mechanism  comes  from  a  compilation  of  six  LGM  modeling  simulations,  four  of 
which  show  an  enhanced  Walker  Circulation.  In  these  four  glacial  simulations,  the 
upper  water  column  is  more  stratified,  and  thermocline  temperatures  show  less 
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cooling  than  surface  temperatures  [DiNezio  etal,  2011].  Less  cooling  in  the  glacial 
thermocline  compared  to  the  surface  is  also  simulated  in  the  eastern  equatorial 
Indian  Ocean  [Otto-Bliesner  et  al,  2009].  The  large  regional  extent  of  the  relatively 
warm  glacial  thermocline  in  these  models  is  consistent  with  our  proxy  data. 
However,  we  note  our  proxy  data  show  yet  warmer  thermocline  temperatures  than 
suggested  by  the  average  of  the  modeling  simulations. 

4.2.2  The  LGM  and  deglacial  -  climatic  implications  based  on  S18Ocaicite 
Though  the  G.  ruber  -  P.  obliquiloculata  Mg/Ca-based  temperature 
differences  suggest  a  deeper  thermocline  (smaller  thermal  gradient)  at  the  LGM, 
these  results  are  inconsistent  with  the  5180Caidte  results.  In  most  continuous  records 
(MD41,  69-3,  MD90  [Steinke  et  al,  2010],  and  MD67  [Bolliet  et  al,  2011],  the 
5180Caicite  difference  between  G.  ruber  and  P.  obliquiloculata  is  approximately  the 
same  at  the  LGM  as  it  is  at  the  modern  (within  0.3%o  or  approximately  1°C  if 
attributed  to  temperature),  suggesting  either  no  change  or  a  slightly  less  stratified 
IndoPacific  Warm  Pool  at  the  LGM.  This  is  in  contrast  to  the  Mg/Ca-based 
temperature  reconstructions  which  suggest  a  more  stratified  LGM,  with  the 
temperature  gradient  approximately  2°C  smaller  than  modern.  MD94  [Yang  and 
ZhiMin,  2009]  in  the  South  China  Sea  is  the  only  record  in  which  the  5180Caidte 
gradient  and  the  Mg/Ca-based  temperature  gradient  show  the  same  trend;  the 
5180Caidte  records  suggest  a  modest  deepening  of  the  thermocline  at  the  LGM,  though 
still  less  than  based  on  the  Mg/Ca-based  temperature. 
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As  discussed  above,  this  apparent  inconsistency  between  the  Mg/Ca-based 
temperature  and  5180Caidte-based  temperature  gradients  may  result  from  more 
depleted  5180  seawater  in  the  thermocline  than  surface.  In  this  case,  a  deeper 
thermocline  would  bathe  P.  obliquiloculata  in  waters  that  were  both  warmer  and 
had  a  higher  5180  seawater  value.  Even  with  a  deepening  of  the  thermocline,  there 
could  be  no  change  in  the  518Ocaicite  of  P.  obliquiloculata  as  the  warming  and  higher 
5180seawater  values  could  offset  each  other. 

4.2.3  The  LGM  and  deglacial  -  potential  biases:  is  there  a  thermal  limit  to  P. 
obliquiloculata  or  is  salinity  influencing  our  Mg/Ca-based  temperature 
estimates? 

We  examine  two  potential  sources  of  bias  in  our  records:  firstly,  whether  the 
lower  thermal  limit  for  P.  obliquiloculata  is  approached  during  the  LGM,  shifting  P. 
obliquiloculata  towards  warmer  conditions  and  secondly,  whether  the  generally 
saltier  conditions  of  the  LGM  could  raise  Mg/Ca  ratios  independent  of  temperature 
[Arbuszewski  etal,  2010;  Mathien-Blard  and  Bassinot,  2009],  We  explore  each  of 
these  possibilities  below. 

In  the  modern  oceans  P.  obliquiloculata  is  an  abundant  species  in  the 
IndoPacific,  but  abundances  quickly  decline  in  cooler  waters  [Be  and  Hutson,  1977; 
Cullen  and  Prell,  1984;  Martinez  etal.,  1998;  Ujiie  and  Ujiie,  2000].  The  relatively 
warm  Mg/Ca-based  thermocline  temperatures  at  most  locations  during  the  LGM  as 
described  above  could  be  due  to  P.  obliquiloculata  migrating  upward  in  the  water 
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column,  switching  its  growth  to  the  warm  season  or  only  growing  during  warm 
years. 

Support  for  a  thermal  limit  come  from  regional  abundance  records.  Within 
the  South  China  Sea,  P.  obliquiloculata  disappears  almost  entirely  during  Heinrich 
Event  1  [Yang  and  ZhiMin,  2009].  In  the  Sulu  Sea  P.  obliquiloculata  disappears 
during  the  last  glacial  period,  even  while  G.  ruber  abundance  remains  high  [Linsley 
and  Thunell,  1990],  suggesting  that  P.  obliquiloculata  has  a  narrower  optimal  habitat 
than  G.  ruber.  Abundance  counts  from  within  the  Makassar  Strait  [Ding  etal,  2006] 
show  similar  results  to  Linsley  and  Thunell  [1990]  suggesting  that  even  within  the 
heart  of  the  IndoPacific  Warm  Pool,  the  lower  thermal  limit  for  P.  obliquiloculata 
was  approached  during  glacial  conditions.  If  we  attribute  changes  in  the  5180Caidte 
records  primarily  to  temperature,  these  data  suggest  P.  obliquiloculata  is  not  biased 
by  cold  temperatures,  but  as  discussed  above,  the  uncertainty  in  interpreting  past 
(or  present)  5180  seawater  values  make  it  difficult  to  rule  out  a  thermal  limit. 

Another  possibility  is  that  our  warm  glacial  thermocline  reconstructions  are 
due  to  the  influence  of  salinity  on  Mg/Ca.  Evidence  is  increasing  that  at  high 
salinities,  Mg/Ca  values  are  elevated  in  G.  ruber  [Arbuszewski  etal,  2010;  Mathien- 
Blard  and  Bassinot,  2009].  Even  though  there  is  some  debate  as  to  the  exact 
mechanisms  of  salinity  influence  on  Mg/Ca,  it  appears  that  at  salinities  below  35 
[Arbuszewski  et  al,  2010]  or  36  [Gibbons  et  al,  in  revision],  the  effect  is  negligible. 
CTD  data  (not  shown)  and  World  Ocean  Atlas  05  data  [Antonov  et  al.,  2006]  suggest 
salinities  at  100  m  at  our  core  sites  are  less  than  35.  The  mean  ocean  was  ~1  saltier 
at  the  LGM  [Adkins  et  al.,  2002],  leaving  thermocline  salinities  at  our  sites  still  below 
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the  likely  threshold  for  a  salinity  influence  on  Mg/Ca  if  this  was  the  only  factor 
affecting  salinity.  However,  within  the  ITF,  the  thermocline  was  likely  even  saltier 
due  to  the  absence  of  freshwater  from  the  Java  Sea  [Linsley  et  al.,  2010],  which 
freshens  the  thermocline  today  [Gordon,  2005].  Thus,  it  is  possible  that  thermocline 
salinities  exceeded  36,  and  contributed  to  the  high  Mg/Ca  during  the  LGM  in  our  ITF 
sites  (136GGC,  and  69-3).  Moreover,  the  effect  of  salinity  on  the  Mg/Ca  ofP. 
obliquiloculata  is  not  established. 

We  can  test  whether  the  salinity  effect  is  a  reasonable  explanation  for  our 
warm  LGM  temperatures.  We  first  assume  that  there  was  no  change  in  the  thermal 
gradient  at  our  sites.  In  other  words,  G.  ruber  and  P.  obliquiloculata  should  have  both 
been  cooler  by  the  same  amount  at  the  LGM.  We  use  the  difference  in  our  Mg/ Ca- 
based  G.  ruber  temperatures  between  20-18  kyr  BP  and  2-0  kyr  BP  to  assign  the 
magnitude  of  thermocline  cooling  at  the  LGM.  We  subtract  this  cooling  from  our 
"modern”  P.  obliquiloculata  temperature,  to  estimate  P.  obliquiloculata' s  LGM 
temperature.  We  then  calculate  predicted  Mg/Ca  values  based  on  these  estimated 
LGM  values  and  the  modified  Anand  [2003]  multispecies  calibration.  We  subtract 
this  predicted  Mg/Ca  from  our  measured  Mg/Ca,  and  use  this  value  in  a  similar 
manner  to  the  "excess  Mg/Ca"  explored  by  Arbuszewski  etal.  [2010],  We  can  use 
the  dataset  of  excess  Mg/Ca  and  salinity  compiled  by  Gibbons  etal.  [in  revision]  to 
predict  the  salinity  associated  with  an  excess  Mg/Ca  value  (Fig  10).  At  salinities 
above  36,  the  relationship  between  excess-Mg/Ca  and  salinity  can  be  described  by: 

Excess  Mg/Ca  =  1.4018*salinity  -  50.449 


117 


If  we  assume,  instead  the  critical  threshold  for  salinity  is  35,  the  relationship 
changes  only  slightly: 

Excess  Mg/Ca  =  1.1032*salinity  -  39.456 

Using  our  excess  Mg/Ca  values,  we  solve  for  salinity.  The  salinities  required 
to  produce  the  LGM  excess  Mg/Ca  are  36.2  at  136GGC  and  69-3  and  36.4  at  MD78 
using  the  first  equation  and  36.0  at  136GGC  and  69-3  and  36.2  at  MD78  using  the 
second  equation.  The  relatively  modest  increase  in  salinity  required  to  affect  Mg/Ca- 
based  temperature  estimates  suggests  this  is  a  potential  concern  for  the  use  of  P. 
obliquiloculata  at  the  LGM,  especially  as  the  sensitivity  of  Mg/ Ca  in  P.  obliquiloculata 
to  salinity  is  unknown. 

4.3  Changes  in  the  upper  water  column  during  the  Holocene  -  The  flooding  of 
the  Sunda  Shelf 

We  suggest  that  during  the  Holocene,  the  Mg/Ca-based  temperature  of  P. 
obliquiloculata  is  a  reliable  proxy  for  thermocline  temperature.  P.  obliquiloculata  is 
relatively  abundant  [Ding  et  ai.,  2006;  Linsiey  and  Thunell,  1990],  suggesting  it  is 
well  within  its  prime  habitat,  and  after  the  disintegration  of  the  ice  sheets  [Clark  and 
Mix,  2002]  salinities  in  the  IndoPacific  Warm  Pool  are  likely  low  enough  where 
Mg/Ca  is  primarily  dependent  on  temperature. 

Previous  work  has  suggested  that  the  flooding  of  the  Sunda  Shelf  that 
occurred  at  10  kyr  BP  cooled  the  thermocline  of  the  ITF  outflow  passages  [Hoibourn 
etal.,  2011;  Xu  et  al,  2008;  Xu  et  al,  2010].  This  hypothesis  is  based  on  the  idea  that 
there  is  a  reduction  of  warm  surface  waters  that  pass  through  the  Makassar  Strait 
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when  the  Sunda  Shelf  floods.  By  contrast,  others  have  argued  based  on  modern  data 
[Fang  etal,  2010;  Qu  etal,  2006]  and  paleo  reconstructions  [Gibbons  et  al,  in 
revision]  that  the  decreased  heat  flow  through  the  Makassar  Strait  would  be 
balanced  by  heat  from  the  ]ava  Sea.  As  a  result,  there  would  be  no  net  change  in  heat 
transport  or  thermocline  stratification  as  a  result  of  the  flooding  of  the  Sunda  Shelf. 

We  test  whether  the  flooding  of  the  Sunda  Shelf  changed  the  thermocline 
thermal  gradient  of  either  the  Makassar  Strait  or  the  outflow  passages  by  computing 
the  Kendall  x  coefficient  [Kendall  and  Gibbons,  1990]  for  each  core.  The  Kendall  x 
coefficient  is  a  nonparamentric  measure  of  trend,  which  can  range  from  -1  to  1.  It  is 
robust  against  both  outliers  and  inaccuracies  in  chronology.  In  neither  69-3  (x  =  - 
0.18;  p  =  0.13),  nor  MD78  (x  =  -0.20;  p  =  0.12)  is  there  a  significant  increase  in 
thermal  gradient  from  10-6  kyr  BP,  which  suggests  that  the  flooding  of  the  Sunda 
Shelf  did  not  change  the  thermal  structure  of  the  outflow  passages. 

We  note  that  the  finding  of  no  change  in  the  thermal  gradient  due  to  the 
flooding  of  the  Sunda  Shelf  is  sensitive  to  the  Mg/ Ca  temperature  equation  used.  If 
we  instead  use  the  equation  suggested  by  Cleroux  et  al.  [2008]  and  Mohtadi  et  al. 
[2011a],  there  is  a  dramatic  increase  in  the  thermal  gradient  at  69-3  (x  =  -0.27;  p  = 
0.02)  from  10-6kyrBP.  The  increase  in  thermal  gradient  at  MD78  is  also  larger, 
though  there  is  still  no  significant  trend  at  MD78  (x  =  -0.19;  p  =  0.12).  In  light  of  the 
relatively  modest  Holocene  temperature  variation  in  the  surface  waters  (Fig  3 
shows  a  2°C  cooling),  the  extremely  large  Holocene  thermocline  temperature 
variation  produced  by  this  equation  (a  Holocene  cooling  of  4.5°C  in  69-3),  suggests 
that  this  calibration  may  overestimate  past  temperature  changes.  Using  this 
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equation,  the  thermal  gradient  during  the  early  Holocene  is  only  2.5°C,  suggesting 
extremely  stratified  conditions  at  69-3.  This  is  difficult  to  reconcile  with  strong 
evidence  of  a  northerly  displaced  ITCZ  [Yuan  et  ah,  2004]  and  increased  monsoon 
upwelling  [Mohtadi  etal,  2011b],  which  would  extend  the  upwelling  season  at  69-3. 
Nevertheless,  the  considerable  uncertainty  in  the  Mg/Ca-temperature  calibration  of 
P.  obliquiloculata  requires  further  work. 

Regardless  of  the  calibration,  however,  the  data  support  a  previously  noted 
larger  Holocene  than  glacial  vertical  thermal  gradient  (a  shoaling  thermocline)  in 
the  outflow  passages,  previously  attributed  to  the  flooding  of  the  Sunda  Shelf  [Xu  et 
al,  2008;  Xu  et  al,  2010].  However,  our  records  demonstrate  that  the  change  in 
thermocline  structure  was  not  uniquely  associated  with  the  flooding  of  the  Sunda 
Shelf  at  ~10  kyr  BP  but  instead  occurs  as  a  gradual  trend  that  began  at  the  LGM  and 
continued  until  modern  (Fig  4).  Moreover,  MD41,  from  the  Sulu  Sea,  should  not  have 
been  directly  affected  by  the  flooding  of  the  Sunda  Shelf  [Sathiamurthy  and  Voris, 
2006]  yet  it  has  a  similar  change  in  thermal  gradient  to  the  other  cores.  The 
observation  that  the  increase  in  thermal  gradient  is  both  long  in  duration  and 
widespread  in  occurrence  suggests  that  it  cannot  be  explained  due  to  a  change  in 
thermocline  temperature  resulting  from  the  flooding  of  the  Sunda  Shelf. 

Conclusion: 

More  work  remains  before  we  can  confidently  use  the  geochemistry  of  P. 
obliquiloculata  as  a  robust  indicator  of  thermocline  conditions.  In  both  our  coretop 
and  downcore  records,  it  is  difficult  to  reconcile  the  trends  from  the  Mg/Ca-based 


120 


temperature  with  the  5180Caidte.  If  we  assume,  as  previous  studies  have,  that  the 
Mg/Ca-based  temperature  gradient  between  G.  ruber  and  P.  obliquiloculata  is  an 
indicator  of  thermocline  structure,  our  results  suggest  that  thermocline 
temperatures  were  warm  at  the  LGM  relative  to  surface  water,  and  gradually 
increased  to  a  peak  in  the  early  Holocene.  These  data  are  consistent  with  an 
enhanced  in  the  Walker  Circulation  during  the  LGM  [DiNezio  et  al,  2011].  Our  data 
are  inconsistent  with  the  hypothesis  that  the  flooding  of  the  Sunda  Shelf  caused  a 
cooling  in  the  thermocline  as  inferred  from  the  Mg/Ca-based  temperature  difference 
of  G.  ruber  and  P.  obliquiloculata.  The  increase  in  thermal  gradient  occurs 
throughout  the  deglacial  and  is  observed  in  a  core  that  is  not  located  in  an  outflow 
passage  of  the  ITF  (MD41). 
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Core  ID,  location 

Core  ID  in 

text 

Lat./Long. 

Water  depth 
(m) 

BJ8-03  13GGC,  Flores  Sea 

13GGC 

7°  24'S,  115° 
13'E 

594 

BJ8-03  23GGC,  S  Makassar  Strait 

23GGC 

5°  12'S,  117° 
29'E 

1189 

BJ8-03  70GGC,  SW  Sulawesi,  Makassar 
Strait 

70GGC 

3°  34'S,  119° 

23 'E 

482 

BJ8-03  136GGC,  Flores  Sea 

136GGC 

5°  56'S,  120° 
14'E 

526 

GeoB  10069-3  Savu  Sea 

69-3 

9°  36'S,  120° 
55'E 

1,250 

MD97-2141,  Sulu  Sea 

MD41 

8°  47'N,  121° 
17'E 

3,633 

Table  1  -  Core  locations  and  depths. 
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Core 

Depth 

(cm) 

G. 

ruber/sac 

C-14 

Age 

(yrs) 

G. 

ruber/sac 
calendar 
age  (yrs) 

Depth 

(cm) 

P.  obliq 
C-14 

Age 

(yrs) 

P.  obliq 
calendar 
age 
(yrs) 

Age 

difference 

(yrs) 

69-3 

88 

1970 

1525 

88 

1910 

1453 

72 

69-3 

468 

(478) 

6800 

7328 

(7566) 

478 

6830 

7353 

(213) 

70GGC 

64 

1310 

853 

64.5 

1310 

853 

0 

70GGC 

321 

7190 

7653 

320.5 

7270 

7731 

-78 

70GGC 

362 

10250 

11223 

361.5 

10100 

11132 

91 

Table  2  -  Difference  in  radiocarbon  dates  between  G.  ruber /G.  sacculifer  and  P. 
obliquiloculata.  Dissolution-prone  foraminifera  have  approximately  the  same  age 
as  the  dissolution-resistant  foraminifera  P.  obliquiloculata.  Due  to  a  mismatch  in 
depth  intervals,  the  69-3  P.  obliquiloculata  date  from  478  cm  was  compared  to  the 
predicted  G.  ruber /sacculifer  age  at  478  cm  based  on  the  age  model.  The  predicted 
value  is  listed  in  parentheses. 


123 


Figure  1  -  Map  of  study  location.  Dots  indicate  locations  of  multicores  used  in 
coretop  study.  Diamonds  are  the  location  of  cores  with  new  P.  obliquiloculata 
downcore  records.  The  square  symbol  is  MD78  [Xu  etal,  2008].  Boxes  indicated 
regional  groupings  for  Figure  3  (from  north  to  south  SB:  Simeulue  Basin;  NB:  Nias 
Basin;  NMB:  Northern  Mentawai  Bain;  SMB:  Southern  Mentawai  Basin;  JB:  Java 
Basin;  LB:  Lombok  Basin)  and  core  69-3  is  located  in  the  Savu  Sea.  Dashed  arrow 
indicates  path  of  freshwater  from  the  South  China  Sea  into  the  Indonesian 
Throughflow. 
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Figure  2  -  Core  top  results.  The  Mg/Ca  based  temperature  estimates  of  G.  ruber 
[Gibbons  etal.,  in  revision]  (blue  triangles)  and  P.  obliquiloculata  (purple  squares). 
Coretop  data  are  plotted  against  CTD  data  from  the  coring  cruises.  Dashed  line  is 
satellite  derived  mean  annual  SST  (http://reason.gsfc.nasa.gov/Giovanni/). 
Temperatures  (both  CTD  and  Mg/Ca  based)  are  plotted  on  the  Y-axis.  Isodepth  lines 
are  plotted  through  temperature  data.  The  Mg/Ca-based  temperature  of  P. 
obliquiloculata  is  consistent  with  temperatures  found  between  80-100m. 
Additionally,  in  areas  of  more  upwelling,  such  as  more  eastern  Indian  Ocean  sites, 
the  reconstructed  temperature  difference  between  G.  ruber  and  P.  obliquiloculata  is 
larger. 
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Figure  3  -  Downcore  records  from  69-3  (Savu  Sea).  G.  ruber  (red)  [Gibbons  etal, 
in  revision]  and  P.  obliquiloculata  (blue)  5180Caicite,  Mg/Ca-based  temperature 
estimates,  and  518Oseawater  estimates.  Note  the  inverted  axes  for  5180.  Gradients  are 
G.  ruber  -  P.  obliquiloculata  with  the  axes  orientation  as  in  Figure  3  (gradients 
corresponding  to  upwelling  conditions  in  the  coretop  data  are  up).  Axes  are  scaled 
so  that  l%o  5180  corresponds  to  4°C. 
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Thermocline  temp  (°C) 


Figure  4  -  Downcore  records  from  70GGC  (Southern  Makassar  Strait).  G.  ruber 
(red)  [Linsley  et  al,  2010]  andP.  obliquiloculata  (blue)  5180Caidte,  Mg/Ca-based 
temperature  estimates,  and  618Oseawater  estimates.  Note  the  inverted  axes  for  S180. 
Gradients  are  G.  ruber  -  P.  obliquiloculata  with  the  axes  orientation  as  in  Figure  3 
(gradients  corresponding  to  upwelling  conditions  in  the  coretop  data  are  up).  Axes 
are  scaled  so  that  l%o  5180  corresponds  to  4°C. 


127 


Thousands  of  years  BP 


Figure  5  -  Downcore  records  from  23GGC  (Southern  Makassar  Strait).  G.  ruber 
(red)  and  P.  obliquiloculata  (blue)  S18Ocaicite,  Mg/Ca-based  temperature  estimates, 
and  5180Seawater  estimates.  Note  the  inverted  axes  for  S180.  Gradients  are  G.  ruber  -  P. 
obliquiloculata  with  the  axes  orientation  as  in  Figure  3  (gradients  corresponding  to 
upwelling  conditions  in  the  coretop  data  are  up).  Axes  are  scaled  so  that  l%o  5180 
corresponds  to  4°C. 
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Figure  6  -  Downcore  records  from  MD41  (Sulu  Sea).  G.  ruber  (red)  and  P. 
obliquiloculata  (blue)  S18Ocaidte,  Mg/Ca-based  temperature  estimates,  and  518Oseawater 
estimates.  Note  the  inverted  axes  for  S180.  Gradients  are  G.  ruber  -  P.  obliquiloculata 
with  the  axes  orientation  as  in  Figure  3  (gradients  corresponding  to  upwelling 
conditions  in  the  coretop  data  are  up).  Axes  are  scaled  so  that  l%o  S180  corresponds 
to  4°C. 
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Figure  7  -  Downcore  records  from  136GGC  (Flores  Sea).  G.  ruber  (red)  and  P. 
obliquiloculata  (blue)  5180Caidte,  Mg/Ca-based  temperature  estimates,  and  518Oseawater 
estimates.  Note  the  inverted  axes  for  S180.  Gradients  are  G.  ruber-  P.  obliquiloculata 
with  the  axes  orientation  as  in  Figure  3  (gradients  corresponding  to  upwelling 
conditions  in  the  coretop  data  are  up).  Axes  are  scaled  so  that  l%o  5180  corresponds 
to  4°C. 
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Thermocline  temp  (°C) 


Figure  8  -  Coretop  summary  -  G.  ruber  (red)  and  P.  obliquiloculata  (blue)  Mg/Ca- 
based  temperature  estimates,  5180Caidte,  and  518Oseawater  estimates  [this  paper, 
Gibbons  etal,  in  revision;  Mohtadi  et  al,  2011a].  The  gradients  between  G.  ruber  and 
P.  obliquiloculata  are  shown  in  the  bottom  panel.  Data  are  pooled  into  the  regions 
indicated  in  Figure  1.  Areas  that  have  increased  upwelling  in  the  modern  ocean, 
such  as  the  Java  Basin  (JB),  have  a  larger  Mg/Ca-based  temperature  gradient  and  a 
more  negative  5180Caidte  gradient.  Areas  that  are  more  stratified  in  the  modern 
ocean,  such  as  the  Nias  Basin  (NB),  show  the  opposite,  a  smaller  Mg/Ca  based 
temperature  gradient  and  a  less  negative  518Ocaidte  gradient.  Axes  are  scaled  so  that 
l%o  5180  corresponds  to  4°C. 
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Figure  9  -  Temperature  and  S18Ocaicite  gradients  from  selected  timeslices.  The 

gradients  are  G.  ruber-  P.  obliquiloculata.  Timeslices  are  2,000  year  averages  and 
are  centered  on  19kyr  BP,  llkyr  BP,  5kyr  BP,  and  1  kyr  BP.  Axes  are  scaled  so  that 
l%o  S180  corresponds  to  4°C.  Error  bars  are  +/-  a.  For  Mg/Ca-based  temperature, 
errors  were  calculated  as  described  in  Chapter  2,  and  a2  Temp  grad  —  CX2  ruber  temp  +  O  2 
obiiq  temp,  cx  of  the  5180Caidte  is  based  on  ~50  G.  ruber  replicates  in  23GGC  and  is 
0.22%o.  It  is  assumed  that  P.  obliquiloculata  has  the  same  error.  For  each  timeslice  cx 
^8180calcite  ave  —  (cx  2si80caidte)/N,  where  N  is  the  number  of  samples,  and  cx2  sisocaicite  grad  = 
CX  2  8180calcite  ave  ruber  +  CX  2  S180calcite  ave  obiiq- 


132 


7  n 


-3  — 

30  32  34  36  38  40 

Salinity 

Figure  10  -  Excess  Mg/Ca  and  salinity  from  G.  ruber .  Lines  were  used  to 
calculate  the  excess  salinity  required  to  cause  our  warm  P.  obliquiloculata  Mg/Ca  - 
based  temperature  reconstruction  at  the  LGM.  Solid  line  shows  the  relationship 
between  excess  Mg/Ca  and  salinity,  assuming  a  salinity  threshold  of  36  (y  =  1.40 18x 
-  50.449;  R2  =  0.62)  and  the  dashed  line  is  the  relationship  assuming  a  salinity 
threshold  of  35  (y  =  1.1032X  -  39.456;  R2  =  0.54). 
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Conclusions 


This  work  has  provided  insight  into  the  dynamics  of  the  hydrology  and 
oceanography  of  the  Indonesian  Throughflow.  The  data  presented  here  strongly 
suggest  that  the  Java  Sea  inflow  does  not  have  a  cooling  effect  on  the  outflow  of  the 
Indonesian  Throughflow.  The  appearance  of  a  cool  tongue  of  Indonesian  water  in 
the  equatorial  Indian  Ocean  is  more  likely  due  to  the  vigorous  upwelling  in  the 
Banda  Sea.  This  work  also  shows  that  the  Intertropical  Convergence  Zone  migrates 
in  concert  with  high  latitude  Northern  Hemisphere  cooling  events,  possibly  in 
response  to  the  decreased  interhemispheric  temperature  gradient.  Additionally,  the 
continued  need  for  ongoing  proxy  development  work  is  clearly  demonstrated.  I 
show  the  tremendous  difficulty  in  reconciling  the  518Ocaiicte  and  Mg/Ca  of  P. 
obliquiloculata  data  variations,  despite  these  proxies  being  in  widespread  use. 

There  are  several  ways  to  move  this  work  forward.  Of  upmost  importance  is 
the  need  for  more  518Oseawater  measurements,  particularly  in  the  subsurface.  Not  only 
would  this  aid  our  understanding  of  modern  hydrology  in  the  region,  it  is  vital  for 
our  interpretation  of  thermocline  (and  benthic)  records. 

Additionally,  more  work  on  the  geochemistry  of  P.  obliquiloculata  is  required. 
Perhaps  there  is  a  species-specific  Mg/Ca  or  5180  temperature  relationship  that  is 
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more  suitable.  Given  the  difficulty  in  simultaneously  assessing  depth  habitat  and 
geochemical  sensitivity  to  temperature,  culturing  work  would  yield  the  most 
straightforward  answer. 

There  is  more  work  to  be  done  on  constraining  the  migration  of  the 
Intertropical  Convergence  Zone.  Within  the  IndoPacific  Warm  Pool,  the  data  show  a 
clear  band  of  no  temporal  change  in  5180  seawater  over  the  past  20,000  years.  I 
speculate  that  during  high  latitude  Northern  Hemisphere  coolings  these  locations 
experience  an  increase  in  boreal  summer  precipitation  and  a  decrease  in  boreal 
winter  precipitation.  The  geographic  coverage  of  the  records  I  use  allows  me  to 
constrain  the  migration  of  the  Intertropical  Convergence  Zone.  In  the  Eastern 
Equatorial  Pacific,  the  southern  most  518Oseawater  record  in  this  study  is  from  a  core 
retrieved  from  1°S.  If  we  could  reconstruct  518Oseawater  at  cores  farther  south,  I  would 
expect  to  find  a  location  where  5180  seawater  does  not  change  during  times  of  high 
latitude  Northern  Hemisphere  cold  events.  Even  farther  south,  I  would  expect  to  see 
depleted  518Oseawater  values.  These  findings  would  confirm  an  ITCZ  mechanism,  as 
we  currently  cannot  rule  out  that  changes  in  the  Eastern  Equatorial  Pacific  are  due 
to  variations  in  moisture  transport  across  the  Isthmus  of  Panama. 

We  can  also  test  whether  changes  in  the  interhemispheric  temperature 
gradient  consistently  respond  to  hydrologic  changes.  Our  longest  record  (69-3) 
shows  a  drying  event  that  is  possibly  consistent  with  Heinrich  Event  3,  but  no 
expression  of  Heinrich  Event  2.  If  this  finding  holds  with  increased  sampling 
resolution  and  better  dating,  then  it  would  suggest  that  changes  in  the 
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interhemispheric  temperature  gradient  did  not  always  occur  during  Heinrich 
Events.  This  is  a  hypothesis  that  can  likely  be  explored  using  existing  data. 

Finally,  the  mechanisms  1  propose  for  temperature  controls  in  the  Indonesian 
Throughflow  can  be  more  thoroughly  examined.  My  hypothesis  requires 
intermediate-depth  water  temperature  changes  to  resemble  Antarctic  air 
temperature  changes.  By  deriving  temperatures  from  benthic  foraminifera  from 
thermocline  and  intermediate  water  depth  cores  in  the  region,  this  oceanic 
transport  mechanism  can  be  tested. 

Most  importantly,  as  we  move  forward  and  generate  new  observational  and 
paleo  data  from  the  Indonesian  Throughflow,  we  should  be  prepared  to  reevaluate 
our  current  hypotheses  for  both  modern  oceanography  and  past  variations. 
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Appendix  A:  All  radiocarbon  dates  generated  for  this  thesis. 
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upper  cal 

(median 

range  (1 

range  (1 

Depth 

C-14  Age 

probability) 

sigma 

sigma) 
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(cm) 

(yr  BP) 

Error  (yr) 

(yr  BP) 

(yr) 

(yr) 

BJ8-03 

13GGC 

32 

925 

30 

526 

547 

497 

100 

2190 

30 

1786 

1831 

1735 

272 

5050 

40 

5398 

5447 

5330 

408 

8220 

45 

8738 

8840 

8631 

505 

9610 

60 

10474 

10546 

10415 

BJ8-03 

23GGC 

7 

30 

25 

30 

73 

2040 

30 

1607 

1657 

1557 
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35 

3723 

3792 

3672 

329 

7160 

45 

7627 

7665 

7581 

409 

9730 

70 

10592 

10662 

10516 

BJ8-03 

70GGC 

64 

1310 

30 

853 

901 

816 

104 

2750 

50 

2472 

2540 

2352 
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4240 

50 

4341 

4413 

4269 

202 

4850 

40 

5155 

5245 

5079 

242 

5410 

35 

5789 

5847 

5737 

280 

6160 

60 

6594 

6666 

6513 
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45 

7653 
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7598 
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8080 

40 

8538 

8589 
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9438 

9325 
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9210 

45 
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394 

12550 

50 
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45 
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BJ8-03 
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25 

92 
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1 
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7440 

35 

7903 

7943 
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8230 

35 

8752 

8842 

8657 
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9740 

50 
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10643 
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55 

16758 

16859 

16675 

313 

29000 

210 

33006 

33358 

32688 

401 

37700 

290 

42063 

42299 

41822 

GeoB 10 

069-3 

8 

755 

35 

396 

449 

336 

65 

1680 

30 

1241 

1278 

1213 

88 

1970 

35 

1525 

1584 

1477 

107 

2260 

30 

1867 

1907 

1821 

143 


122 

2360 

30 

1982 

2032 

1931 

188 

2910 

35 

2698 

2738 

2670 

268 

3950 

35 

3939 

3990 

3872 

378 

5140 

35 

5511 

5557 

5473 

468 

6800 

45 

7328 

7380 

7279 

593 

9450 

65 

10307 

10378 

10221 

644 

10600 

50 

11882 

11970 

11747 

693 

11700 

55 

13194 

13250 

13132 

724 

13200 

75 

15215 

15481 

15007 

753 

16050 

75 

18789 

18867 

18701 

786 

19050 

90 

22267 

22391 

22159 

818 

21900 

110 

25723 

25968 

25524 

894 

27100 

160 

31177 

31268 

31085 

953 

31000 

240 

35043 

35228 

34766 
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Appendix  B:  Coretop  data  generated  for  this  thesis. 


Sediment 

core 

Latitude 

Longitude 

water 

depth 

[°N] 

[°E] 

[m] 

GeoB 

114  MC 

3.49 

95.33 

1535 

118  MC 

3.52 

96.31 

804 

89  MC 

2.78 

96.42 

916 

139  MC 

1.76 

96.77 

1854 

60  MC 

1.45 

98.05 

551 

147  MC 

-0.69 

98.07 

1052 

41  MC 

0.34 

98.13 

674 

10027-3 

-0.81 

99.65 

875 

31  MC 

-1.30 

99.72 

1734 

10033-3 

-1.56 

99.95 

1756 

10029-3 

-1.50 

100.13 

974 

10034-3 

-4.16 

101.50 

995 

03  MC 

-4.70 

101.96 

1707 

02  MC 

-5.48 

103.01 

1972 

10041-3 

-6.27 

103.09 

1540 

10038-3 

-5.94 

103.25 

1891 

10039-3 

-5.87 

103.29 

1799 

10036-3 

-5.34 

103.66 

1502 

10042-2 

-7.11 

104.64 

2457 

10049-5 

-8.78 

110.50 

1288 

10058-1 

-8.68 

112.64 

1103 

10065-9 

-9.22 

118.89 

1284 

10069-4 

-9.60 

120.92 

1249 

BJ8-03 

MCI  6 

-7.32 

115.18 

409 

MC1 1 

-7.39 

115.21 

590 

MC84 

-1.40 

117.53 

401 

MC31 

-3.88 

119.45 

459 

MC125 

-5.89 

120.29 

427 

Mg/C  a 
ruber 

6180 

ruber 

Mg/Ca 

obliq 

S180  obliq 

[mmol/mol] 

[%o  PDB] 

[mmol/mol] 

[%o  PDB] 

4.88 

2.40 

5.21 

3.01 

4.86 

3.08 

5.30 

2.84 

4.92 

3.10 

5.07 

3.11 

4.79 

3.29 

5.48 

-3.13 

2.66 

-2.31 

4.59 

2.08 

4.96 

-3.18 

-2.34 

5.31 

-3.27 

2.51 

-2.46 

4.67 

-3.00 

2.40 

-2.14 

4.64 

4.90 

1.97 

4.69 

-2.68 

2.61 

-2.15 

4.96 

-3.01 

-2.18 

5.03 

-3.19 

2.23 

-2.34 

5.10 

-3.23 

2.32 

-2.35 

4.43 

-3.19 

-2.34 

4.55 

-2.87 

1.82 

-2.03 

4.44 

-3.09 

2.04 

-1.85 

4.37 

-2.91 

2.31 

-1.62 

4.34 

-2.92 

1.92 

-2.17 

5.06 

-2.87 

2.55 

-1.96 

4.89 

-2.87 

2.86 

-1.67 

5.27 

-3.17 

2.79 

5.14 

-2.78 

2.80 

-2.27 

4.94 

-2.81 

3.15 

-1.85 
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Appendix  C:  Downcore  data  generated  for  this  thesis. 


GeoB10069-3 

G.  ruber 

Mg/Ca 

Depth 

5lsO 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

3 

4.10 

0 

-2.760 

8 

4.68 

1 

-2.870 

13 

3.84 

2 

-2.990 

18 

3.72 

3 

-1.920 

23 

3.85 

4 

-3.160 

28 

3.83 

5 

-2.930 

33 

3.77 

6 

-2.770 

38 

4.28 

7 

-2.690 

43 

3.81 

8 

-2.980 

48 

3.99 

9 

-2.360 

53 

3.93 

10 

-2.690 

58 

4.04 

11 

-2.580 

63 

4.36 

12 

-2.760 

68 

3.93 

13 

-2.880 

73 

4.07 

14 

-2.600 

78 

4.03 

15 

-2.780 

83 

4.09 

16 

-2.880 

88 

3.70 

17 

-2.790 

93 

3.83 

18 

-2.760 

98 

3.89 

19 

-2.850 

103 

4.12 

20 

-2.720 

108 

4.01 

21 

-2.810 

113 

4.10 

22 

-2.540 

118 

4.17 

23 

-2.790 

123 

3.99 

24 

-2.870 

128 

4.28 

25 

-2.860 

133 

4.11 

26 

-2.730 

138 

4.27 

27 

-2.790 

143 

3.81 

28 

-2.460 

148 

4.16 

29 

-2.660 

153 

3.93 

30 

-2.500 

158 

3.87 

31 

-2.610 

163 

4.12 

32 

-2.650 

168 

3.84 

33 

-2.750 

173 

3.87 

34 

-2.760 

178 

3.81 

35 

-2.690 

183 

3.84 

36 

-2.700 

188 

4.28 

37 

-2.710 

193 

4.13 

38 

-2.540 

198 

3.90 

39 

-2.740 

203 

3.95 

40 

-2.720 

208 

4.11 

41 

-2.610 

213 

3.86 

42 

-2.710 

218 

4.26 

43 

-2.610 
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223 

4.34 

44 

-2.680 

228 

3.83 

45 

-2.620 

233 

4.06 

46 

-2.580 

238 

3.96 

47 

-2.710 

243 

4.02 

48 

-2.710 

248 

4.08 

49 

-2.720 

253 

4.13 

50 

-2.720 

258 

4.21 

51 

-2.740 

263 

3.89 

52 

-2.600 

268 

3.93 

53 

-2.780 

273 

4.16 

54 

-2.610 

278 

4.00 

55 

-2.510 

283 

4.09 

56 

-2.890 

288 

4.29 

57 

-2.830 

293 

4.39 

58 

-2.730 

298 

4.17 

59 

-2.740 

303 

4.19 

60 

-2.770 

308 

3.83 

61 

-2.630 

313 

4.23 

62 

-2.720 

318 

3.83 

63 

-2.580 

323 

4.12 

64 

-2.670 

328 

4.32 

65 

-2.870 

333 

4.09 

66 

-2.660 

338 

4.13 

67 

-2.660 

343 

4.16 

68 

-2.690 

348 

4.24 

69 

-2.420 

353 

4.29 

70 

-2.640 

358 

4.19 

71 

-2.280 

363 

4.18 

72 

-2.610 

368 

4.36 

73 

-2.690 

373 

4.20 

74 

-2.580 

378 

4.18 

75 

-2.560 

383 

4.07 

76 

-2.640 

388 

4.02 

77 

-2.820 

393 

4.19 

78 

-2.870 

398 

4.36 

79 

-2.580 

403 

4.13 

80 

-2.740 

408 

4.21 

81 

-2.790 

413 

4.43 

82 

-2.620 

418 

4.20 

83 

-2.660 

423 

4.16 

84 

-2.580 

428 

4.49 

85 

-2.570 

433 

4.31 

86 

-2.640 

438 

4.19 

87 

-2.640 

443 

4.17 

88 

-2.770 

449 

4.52 

89 

-2.800 

453 

4.17 

90 

-2.710 

458 

4.33 

91 

-2.670 

463 

4.09 

92 

-2.960 

148 


468 

4.27 

93 

-2.560 

473 

4.22 

94 

-2.880 

478 

4.16 

95 

-2.660 

483 

4.34 

96 

-2.520 

488 

4.29 

97 

-2.680 

493 

4.18 

98 

-2.410 

498 

3.93 

99 

-2.710 

503 

4.15 

100 

-2.730 

508 

4.28 

101 

-2.820 

513 

4.40 

102 

-2.710 

518 

4.34 

103 

-2.310 

523 

4.14 

104 

-2.630 

528 

4.14 

105 

-2.730 

533 

4.43 

106 

-2.770 

538 

4.16 

107 

-2.500 

543 

4.14 

108 

-2.790 

548 

4.31 

109 

-2.700 

553 

4.31 

110 

-2.460 

557 

4.40 

111 

-2.590 

563 

4.35 

112 

-2.750 

568 

4.43 

113 

-2.620 

573 

4.39 

114 

-2.880 

578 

4.37 

115 

-2.560 

583 

4.65 

116 

-2.920 

588 

4.72 

117 

-2.610 

593 

4.16 

118 

-2.640 

598 

4.84 

119 

-2.730 

602 

4.52 

120 

-2.540 

603 

4.59 

121 

-2.580 

603 

4.51 

122 

-2.660 

604 

4.35 

123 

-2.640 

606 

4.25 

124 

-2.680 

608 

4.28 

125 

-2.560 

608 

4.73 

126 

-2.800 

610 

4.76 

127 

-2.830 

612 

4.22 

128 

-2.920 

613 

4.63 

129 

-2.940 

613 

4.69 

130 

-2.670 

614 

4.11 

131 

-2.770 

616 

4.48 

132 

-2.540 

618 

4.65 

133 

-2.720 

618 

4.21 

138 

-2.610 

620 

4.33 

143 

-2.520 

622 

4.29 

148 

-2.630 

623 

4.35 

153 

-2.570 

623 

4.44 

158 

-2.840 

624 

4.17 

163 

-2.620 

626 

4.29 

168 

-2.540 

628 

4.34 

173 

-2.770 

149 


628 

4.10 

178 

-2.630 

630 

4.47 

183 

-2.660 

632 

4.43 

188 

-2.450 

633 

4.29 

193 

-2.740 

633 

4.41 

198 

-2.710 

634 

4.49 

203 

-2.630 

636 

4.70 

208 

-2.780 

638 

4.30 

218 

-2.730 

638 

4.52 

223 

-2.590 

640 

4.42 

228 

-2.500 

642 

4.66 

233 

-2.810 

643 

4.51 

238 

-2.650 

643 

4.28 

243 

-2.670 

644 

4.16 

248 

-2.830 

646 

4.32 

253 

-2.420 

648 

4.16 

258 

-2.820 

650 

4.39 

263 

-2.730 

652 

4.11 

268 

-2.860 

653 

4.08 

273 

-2.590 

654 

4.61 

278 

-2.790 

656 

4.29 

283 

-2.670 

658 

3.93 

288 

-2.690 

660 

3.94 

293 

-2.670 

662 

4.57 

298 

-2.780 

663 

3.80 

303 

-2.670 

664 

3.91 

308 

-2.930 

666 

3.99 

313 

-2.670 

668 

3.74 

318 

-2.900 

670 

4.20 

323 

-2.740 

672 

3.81 

328 

-2.810 

673 

3.86 

333 

-3.050 

674 

4.13 

338 

-2.790 

676 

4.09 

343 

-2.730 

678 

3.89 

348 

-2.730 

680 

4.38 

353 

-2.830 

682 

3.90 

358 

-2.610 

683 

3.95 

363 

-2.640 

684 

4.17 

368 

-2.700 

686 

3.91 

373 

-2.640 

688 

4.00 

378 

-2.390 

692 

3.92 

383 

-2.610 

693 

3.73 

388 

-2.660 

694 

3.84 

393 

-2.600 

696 

4.03 

398 

-2.570 

698 

4.18 

403 

-2.880 

700 

3.94 

408 

-2.830 

702 

3.99 

413 

-2.740 

703 

4.05 

418 

-2.780 

704 

3.76 

423 

-2.600 

150 


706 

3.80 

428 

-2.730 

708 

3.90 

433 

-2.660 

710 

3.87 

438 

-2.830 

712 

3.82 

443 

-2.600 

713 

3.83 

448 

-2.740 

714 

3.67 

453 

-2.610 

716 

3.90 

458 

-2.520 

718 

3.70 

463 

-2.590 

720 

3.70 

468 

-2.530 

111 

3.66 

473 

-3.020 

723 

3.93 

478 

-2.920 

724 

3.78 

483 

-2.520 

726 

3.71 

488 

-2.620 

728 

3.94 

493 

-2.430 

730 

3.49 

498 

-2.290 

732 

3.55 

503 

-2.930 

733 

3.57 

508 

-2.870 

734 

3.92 

513 

-2.560 

736 

3.80 

518 

-2.580 

738 

4.01 

523 

-2.650 

740 

3.72 

528 

-2.610 

742 

3.76 

533 

-2.550 

743 

3.41 

538 

-2.480 

744 

3.67 

543 

-2.610 

746 

3.50 

548 

-2.320 

748 

3.41 

553 

-2.870 

750 

3.61 

558 

-2.600 

752 

3.74 

563 

-2.380 

753 

3.34 

568 

-2.440 

754 

3.38 

573 

-2.430 

756 

3.53 

578 

-2.560 

758 

3.55 

583 

-2.500 

760 

3.71 

588 

-2.350 

762 

3.54 

593 

-2.640 

763 

3.63 

598 

-2.300 

764 

3.41 

602 

-2.340 

766 

3.53 

603 

-2.510 

768 

3.32 

604 

-2.340 

770 

3.66 

606 

-2.260 

772 

3.54 

608 

-2.350 

773 

3.39 

610 

-2.230 

774 

3.34 

612 

-2.320 

776 

3.21 

613 

-2.440 

780 

3.41 

614 

-2.130 

782 

3.65 

616 

-2.130 

783 

3.15 

618 

-2.090 

783 

3.48 

620 

-2.190 

784 

3.34 

622 

-2.320 

786 

3.70 

623 

-2.160 

151 


788 

3.44 

624 

-2.140 

790 

3.56 

626 

-2.160 

792 

3.57 

628 

-2.100 

793 

3.47 

630 

-2.150 

794 

3.57 

632 

-2.260 

796 

3.51 

633 

-1.900 

798 

3.29 

634 

-2.020 

800 

3.45 

636 

-2.110 

802 

3.42 

638 

-1.730 

803 

3.14 

640 

-2.030 

804 

3.60 

642 

-1.900 

806 

3.51 

643 

-1.900 

808 

3.40 

644 

-1.870 

810 

3.40 

646 

-1.960 

812 

3.62 

648 

-1.870 

813 

3.56 

650 

-1.720 

814 

3.61 

652 

-1.700 

816 

3.57 

653 

-1.720 

818 

3.63 

654 

-1.770 

820 

3.58 

656 

-1.660 

822 

3.74 

658 

-1.520 

823 

3.42 

660 

-1.530 

824 

3.74 

662 

-1.670 

826 

3.34 

663 

-1.630 

828 

3.39 

664 

-1.910 

830 

3.86 

666 

-1.630 

832 

3.62 

668 

-1.740 

833 

3.37 

670 

-1.920 

834 

3.26 

672 

-1.840 

836 

3.42 

673 

-1.770 

838 

3.18 

674 

-1.940 

840 

3.38 

676 

-1.920 

842 

3.47 

678 

-1.660 

843 

3.33 

680 

-1.620 

844 

3.37 

682 

-1.720 

846 

3.24 

683 

-1.660 

848 

3.49 

684 

-1.660 

850 

3.33 

686 

-1.430 

852 

3.57 

688 

-1.550 

853 

3.56 

690 

-1.570 

854 

3.60 

692 

-1.890 

856 

3.43 

693 

-1.780 

858 

3.69 

694 

-1.620 

860 

3.48 

696 

-2.020 

862 

3.19 

698 

-1.440 

863 

3.49 

700 

-1.700 

864 

3.22 

702 

-1.830 

866 

3.26 

703 

-1.580 

868 

3.44 

704 

-1.810 

152 


870 

3.48 

706 

-1.750 

872 

3.63 

708 

-1.510 

873 

3.60 

710 

-1.800 

874 

3.51 

712 

-1.420 

876 

3.52 

713 

-1.650 

878 

3.55 

714 

-1.350 

880 

3.54 

716 

-1.500 

882 

3.48 

718 

-1.360 

883 

3.56 

720 

-1.420 

884 

3.36 

722 

-1.360 

888 

3.51 

723 

-1.620 

890 

3.34 

724 

-1.560 

892 

3.40 

726 

-1.360 

893 

3.28 

728 

-1.310 

894 

3.44 

730 

-1.070 

896 

3.48 

732 

-0.960 

898 

3.63 

733 

-1.200 

900 

3.44 

734 

-1.240 

902 

3.39 

736 

-1.480 

903 

3.54 

738 

-1.210 

904 

3.52 

740 

-1.140 

906 

3.33 

742 

-1.070 

908 

3.26 

743 

-1.100 

910 

3.62 

744 

-1.120 

912 

3.51 

746 

-1.060 

913 

3.51 

748 

-0.930 

914 

3.34 

750 

-1.070 

916 

3.37 

752 

-0.950 

918 

3.34 

753 

-0.950 

920 

3.49 

754 

-0.990 

922 

3.36 

756 

-0.800 

923 

3.34 

758 

-1.020 

924 

3.45 

760 

-0.820 

926 

3.09 

762 

-1.220 

928 

3.34 

763 

-1.170 

930 

3.59 

764 

-1.100 

932 

3.57 

766 

-1.050 

933 

3.35 

768 

-1.110 

934 

3.47 

770 

-0.840 

936 

3.25 

772 

-1.220 

938 

3.48 

773 

-1.310 

940 

3.26 

774 

-1.480 

942 

3.13 

776 

-0.930 

943 

3.43 

778 

-1.330 

944 

3.28 

780 

-1.100 

946 

3.43 

782 

-1.240 

948 

3.63 

783 

-1.010 

950 

3.45 

784 

-1.340 

952 

3.51 

786 

-1.200 

153 


953 

3.49 

788 

954 

3.81 

790 

956 

3.56 

792 

793 

794 
796 
798 
800 
802 

803 

804 
806 
808 
810 
812 

813 

814 
816 
818 
822 

823 

824 
826 
828 
830 

832 

833 

834 
836 
838 
840 

842 

843 

844 
846 
848 
850 

852 

853 

854 
856 
858 
860 
862 

863 

864 
866 
868 
870 


-0.930 

-1.400 

-1.010 

-1.120 

-0.990 

-1.040 

-1.040 

-1.190 

-1.060 

-1.060 

-1.170 

-1.100 

-1.240 

-1.220 

-1.090 

-0.880 

-1.420 

-1.250 

-1.410 

-1.220 

-1.230 

-0.910 

-1.250 

-1.550 

-0.900 

-0.850 

-1.080 

-1.260 

-1.210 

-1.210 

-1.140 

-1.040 

-1.440 

-1.130 

-1.230 

-1.150 

-1.150 

-1.190 

-0.980 

-1.090 

-1.320 

-1.150 

-1.370 

-1.430 

-1.380 

-1.280 

-1.360 

-1.200 

-1.440 


154 


872 

-1.360 

873 

-1.330 

874 

-1.360 

876 

-1.500 

878 

-1.340 

880 

-1.540 

882 

-0.980 

883 

-1.330 

884 

-1.370 

886 

-1.380 

888 

-1.420 

890 

-1.370 

892 

-1.250 

893 

-1.400 

894 

-1.210 

896 

-1.350 

898 

-1.140 

900 

-1.130 

902 

-1.230 

903 

-0.900 

904 

-1.190 

906 

-1.130 

908 

-1.210 

910 

-1.350 

912 

-1.100 

913 

-1.320 

914 

-1.120 

916 

-1.100 

918 

-1.390 

920 

-1.600 

922 

-1.320 

923 

-1.460 

924 

-1.460 

926 

-1.340 

928 

-1.260 

930 

-1.350 

932 

-1.330 

933 

-1.470 

934 

-1.260 

936 

-1.500 

938 

-1.300 

940 

-1.570 

942 

-1.320 

943 

-1.410 

944 

-1.350 

946 

-1.330 

948 

-1.610 

950 

-1.660 

952 

-1.370 

155 


953 

-1.630 

954 

-1.280 

956 

-1.450 

GeoB10069-3 

P.  obliquiloculata 

5lsO 

Mg/Ca 

Depth 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

3 

4.10 

1 

-2.037 

8 

4.68 

2 

-1.691 

13 

3.84 

3 

-1.977 

18 

3.72 

4 

-1.847 

23 

3.85 

6 

-2.061 

28 

3.83 

7 

-1.548 

33 

3.77 

8 

-1.935 

38 

4.28 

9 

-1.666 

43 

3.81 

11 

-1.727 

48 

3.99 

13 

-1.779 

53 

3.93 

14 

-1.388 

58 

4.04 

15 

-1.956 

63 

4.36 

16 

-1.934 

68 

3.93 

17 

-1.879 

73 

4.07 

18 

-1.864 

78 

4.03 

19 

-1.748 

83 

4.09 

20 

-2.506 

88 

3.70 

21 

-1.413 

93 

3.83 

22 

-1.531 

98 

3.89 

23 

-1.990 

103 

4.12 

24 

-1.825 

108 

4.01 

25 

-1.426 

113 

4.10 

26 

-1.647 

118 

4.17 

27 

-1.530 

123 

3.99 

28 

-1.558 

128 

4.28 

29 

-1.563 

133 

4.11 

30 

-1.331 

138 

4.27 

31 

-1.378 

143 

3.81 

32 

-1.817 

148 

4.16 

33 

-1.375 

153 

3.93 

34 

-1.516 

158 

3.87 

35 

-1.000 

163 

4.12 

36 

-1.434 

173 

3.87 

37 

-1.229 

178 

3.81 

38 

-1.717 

183 

3.84 

40 

-1.227 

188 

4.28 

41 

-1.533 

193 

4.13 

42 

-1.311 

198 

3.90 

43 

-1.386 

203 

3.95 

45 

-1.216 

208 

4.11 

46 

-0.724 

156 


213 

3.86 

47 

-1.165 

218 

4.26 

48 

-1.586 

223 

4.34 

49 

-1.379 

228 

3.83 

50 

-1.327 

233 

4.06 

51 

-1.404 

238 

3.96 

52 

-1.503 

243 

4.02 

53 

-1.265 

248 

4.08 

54 

-1.413 

253 

4.13 

55 

-1.414 

258 

4.21 

56 

-2.154 

263 

3.89 

57 

-0.862 

268 

3.93 

58 

-2.047 

273 

4.16 

59 

-1.570 

278 

4.00 

60 

-1.312 

283 

4.09 

61 

-2.136 

288 

4.29 

62 

-1.694 

293 

4.39 

63 

-1.725 

298 

4.17 

65 

-1.464 

303 

4.19 

67 

-1.441 

308 

3.83 

68 

-1.729 

313 

4.23 

70 

-1.379 

318 

3.83 

70 

-1.388 

323 

4.12 

71 

-1.467 

328 

4.32 

72 

-1.447 

333 

4.09 

73 

-1.784 

338 

4.13 

74 

-1.746 

343 

4.16 

75 

-1.553 

348 

4.24 

76 

-2.040 

353 

4.29 

77 

-1.584 

358 

4.19 

78 

-2.094 

363 

4.18 

79 

-1.629 

368 

4.36 

80 

-1.730 

373 

4.20 

81 

-1.454 

378 

4.18 

82 

-1.705 

383 

4.07 

83 

-1.671 

388 

4.02 

84 

-1.785 

393 

4.19 

85 

-1.635 

398 

4.36 

86 

-1.879 

403 

4.13 

87 

-1.611 

408 

4.21 

88 

-1.484 

413 

4.43 

89 

-1.214 

418 

4.20 

90 

-1.954 

423 

4.16 

91 

-1.349 

428 

4.49 

92 

-1.501 

433 

4.31 

93 

-1.623 

438 

4.19 

94 

-1.848 

443 

4.17 

95 

-1.679 

449 

4.52 

96 

-1.453 

453 

4.17 

97 

-1.651 

157 


458 

4.33 

98 

-1.757 

463 

4.09 

99 

-1.385 

468 

4.27 

100 

-1.738 

473 

4.22 

101 

-1.473 

478 

4.16 

102 

-1.521 

483 

4.34 

103 

-1.475 

488 

4.29 

104 

-1.798 

493 

4.18 

105 

-1.949 

498 

3.93 

106 

-1.563 

503 

4.15 

107 

-1.842 

508 

4.28 

108 

-1.815 

513 

4.40 

109 

-1.794 

518 

4.34 

110 

-1.758 

523 

4.14 

111 

-1.378 

528 

4.14 

112 

-1.806 

533 

4.43 

113 

-1.931 

538 

4.16 

114 

-1.725 

543 

4.14 

115 

-1.588 

548 

4.31 

116 

-2.145 

553 

4.31 

117 

-1.510 

563 

4.35 

118 

-1.975 

568 

4.43 

119 

-1.979 

573 

4.39 

120 

-1.824 

578 

4.37 

121 

-1.569 

583 

4.65 

122 

-2.015 

588 

4.72 

123 

-1.925 

593 

4.16 

124 

-2.134 

598 

4.84 

125 

-2.053 

602 

4.52 

126 

-1.812 

603 

4.55 

127 

-1.629 

604 

4.35 

128 

-1.815 

606 

4.25 

129 

-1.797 

610 

4.76 

130 

-2.265 

612 

4.22 

131 

-1.962 

613 

4.66 

132 

-1.823 

614 

4.11 

133 

-1.856 

616 

4.48 

138 

-1.812 

618 

4.43 

143 

-1.542 

620 

4.33 

148 

-1.549 

622 

4.29 

153 

-1.583 

623 

4.40 

158 

-1.662 

624 

4.17 

163 

-1.790 

626 

4.29 

168 

-1.266 

628 

4.22 

173 

-1.528 

630 

4.47 

178 

-1.650 

632 

4.43 

183 

-1.589 

633 

4.35 

188 

-1.672 

634 

4.49 

193 

-1.553 

636 

4.70 

198 

-1.674 

158 


638 

4.41 

203 

-1.494 

640 

4.42 

208 

-1.657 

642 

4.66 

213 

-1.518 

643 

4.40 

218 

-1.536 

644 

4.16 

223 

-1.641 

646 

4.32 

228 

-1.590 

648 

4.16 

233 

-1.558 

650 

4.39 

238 

-1.147 

652 

4.11 

243 

-1.443 

653 

4.08 

248 

-1.315 

654 

4.61 

253 

-1.603 

656 

4.29 

258 

-1.511 

658 

3.93 

263 

-1.420 

660 

3.94 

268 

-1.547 

662 

4.57 

273 

-1.786 

663 

3.80 

278 

-1.765 

664 

3.91 

283 

-1.569 

666 

3.99 

288 

-1.423 

668 

3.74 

293 

-1.448 

670 

4.20 

298 

-1.355 

672 

3.81 

303 

-1.551 

673 

3.86 

308 

-1.161 

674 

4.13 

313 

-1.687 

676 

4.09 

318 

-2.002 

678 

3.89 

323 

-1.488 

680 

4.38 

328 

-1.467 

682 

3.90 

333 

-1.721 

683 

3.95 

338 

-1.806 

684 

4.17 

343 

-1.699 

686 

3.91 

348 

-1.763 

688 

4.00 

353 

-1.933 

692 

3.92 

358 

-1.591 

693 

3.73 

363 

-1.831 

694 

3.84 

368 

-1.642 

696 

4.03 

373 

-1.385 

698 

4.18 

378 

-1.285 

700 

3.94 

383 

-1.298 

702 

3.99 

388 

-1.801 

703 

4.05 

393 

-1.549 

704 

3.76 

398 

-1.658 

706 

3.80 

403 

-1.365 

708 

3.90 

408 

-1.584 

710 

3.87 

413 

-1.892 

712 

3.82 

418 

-1.576 

713 

3.83 

423 

-1.518 

714 

3.67 

428 

-1.517 

716 

3.90 

433 

-1.709 

718 

3.70 

438 

-1.719 

720 

3.70 

443 

-1.482 
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722 

3.66 

449 

-1.830 

723 

3.93 

453 

-1.757 

724 

3.78 

458 

-1.523 

726 

3.71 

463 

-1.592 

728 

3.94 

468 

-1.647 

730 

3.49 

473 

-1.562 

732 

3.55 

478 

-1.977 

733 

3.57 

483 

-1.700 

734 

3.92 

488 

-1.449 

736 

3.80 

493 

-1.461 

738 

4.01 

498 

-1.545 

740 

3.72 

503 

-1.460 

742 

3.76 

508 

-1.446 

743 

3.41 

513 

-1.558 

744 

3.67 

518 

-1.455 

746 

3.50 

523 

-1.473 

748 

3.41 

528 

-1.721 

750 

3.61 

533 

-1.523 

752 

3.74 

538 

-1.349 

753 

3.34 

543 

-1.617 

754 

3.38 

548 

-1.426 

756 

3.53 

553 

-1.403 

758 

3.55 

557 

-1.390 

760 

3.71 

563 

-1.032 

762 

3.54 

568 

-1.321 

763 

3.63 

573 

-1.359 

764 

3.41 

578 

-1.601 

766 

3.53 

583 

-1.203 

768 

3.32 

588 

-1.325 

770 

3.66 

593 

-1.516 

772 

3.54 

598 

-1.144 

773 

3.39 

602 

-0.958 

774 

3.34 

603 

-0.938 

776 

3.21 

604 

-0.906 

780 

3.41 

606 

-1.033 

782 

3.65 

608 

-1.177 

783 

3.15 

610 

-1.154 

783 

3.48 

612 

-0.990 

784 

3.34 

613 

-1.061 

786 

3.70 

614 

-0.784 

788 

3.44 

616 

-1.115 

790 

3.56 

618 

-0.931 

792 

3.57 

620 

-1.228 

793 

3.47 

622 

-0.760 

794 

3.57 

623 

-0.879 

796 

3.51 

626 

-1.095 

798 

3.29 

628 

-0.597 

800 

3.45 

630 

-1.003 

802 

3.42 

632 

-0.794 
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803 

3.14 

633 

-0.625 

804 

3.60 

634 

-0.675 

806 

3.51 

636 

-0.743 

808 

3.40 

638 

-0.710 

810 

3.40 

640 

-0.208 

812 

3.62 

642 

-0.322 

813 

3.56 

643 

-0.185 

814 

3.61 

644 

-0.253 

816 

3.57 

646 

-0.711 

818 

3.63 

648 

-0.246 

820 

3.58 

650 

-0.086 

822 

3.74 

652 

-0.278 

823 

3.42 

653 

-0.025 

824 

3.74 

654 

-0.252 

826 

3.34 

656 

-0.226 

828 

3.39 

658 

-0.411 

830 

3.86 

660 

-0.056 

832 

3.62 

662 

-0.262 

833 

3.37 

663 

-0.018 

834 

3.26 

664 

-0.155 

836 

3.42 

666 

-0.216 

838 

3.18 

668 

-0.310 

840 

3.38 

670 

-0.212 

842 

3.47 

672 

-0.419 

843 

3.33 

673 

-0.225 

844 

3.37 

674 

-0.134 

846 

3.24 

676 

-0.522 

848 

3.49 

678 

-0.291 

850 

3.33 

683 

-0.172 

852 

3.57 

688 

-0.174 

853 

3.56 

693 

-0.361 

854 

3.60 

698 

-0.448 

856 

3.43 

703 

-0.212 

858 

3.69 

708 

0.345 

860 

3.48 

713 

-0.011 

862 

3.19 

716 

-0.345 

863 

3.49 

718 

0.148 

864 

3.22 

720 

0.080 

866 

3.26 

722 

-0.074 

868 

3.44 

723 

0.102 

870 

3.48 

724 

0.104 

872 

3.63 

728 

0.239 

873 

3.60 

730 

-0.163 

874 

3.51 

732 

0.118 

876 

3.52 

733 

0.031 

878 

3.55 

734 

0.299 

880 

3.54 

738 

0.242 

882 

3.48 

740 

0.090 

883 

3.56 

742 

-0.164 
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884 

3.36 

743 

0.285 

888 

3.51 

744 

0.351 

890 

3.34 

746 

-0.306 

892 

3.40 

748 

-0.114 

893 

3.28 

750 

0.353 

894 

3.44 

752 

0.501 

896 

3.48 

753 

0.469 

898 

3.63 

754 

0.258 

900 

3.44 

756 

0.181 

902 

3.39 

758 

0.283 

903 

3.54 

760 

0.199 

904 

3.52 

762 

0.213 

906 

3.33 

763 

-0.062 

908 

3.26 

764 

-0.025 

910 

3.62 

768 

0.158 

912 

3.51 

773 

-0.065 

913 

3.51 

783 

0.001 

914 

3.34 

783 

-0.116 

916 

3.37 

788 

0.167 

918 

3.34 

793 

-0.030 

920 

3.49 

798 

0.003 

922 

3.36 

803 

-0.153 

923 

3.34 

808 

-0.009 

928 

3.34 

813 

0.103 

933 

3.35 

818 

-0.051 

938 

3.48 

823 

-0.121 

943 

3.43 

828 

0.151 

948 

3.63 

833 

-0.069 

953 

3.49 

838 

0.117 

843 

0.296 

848 

0.121 

853 

0.247 

858 

0.020 

863 

-0.205 

868 

-0.250 

873 

-0.052 

878 

-0.148 

883 

0.300 

888 

-0.069 

893 

-0.049 

898 

0.124 

903 

0.211 

908 

-0.162 

913 

-0.030 

918 

-0.015 

923 

-0.198 

928 

-0.129 

933 

-0.083 

938 

-0.006 
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943 

-0.140 

948 

-0.365 

953 

0.046 

BJ8-03 

P. 

70GGC 

obliquilocula  ta 

5lsO 

Mg/Ca 

Depth 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

7.5 

2.51 

7.5 

-2.168 

9.5 

2.98 

8.5 

-1.746 

13.5 

2.85 

13.5 

-1.992 

15.5 

2.37 

21.5 

-2.143 

17.5 

2.75 

29.5 

-1.987 

21.5 

2.37 

15.5 

-1.969 

23.5 

2.56 

23.5 

-1.953 

24.5 

2.97 

24.5 

-2.269 

25.5 

2.44 

31.5 

-2.172 

29.5 

3.08 

32.5 

-1.948 

31.5 

2.75 

37.5 

-2.365 

32.5 

2.37 

39.5 

-2.043 

33.5 

2.91 

45.5 

-2.291 

37.5 

2.60 

47.5 

-2.244 

39.5 

2.65 

48.5 

-2.137 

41.5 

2.75 

53.5 

-2.070 

45.5 

2.78 

56.5 

-1.869 

47.5 

2.72 

61.5 

-1.267 

48.5 

2.68 

69.5 

-2.025 

49.5 

2.60 

72.5 

-2.019 

53.5 

2.86 

77.5 

-2.210 

56.5 

2.89 

80.5 

-2.148 

57.5 

2.77 

85.5 

-2.053 

61.5 

2.80 

88.5 

-2.091 

65.5 

2.67 

93.5 

-1.923 

69.5 

2.59 

101.5 

-1.933 

72.5 

2.49 

109.5 

-1.724 

73.5 

2.62 

125.5 

-1.572 

77.5 

2.35 

133.5 

-2.028 

80.5 

2.94 

136.5 

-2.068 

81.5 

2.89 

141.5 

-1.854 

85.5 

2.62 

144.5 

-2.096 

88.5 

2.45 

149.5 

-1.912 

89.5 

3.18 

152.5 

-1.726 

93.5 

2.77 

157.5 

-2.153 

97.5 

2.55 

160.5 

-2.016 

101.5 

2.54 

165.5 

-2.124 

105.5 

2.73 

168.5 

-1.922 

109.5 

2.83 

173.5 

-1.939 

113.5 

2.84 

181.5 

-1.985 

117.5 

2.83 

189.5 

-2.083 
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121.5 

2.55 

192.5 

-1.848 

125.5 

2.73 

197.5 

-1.642 

129.5 

2.49 

205.5 

-1.989 

133.5 

2.80 

208.5 

-1.854 

136.5 

2.82 

213.5 

-2.036 

137.5 

2.80 

221.5 

-1.677 

141.5 

2.74 

229.5 

-1.999 

144.5 

2.70 

232.5 

-1.991 

145.5 

2.73 

237.5 

-1.735 

149.5 

3.38 

240.5 

-1.765 

152.5 

2.73 

245.5 

-2.102 

153.5 

2.81 

248.5 

-1.914 

157.5 

2.87 

253.5 

-1.962 

160.5 

2.65 

256.5 

-1.643 

161.5 

2.67 

261.5 

-1.864 

165.5 

2.97 

264.5 

-1.969 

168.5 

2.90 

269.5 

-2.125 

169.5 

2.68 

272.5 

-1.926 

173.5 

2.93 

277.5 

-2.064 

177.5 

2.73 

285.5 

-2.096 

181.5 

2.66 

288.5 

-1.889 

185.5 

3.00 

293.5 

-1.946 

189.5 

2.82 

296.5 

-1.803 

192.5 

2.63 

300.5 

-1.826 

193.5 

2.90 

301.5 

-2.008 

197.5 

2.60 

302.5 

-2.007 

201.5 

3.06 

303.5 

-1.791 

205.5 

2.70 

306.5 

-1.824 

208.5 

2.69 

307.5 

-1.715 

209.5 

2.93 

308.5 

-1.790 

213.5 

3.18 

310.5 

-1.609 

217.5 

2.93 

311.5 

-1.904 

221.5 

2.72 

314.5 

-1.682 

225.5 

3.11 

315.5 

-1.715 

229.5 

2.95 

316.5 

-1.699 

232.5 

2.98 

318.5 

-1.628 

233.5 

2.93 

319.5 

-1.891 

237.5 

2.93 

322.5 

-1.895 

240.5 

2.57 

323.5 

-1.829 

245.5 

2.64 

324.5 

-1.666 

248.5 

2.80 

326.5 

-2.090 

249.5 

2.83 

327.5 

-1.676 

253.5 

2.71 

330.5 

-1.814 

256.5 

2.80 

331.5 

-1.714 

257.5 

3.38 

332.5 

-1.718 

261.5 

2.60 

334.5 

-1.668 

264.5 

2.88 

335.5 

-1.684 

265.5 

3.11 

338.5 

-1.724 

269.5 

2.86 

339.5 

-1.456 
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272.5 

2.74 

340.5 

-1.388 

273.5 

2.85 

342.5 

-1.599 

277.5 

2.71 

343.5 

-1.520 

281.5 

2.77 

344.5 

-1.533 

285.5 

2.71 

346.5 

-1.520 

288.5 

2.44 

347.5 

-1.821 

289.5 

2.90 

348.5 

-1.418 

293.5 

2.86 

349.5 

-1.443 

296.5 

2.71 

351.5 

-1.782 

297.5 

3.13 

352.5 

-1.564 

300.5 

2.75 

354.5 

-1.256 

301.5 

2.73 

355.5 

-1.225 

302.5 

2.92 

356.5 

-0.805 

303.5 

2.91 

358.5 

-1.432 

305 

3.02 

359.5 

-1.122 

306.5 

2.78 

360.5 

-1.005 

307.5 

2.73 

362.5 

-1.042 

308.5 

2.88 

363.5 

-1.109 

310.5 

2.63 

364.5 

-1.037 

311.5 

2.77 

366.5 

-0.968 

314.5 

2.86 

367.5 

-0.840 

315.5 

2.80 

368.5 

-0.748 

316.5 

2.88 

370.5 

-0.645 

318.5 

3.02 

371.5 

-0.606 

319.5 

2.98 

372.5 

-0.605 

322.5 

2.69 

374.5 

-0.597 

323.5 

2.70 

375.5 

-0.510 

324.5 

2.84 

376.5 

-0.399 

326.5 

2.86 

377.5 

-0.613 

327.5 

3.11 

379.5 

-0.791 

329.5 

2.88 

380.5 

-0.683 

330.5 

2.78 

382.5 

-0.649 

331.5 

2.61 

383.5 

-0.769 

332.5 

2.74 

384.5 

-0.620 

334.5 

2.82 

386.5 

-0.559 

335.5 

2.88 

387.5 

-0.504 

336.5 

2.89 

388.5 

-0.476 

338.5 

2.89 

389.5 

-0.636 

339.5 

2.89 

391.5 

-0.516 

340.5 

2.68 

392.5 

-0.360 

342.5 

2.80 

394.5 

-0.450 

343.5 

2.74 

395.5 

-0.365 

344.5 

3.00 

346.5 

-0.439 

345.5 

2.94 

398.5 

-0.443 

346.5 

2.81 

399.5 

-0.352 

347.5 

2.90 

400.5 

-0.583 

348.5 

2.90 

402.5 

-0.459 

349.5 

3.01 

403.5 

-0.321 

351.5 

2.91 

404.5 

-0.545 
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352.5 

2.75 

353.5 

3.20 

354.5 

2.93 

355.5 

2.98 

356.5 

2.83 

358.5 

2.97 

359.5 

2.92 

360.5 

2.71 

361.5 

2.80 

362.5 

2.60 

366.5 

2.50 

367.5 

2.82 

368.5 

2.70 

369.5 

2.50 

370.5 

2.71 

371.5 

2.28 

372.5 

2.40 

374.5 

2.69 

375.5 

2.55 

376.5 

2.44 

377.5 

2.41 

377.5 

2.82 

379.5 

2.67 

380.5 

2.35 

382.5 

2.72 

383.5 

2.84 

384.5 

2.63 

385.5 

2.62 

386.5 

2.56 

387.5 

2.50 

388.5 

2.45 

389.5 

2.49 

391.5 

2.41 

392.5 

2.43 

393.5 

2.73 

394.5 

2.36 

395.5 

2.34 

396.5 

2.49 

398.5 

2.79 

399.5 

2.56 

400.5 

2.46 

401.5 

2.48 

402.5 

2.45 

403.5 

2.69 

404.5 

2.37 

406.5 

2.38 

BJ8-03 

23GGC  G.  ruber 
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Mg/Ca 

Depth 

5lsO 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

5.5 

4.78 

5.5 

-3.392 

8.5 

4.38 

5.5 

-3.477 

24.5 

4.46 

6.5 

-3.049 

40.5 

4.68 

6.5 

-3.149 

56.5 

4.68 

8.5 

-3.316 

72.5 

4.56 

8.5 

-3.281 

88.5 

4.48 

16.5 

-3.143 

104.5 

4.58 

16.5 

-3.154 

120.5 

4.71 

24.5 

-3.030 

136.5 

4.69 

24.5 

-2.915 

152.5 

4.58 

32.5 

-3.053 

168.5 

4.43 

32.5 

-3.142 

200.5 

4.50 

40.5 

-2.761 

216.5 

4.57 

40.5 

-2.901 

232.5 

4.55 

48.5 

-3.075 

248.5 

4.80 

48.5 

-3.217 

264.5 

4.66 

56.5 

-3.178 

280.5 

4.67 

56.5 

-3.058 

296.5 

4.61 

64.5 

-2.787 

312.5 

4.69 

64.5 

-3.025 

328.5 

5.11 

72.5 

-3.240 

344.5 

4.70 

72.5 

-3.188 

360.5 

4.57 

80.5 

-3.179 

376.5 

4.87 

80.5 

-3.275 

392.5 

4.69 

88.5 

-3.019 

408.5 

4.60 

88.5 

-3.031 

96.5 

-3.019 

96.5 

-3.141 

104.5 

-3.058 

104.5 

-3.108 

112.5 

-3.062 

112.5 

-3.232 

120.5 

-3.004 

120.5 

-3.270 

128.5 

-3.053 

128.5 

-2.710 

136.5 

-3.082 

136.5 

-3.238 

144.5 

-3.043 

144.5 

-3.047 

152.5 

-3.074 

152.5 

-3.412 

160.5 

-3.263 

160.5 

-3.409 

168.5 

-2.890 

168.5 

-2.883 

176.5 

-3.289 
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176.5 

-2.899 

184.5 

-2.868 

184.5 

-3.174 

192.5 

-3.216 

192.5 

-2.976 

200.5 

-2.996 

200.5 

-3.243 

208.5 

-3.184 

208.5 

-3.302 

216.5 

-2.888 

216.5 

-3.291 

224.5 

-3.363 

224.5 

-3.197 

232.5 

-3.228 

232.5 

-3.089 

240.5 

-2.941 

240.5 

-3.241 

248.5 

-3.196 

248.5 

-3.325 

256.5 

-3.294 

256.5 

-2.965 

264.5 

-2.988 

264.5 

-3.139 

272.5 

-3.346 

272.5 

-3.040 

280.5 

-3.172 

280.5 

-3.080 

288.5 

-3.023 

288.5 

-2.961 

296.5 

-3.098 

296.5 

-2.659 

304.5 

-3.068 

304.5 

-3.031 

312.5 

-3.150 

312.5 

-2.835 

324 

-2.832 

324 

-2.853 

328.5 

-2.941 

328.5 

-2.841 

336.5 

-2.895 

336.5 

-3.168 

344.5 

-3.017 

344.5 

-2.717 

352.5 

-2.827 

352.5 

-2.711 

360.5 

-2.953 

360.5 

-2.601 

368.5 

-3.120 

368.5 

-3.251 

168 


376.5 

-2.676 

376.5 

-2.763 

384.5 

-2.779 

384.5 

-2.608 

392.5 

-2.617 

392.5 

-2.927 

400.5 

-2.603 

400.5 

-2.919 

408.5 

-2.398 

408.5 

-2.753 

416.5 

-2.498 

416.5 

-2.432 

BJ8-03 

P. 

23GGC 

obliquilocula  ta 

5lsO 

Mg/Ca 

Depth 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

16.5 

2.95 

8.5 

-2.497 

24.5 

2.93 

16.5 

-2.394 

32.5 

2.66 

24.5 

-2.004 

40.5 

2.45 

32.5 

-2.268 

48.5 

2.97 

40.5 

-2.237 

56.5 

2.60 

48.5 

-2.088 

64.5 

2.59 

56.5 

-2.125 

72.5 

2.67 

64.5 

-2.492 

80.5 

2.79 

72.5 

-2.353 

88.5 

2.79 

80.5 

-2.418 

96.5 

2.72 

88.5 

-1.938 

104.5 

2.76 

96.5 

-2.028 

120.5 

2.61 

104.5 

-2.144 

128.5 

2.64 

112.5 

-2.119 

144.5 

2.78 

120.5 

-2.021 

152.5 

2.69 

128.5 

-1.958 

160.5 

2.89 

136.5 

-2.103 

168.5 

2.63 

144.5 

-2.023 

176.5 

2.82 

152.5 

-1.941 

184.5 

3.27 

160.5 

-2.060 

192.5 

2.66 

168.5 

-1.923 

200.5 

2.77 

176.5 

-1.905 

208.5 

2.63 

184.5 

-2.243 

216.5 

2.68 

192.5 

-1.966 

224.5 

2.71 

200.5 

-2.019 

232.5 

2.90 

208.5 

-2.014 

240.5 

2.82 

216.5 

-2.509 

248.5 

2.72 

224.5 

-2.128 

264.5 

2.80 

232.5 

-2.118 

272.5 

2.54 

240.5 

-2.162 

280.5 

2.89 

248.5 

-1.601 

288.5 

2.82 

256.5 

-2.206 

296.5 

2.69 

264.5 

-1.819 

304.5 

2.74 

272.5 

-2.206 

169 

312.5 

2.90 

280.5 

-2.759 

320.5 

2.93 

288.5 

-2.097 

328.5 

3.91 

296.5 

-2.091 

336.5 

2.99 

304.5 

-1.734 

344.5 

2.61 

312.5 

-1.754 

360.5 

2.67 

328.5 

-1.795 

368.5 

2.69 

336.5 

-1.985 

384.5 

2.78 

344.5 

-1.896 

392.5 

2.78 

352.5 

-1.772 

400.5 

2.96 

360.5 

-1.514 

408.5 

2.86 

368.5 

-1.897 

416.5 

3.16 

376.5 

-2.034 

384.5 

-1.448 

392.5 

-1.724 

400.5 

-1.346 

408.5 

-1.458 

416.5 

-1.355 

P. 

MD97-2141 

obliquilocula  ta 

5lsO 

Mg/Ca 

Depth 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

0.5 

2.32 

0.5 

-0.671 

20 

2.24 

20 

-1.533 

40 

2.30 

40 

-1.883 

60 

2.58 

51 

-1.623 

80 

2.65 

60 

-1.881 

100 

2.55 

80 

-1.427 

121 

2.86 

100 

-1.484 

140 

2.87 

110 

-1.361 

144 

2.87 

121 

-1.171 

145 

2.76 

140 

-0.896 

151 

2.53 

144 

-1.234 

155 

2.45 

145 

-1.148 

160 

2.38 

151 

-0.910 

165 

2.56 

155 

-0.978 

170 

2.66 

160 

-1.353 

180 

2.80 

165 

-1.204 

200 

2.68 

168 

-1.045 

220 

2.55 

170 

-1.091 

250 

2.24 

180 

-0.472 

260 

2.99 

200 

-0.560 

290 

2.35 

220 

-0.632 

299 

2.46 

230 

-1.462 

320 

2.66 

240 

-0.123 

329 

2.29 

245 

-0.290 

340 

2.53 

250 

-0.888 

349 

2.02 

260 

-0.367 

360 

2.16 

265 

-0.064 

371 

2.33 

280 

-1.426 

170 


380 

2.05 

290 

-0.248 

401 

2.73 

299 

-0.333 

419 

2.55 

320 

-0.921 

430 

2.28 

329 

-0.180 

441 

2.55 

340 

-0.248 

450 

2.26 

349 

-0.519 

460 

2.48 

360 

-0.292 

371 

-1.058 

380 

-1.173 

401 

-0.348 

419 

-0.488 

430 

-0.323 

441 

-0.463 

450 

-0.582 

460 

-2.210 

BJ8-03 

136GGC 

G.  ruber 

Mg/Ca 

Depth 

5lsO 

Depth  (cm) 

(mmol/mol) 

(cm) 

%o  PDB 

2.5 

4.94 

2.5 

-2.886 

3.5 

4.66 

3.5 

-2.939 

6.5 

5.59 

8.5 

-2.798 

8.5 

4.54 

16.5 

-2.825 

10.5 

4.47 

24.5 

-2.777 

12.5 

5.26 

32.5 

-2.851 

14.5 

4.83 

40.5 

-2.898 

16.5 

4.49 

48.5 

-2.636 

18.5 

4.55 

56.5 

-3.059 

20.5 

4.67 

64.5 

-2.885 

22.5 

4.88 

72.5 

-3.217 

24.5 

4.70 

80.5 

-3.253 

26.5 

4.52 

88.5 

-3.243 

28.5 

4.65 

96.5 

-2.821 

32.5 

4.66 

104.5 

-3.105 

34.5 

4.70 

112.5 

-3.118 

36.5 

4.61 

120.5 

-2.921 

38.5 

4.32 

128.5 

-3.045 

40.5 

4.76 

136.5 

-2.958 

42.5 

4.62 

144.5 

-3.229 

44.5 

4.85 

152.5 

-3.132 

46.5 

4.60 

160.5 

-2.986 

48.5 

4.62 

168.5 

-3.283 

50.5 

4.76 

176.5 

-3.157 

56.5 

4.57 

184.5 

-2.762 

64.5 

4.43 

192.5 

-2.736 

66.5 

4.62 

200.5 

-2.404 

68.5 

4.67 

208.5 

-2.378 

70.5 

4.41 

216.5 

-2.174 

80.5 

4.49 

224.5 

-1.432 

171 

88.5 

4.76 

232.5 

-2.093 

96.5 

4.29 

240.5 

-1.712 

104.5 

4.65 

248.5 

-1.502 

112.5 

4.19 

256.5 

-1.363 

120.5 

4.68 

264.5 

-1.524 

128.5 

4.34 

272.5 

-1.461 

132.5 

4.45 

280.5 

-1.719 

134.5 

4.49 

288.5 

-1.590 

136.5 

3.95 

296.5 

-2.050 

138.5 

4.64 

304.5 

-1.744 

140.5 

4.50 

312.5 

-1.625 

142.5 

4.50 

320.5 

-1.736 

144.5 

4.74 

328.5 

-1.809 

146.5 

4.55 

336.5 

-1.694 

148.5 

4.35 

344.5 

-1.643 

150.5 

4.51 

352.5 

-1.490 

152.5 

4.53 

360.5 

-1.695 

154.5 

4.73 

368.5 

-1.917 

156.5 

4.30 

376.5 

-1.562 

158.5 

4.50 

384.5 

-1.442 

160.5 

4.38 

392.5 

-1.263 

162.5 

4.55 

400.5 

-1.991 

164.5 

4.44 

166.5 

4.40 

168.5 

4.34 

170.5 

4.21 

172.5 

4.87 

174.5 

4.62 

176.5 

4.26 

178.5 

4.86 

180.5 

4.81 

182.5 

4.49 

184.5 

4.98 

186.5 

4.62 

188.5 

4.72 

190.5 

4.61 

192.5 

4.36 

194.5 

4.21 

196.5 

4.42 

198.5 

4.08 

200.5 

4.62 

202.5 

4.20 

204.5 

4.47 

206.5 

4.14 

208.5 

4.33 

210.5 

4.44 

212.5 

4.63 

216.5 

4.21 

218.5 

4.30 
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220.5 

3.89 

222.5 

4.10 

226.5 

4.06 

228.5 

4.08 

230.5 

4.11 

232.5 

3.98 

234.5 

4.18 

236.5 

4.11 

238.5 

3.91 

240.5 

3.87 

242.5 

3.75 

244.5 

3.67 

246.5 

3.58 

250.5 

4.03 

252.5 

3.84 

256.5 

3.81 

264.5 

3.71 

272.5 

3.78 

280.5 

3.48 

286.5 

3.45 

296.5 

3.57 

304.5 

3.44 

312.5 

3.64 

320.5 

3.66 

328.5 

3.34 

336.5 

3.69 

344.5 

3.63 

352.5 

3.50 

360.5 

3.76 

368.5 

3.56 

376.5 

3.39 

384.5 

3.43 

392.5 

3.31 

400.5 

3.10 

BJ8-03 

P. 

136GGC 

obliquilocula  ta 

Mg/Ca 

Depth 

Depth  (cm) 

(mmoi/mol) 

(cm) 

2.5 

3.47 

3.5 

3.5 

2.98 

8.5 

14.5 

2.54 

16.5 

18.5 

2.61 

18.5 

20.5 

2.43 

20.5 

22.5 

2.52 

22.5 

34.5 

2.47 

24.5 

38.5 

2.70 

32.5 

40.5 

2.65 

34.5 

56.5 

3.07 

38.5 

64.5 

3.09 

40.5 

5lsO 
%o  PDB 

-2.286 

-2.313 

-2.443 

-2.591 

-2.592 

-2.315 

-2.247 

-2.093 

-2.429 

-2.563 

-2.032 
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66.5 

2.54 

48.5 

-1.947 

68.5 

2.61 

56.5 

-2.064 

70.5 

2.82 

64.5 

-2.065 

72.5 

2.66 

68.5 

-2.526 

80.5 

2.60 

70.5 

-2.113 

88.5 

2.70 

72.5 

-2.189 

96.5 

2.75 

80.5 

-1.827 

104.5 

2.85 

88.5 

-1.509 

136.5 

2.82 

96.5 

-1.805 

144.5 

3.05 

104.5 

-2.319 

152.5 

3.09 

112.5 

-2.190 

160.5 

3.29 

120.5 

-2.100 

168.5 

3.22 

128.5 

-1.818 

184.5 

2.73 

136.5 

-1.406 

188.5 

2.64 

144.5 

-1.684 

190.5 

2.97 

152.5 

-2.264 

194.5 

2.46 

160.5 

-1.501 

196.5 

2.30 

168.5 

-2.298 

198.5 

2.40 

176.5 

-2.063 

202.5 

3.09 

188.5 

-2.254 

204.5 

2.72 

190.5 

-1.916 

206.5 

2.58 

192.5 

-1.440 

208.5 

2.71 

196.5 

-1.353 

210.5 

2.71 

198.5 

-1.315 

212.5 

2.65 

200.5 

-1.925 

214.5 

2.71 

202.5 

-1.424 

216.5 

2.94 

204.5 

-1.420 

218.5 

2.51 

206.5 

-1.342 

220.5 

2.39 

208.5 

-1.586 

222.5 

2.57 

210.5 

-0.863 

224.5 

3.09 

212.5 

-1.727 

226.5 

2.05 

214.5 

-1.408 

228.5 

2.46 

216.5 

-0.998 

230.5 

2.45 

218.5 

-0.625 

232.5 

2.74 

220.5 

-0.533 

234.5 

2.27 

222.5 

-1.547 

236.5 

2.26 

224.5 

-0.731 

242.5 

2.54 

226.5 

-1.137 

244.5 

2.33 

228.5 

-1.715 

246.5 

2.34 

230.5 

-1.538 

248.5 

2.55 

234.5 

-1.026 

250.5 

2.50 

236.5 

-0.631 

252.5 

2.89 

240.5 

-0.230 

254.5 

2.43 

242.5 

-0.383 

258.5 

2.32 

244.5 

-0.057 

258.5 

2.13 

246.5 

-0.276 

260.5 

2.75 

248.5 

0.026 

262.5 

2.29 

250.5 

-0.417 

264.5 

2.28 

252.5 

0.029 
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266.5 

2.57 

254.5 

-0.395 

268.5 

2.47 

256.5 

0.204 

258.5 

-0.190 

260.5 

-0.162 

262.5 

0.235 

264.5 

0.127 

266.5 

-0.322 

268.5 

-0.019 

272.5 

-0.009 
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